AFWAL-TR-80-3141 
Part  I 


(  <f 


LIVE 


A 


r-l  '< 


INVESTIGATION  OF  HIGH-ANGLE-OF-ATTACK 
^  MANEUVER-LIMITING  FACTORS 
Part  I:  Analysis  and  Simulation 


DONALD  E.  JOHNSTON 
DAVID  G.  MITCHELL 
THOMAS  T,  MYERS 

SYSTEMS  TECHNOLOGY,  INC. 

13766  SOUTH  HAWTHORNE  BOULEVARD 
HAWTHORNE,  CALIFORNIA  90250 

DECEMBER  1980 


;  V‘ 


TECHNICAL  REPORT  AFWAL-TR-80-3141,  PART  I 
Final  Report  for  Period  May  1976  —  July  1980 


Approved  for  public  release;  distribution  unlimited. 


FLIGHT  DYNAMICS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


81 


21 


NOTICE 


Vhen  Government  drawings,  specifications ,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Govern¬ 
ment  pwourement  operation,  the  United,  States  Government  thereby  incurs 
no  responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
government  may  have  formulated,  furnished,  or  in  any  my  supplied  the 
said  drawings,  specifications ,  or  other  data,  is  not  to  be  regarded  by 
implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented,  invention  that  may  in  any  my  be 
related  thereto . 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs 
(ASD/PA)  and  is  relsasable  to  the  National  Technical  Information  Service 
( NTIS ),  At  NTIS,  it  will  be  available  to  the  general  public,  including 
foreign  nations . 

This  technical  report  has  been  reviewed  and  is  approved  for  publica¬ 
tion . 


MICHAEL  E.  BISE,  Project  Engineer 
Control  Dynamic 3  Branch 
Flight  Control  Division 


FOR  THE  COMMANDER 


ROBERT  C.  ETTINGERjfuo lone l,  USAF 
Chief ,  Flight  Control  Division 


R.  0.  ANDERSON,  Chief 
Control  Dynamics  Branch 
Flight  Control  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organis¬ 
ation  please  notify  AFVAL/FIGC,  V-PAFB,  OH  4S433  to  help  us  maintain  a 
current  mailing  list 

Copies  of  this  report  should  not  be  returned  unless  return  is  re¬ 
quired  by  security  considerations,  contractual  obligations,  or  notice  on 
a  specific  dooument. 

Al  R  FORCE/S6780/13  July  1981  -  4*0 


UNCLASSIFIED 


1)  TR-  M-Hm-PT-l 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whtn  D alt  ErUtrad) 


REPORT  DOCUMENTATION  PAGE 


I  T-  -REPORT  HUMBER 


5iPC  READ  INSTRUCTIONS 

_ _ BEFORE  COMPLETING  FORM 

2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  HUMBER 


* 1  AFWAI/TR-80-3 141,  Part  I  \fy  ^  AXC  ^  L  H 

^  4.  TITLE  (and  Sybtltta)  JJOUi  HMHtM 

y  I  ’^iVESTIGATION  OFJIIGHjANGLE-OF-ATTACK  pQ  T  Final  llep€?t 

3/  VJ^May  1976 

i  Fart  1.  Analysis  and  Siinlilation*  |*  r ,n r : nu  1 '/ V 'h 


^EUVER^jgg^l 
Part  I.  Analysis 


5luly  198S 


__  ^Cauthosto  j  ■  y  ~ 

rf/T\  |  Donald  E. /Johnston /David  G.l Mitchell/"' 
/(jj  andj  Thomas  T^/ Myers  •  y ' 

«:  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Systems  Technology,  Inc.  _/ 

13766  South  Hawthorne  Boulevard,/  /,  / 
Hawthorne,  California  90250 

II.  CONTRnl  I  INrt  OFFICE  NAME  AND  ADDRESS  ' 

AFWAL/FIGC 

Air  Force  Wright  Aeronautical  Laboratories 
Air  Force  Systems  Command 
Wright-Patterson  AFB.  OH  45433 _ 


I.  CONTRACT  OR  GRANT  NUMBER(i) 

/.^  F33615-76-C-3^72 

KCgj 

10.  PRO0P*.M  ELEMENT,  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 


PE  62201F 

~  U.  REPORT  OATE 

M  .. . 

...  <*.  UUttER  ClF.P.ACSS 

203 


<&)i  ax 


U.  MONITORINO  AGENCY  NAME  A  ADDRESUf//  dlitirmnt  from  Controlling  Otlic *)  IS.  SECURITY  CLASS,  ( ol  thlt  roport) 

fiL  \4  lL/tA  /  a/  >1  i  UNCLASSIFIED 


(It*) ltd*,  /  (ft)  4$  /  UNCLASSIFIED 

-  .  '  /  V  ^  IS*.  DECLASSIFY  ATI  ON?  DOWN  GRADING 

■— ' — I  .  schedule 

TT  DISTRIBUTION  STATEMENT  (at  ihfi  RAjiorl) 

Approved  for  public  release; , distribution  unlimited. 

•  \jj yyi  -  i'h-  j  l1  (  J  •  i  \  ' 

I7.~0!ST  RIBUt"|ON  ST  AT  EM  ENT  (of  onltfd  IfTflTocVjO^  Udlfhtonllrom  Rtport) 


MB.  SUPPLEMENTARY  NOTES 


IB,  KEY  WORDS  (Con/lmi«  on  r«v«ra«  «fd«  // n«c»«««ry  one/  Identity  by  bfoc*  numb#0 

High  AOA,  flying  qualities,  stall/departuxc,  simulation 


BO.  \fl^TRACT  (Conf/nu*  on  rovf  «ld»  1/  n«c»««<ry  and  Idonttly  by  blocft  nutnbor)  —  - 

The  high-nngle-of-attack,  low-speed  stall/departure  characteristics  of  the 
F-4J  and  F-1.4A  are  analyzed,  using  a  six-degree-of-f reedotn  mathematical  model 
with  nonlinear  aerodynamics.  Cause-effect  relationships  are  investigated  for 
maneuver  limiting  factors  including  wing  rock,  nose  slice,  and  rolling  depar¬ 
tures.  CrosB-derlvatives  of  N^,  and  MB  alter  key  transfer  function, 

(continued) 


DD  1  JAN  ^71  1  473  EDITION  OF  I  NOY  «8  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OK  THIS  PAOE  fl«,.n  Dai'a  gnTlffll 


r-.uii 'tiVU Xli, ,i  j.  *  X*b. 


UNCLASSIFII D 


SECURITY  CLASSIFICATION  Of  THIS  PAGEfKTi»n  P«U  EnlinilJ 

^parameters.  A  piloted  simulation  validates  analytic  predictions  and  demon¬ 
strates  that  departure  warning,  susceptibility,  and  severity  are  strongly 
influenced  by  the  static  cross-derivatives.  A  connection  between  roll  numera¬ 
tor  parameter  values  and  pilot  perception  of  departure  susceptibility  and 
severity  is  identified.  Potential  modifications  for  the  high  AOA  sections  of 
the  MIL-F-8785B  Flying  Qualities  Specification  are  proposed:  a  criterion  for 
the  real  part  of  the  roll  numerator  root,  further  recommendations  for  minimiz¬ 
ing  departure  susceptibility  and  certain  sideslip  influences,  and  a  flying 
quality  rating  form  for  assessing  departure  and  recovery  characteristics.  ^ 

^Part  I,  Analysis  and  Simulation,  presents  a  summary  of  the  complete  in¬ 
vestigation  and  results.  Part  II,  Piloted  Simulation  Assessment  of  Bihrle 
Departure  Criteria,  presents  a  detailed  comparison  of  analytical  prediction 
and  piloted  simulation  results  for  a  specific  set  of  programmed  control  de¬ 
flections.  Part  III,  Appendices - Aerodynamic  Models,  contains  the  detailed 

aerodynamic  models  employed  in  the  F-4J  and  F-14A  high-angle-of-attack  analy¬ 
sis  and  validation  and  the  equations  of  motion,  aerodynamic  models,  control 
system  configurations,  etc.,  employed  in  the  piloted  simulation. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  of  TU,f  PAOEfHTi.n  D«l»  Knttrtil) 


FOREWORD 


4 


\ 


d 


1 


This  research  was  sponsored  by  the  Air  Force  Wright  Aeronautical 
Laboratories,  Air  Force  Systems  Command,  under  Contract  F33615-76-C- 
3072,  Work  Unit  24030514,  Mr.  Gary  K.  Hellmann  was  the  initial  project 
monitor.  This  responsibility  was  later  transferred  to  Mr.  Michael  E. 
Bise  (AFWAL/FIGC).  Support  of  the  piloted  simulation  was  also  provided 
by  the  Naval  Air  Development  Center  where  Mr.  Mark  Stifel  served  as 
project  monitor.  The  analytic  work  was  performed  at  Systems  Technology, 
Inc.,  Hawthorne,  California.  The  work  was  performed  during  the  period 
13  May  1976  through  30  July  1980.  The  STI  Technical  Director  was 
Mr.  I.  L.  Aahkenas.  Mr.  D.  E.  Johnston  was  Principal  Investigator  and 
STI  Project  Engineer.  The  piloted  simulation  was  accomplished  at  the 
McDonnell  Aircraft  Co.,  St.  Louis,  Missouri.  The  report  manuscript  was 
submitted  in  September  1980. 

The  Authors  wish  to  express  acknowledgment  and  thanks  to  their  many 
coworkers  for  contributions,  both  general  and  detailed,  in  the  program; 
at  STI,  Mr.  G.  L.  Teper  for  invaluable  aid  in  accomplishing  the  digital 
simulation,  Mr.  T.  T.  Myers  for  development  and  validation  of  the  F-4 
and  F-14  aerodynamic  models  and  for  most  of  the  initial  analytic  sup¬ 
port,  and  Mr.  R.  H.  Hoh  for  checkout  and  accomplishment  of  the  piloted 
simulation.  Mr.  Hoh  also  served  as  one  of  the  subject  pilots.  At 
MCAIR,  Mr.  H.  Passmore  directed  setup  and  operation  of  the  piloted 
simulation.  Special  thanks  are  due  to  Lt.  Col.  R.  M.  Cooper,  Maj.  J.  A. 
Fain,  Jr.,  and  Maj.  J.  Jannarone  of  the  6510th  Test  Wing  and  Maj.  P. 
Tackabury  of  the  Test  Pilot  School,  Edwards  Air  Force  Base,  for  their 
contribution  in  refining  the  high  angle  of  attack  flying  quality  rating 
scale  and  their  professional  approach  in  accomplishing  the  sometimes 
tedious  simulation  experiments. 

Finally,  appreciation  is  extended  to  Mr.  P.  Kelly  of  Grumman  Air¬ 
craft  Co.  and  Mr.  M.  Humphreys  of  the  Naval  Air  Test  Center  for  support 
in  obtaining  much  of  the  F-14  data  and  information,  to  Mr.  R.  Wood  of 
the  Air  Force  Flight  Test  Center  for  invaluable  comments  and  suggestions 
concerning  the  flying  quality  rating  scale,  and  to  Mr.  R.  Woodcock 
(AFWAL/FIGC)  for  his  careful  critique  and  editorial  refinement  of  this 
final  report. 


Aooesslon  Foi 

—NtTs  GRAltl 
T3TIC  TAB 
Unannounced 
Justiflcatio 

r  y 

« 

a 

n 

Hv  _  — I 

I  Distribution 
1  Aral lab tilt 
I  AvuL'i  '( 

Hi  1st  Spec 

\n 

/  . . 

y  Cortes 
ind/or 

lei 

L_  . 

Y 


iii 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . . . 

A.  Background . . . . 

B.  Technical  Approach. . 

II.  AIRCRAFT  HIGH  AQA  DYNAMIC  CHARACTERISTICS 

A.  F-4J  Analysis . . . 

B.  F-14A  Analysis . . 

III.  AERODYNAMIC  CROSS-COUPLING  INFLUENCES.... 

A.  Characteristic  Modes . . . . 

B.  n|  Numerator . 

C.  N§  .  Numerator . . 

°stab 

D.  Summary . . . 


Page 

1 

1 

3 

5 

5 

42 

66 

66 

70 

80 

87 


IV.  CONFIGURATIONS  SELECTED  FOR  SIMULATION .  90 

A.  Aerodynamic  Configurations. . 91 

B<  Augmented  Configurations .  103 

C.  Summary  of  Configurations  and  Predicted 

High  AOA  Characteristics . 106 

V.  PILOTED  SIMULATION .  108 

A.  Simulator.......... . 109 

B.  Head-Up  Display .  109 

C.  Flight  Control  System . Ill 

D.  Equations  of  Motion. . Ill 

E.  Aerodynamic  Models . 113 

F.  Data  Recorded . 119 

G.  Tasks  and  Maneuvers . . . 119 

H.  Flying  Quality  Rating  Scales..... .  121 

VI.  HIGH  AOA  FLYING  QUALITIES  ASSESSMENTS .  127 

A.  Ratings  Scale  and  Questionnaire  Assessment..............  129 

B.  Configuration  Assessments. . .  135 

C.  Summary .  152 


v 


IrfECEDlNQ  PiflB  BLANK -NCT  FXitm 


ZABLE  OF  OOTiaHTS  (Concluded) 


Page 

VII.  DEPARTURE  SUSCEPTIBILITY  CRITERION....,  . .  155 

A.  Closed-Loop  Analysis. ...... »  . .  156 

B.  Possible  Departure  Susceptibility  Criteria . .  165 

C.  Summary . . . .  172 

VIII.  CONCLUSIONS  AND  RECOMMENDATIONS . 173 

A.  Key  Maneuver-Limiting  Factors . . . .  173 

B.  Additional  Findings  and  Conclusions . . . .  175 

C.  Proposed  MIL-F-8785C  Modifications . . .  177 

REFERENCES . 180 


vi 


•% 


LIST  OF  ILLUSTRATIONS 

Pajje 


1.  Basic  F-4  Trimmed  Lift  Curve.... .  6 

2.  F-4J  Flight  Control  System......... . .  8 

3.  F-4E  Flight  Test  Pull-Up  and  Steady  Wing  Rock... . . .  11 

4.  F-4J  Model  Wing  Rock  at  a  *  23  deg,  0  ■  0; 

Pitch  and  Yaw  SAS  On . 7 . 7 .  12 

5.  Departure  During  Left  Turn. . . . . .  13 

6.  F-4J  Pitch  Coefficient  Variation  with  a,  0 . 15 

7.  F-4J  Roll  Coefficient  Variation  with  a,  0. . 15 

8.  F-4J  Yaw  Coefficient  Variation  with  a,  0. . . 16 

9.  F-4J  Lateral-Longitudinal  Short-Period  Root 

Migration  with  a;  0O  ■  0 . . . .  17 

10.  F-4E  Stall  Approach  Character  ’  sties ....  . . . .  18 

11.  F-4J  6  DOF  Linearized  Equation;  Lateral-Longitudinal 

Root  Migration  with  a  and  0 . . . .  19 

12.  Open-Loop  Nonlinear  Airframe  Response  to 

Pulse  Aileron,  aQ  -  21  deg,  0g  ■  0  deg . .  20 

13.  Matrix  for  Coupled  Equations  with  0Q  j*  0 . . .  21 

14.  F-4J  Eigenvectors;  aQ  ■  21  deg,  0Q  ■  1.5  deg . 23 

15.  F-4J  Vector  Polygons;  aQ  ■  21  deg,  0Q  ■  1.5  deg .  24 

16.  Open-Loop  Wing  Rock  with  Cmn  Small; 

aQ  -  21  deg,  0O  -  0  deg....? .  26 

17.  Root  Migration  with  AOA;  F-4J . . . .  29 

18.  Root  Migration  with  AOA;  F-4J .  32 

A 

19.  Effect  of  Lateral-Longitudinal  Coupling  on  Ng  /A 

Roots;  aQ  *  23  deg,  0Q  ■  5.5  deg . .  34 

20.  System  Survey  for  0  +  <5stab  Closure;  a0  ■  23  deg, 

0O  -  5.5  deg .  35 


vii 


LIST  OF  ILLUSTRATIONS  (Continued) 

0 

N?  First-Order  Root  Migrations  with  AOA;  F-4J . . . 
3tdD 

F-4J  Static  Aerodynamic  Cross-Coupling 
Derivatives;  30  -  5.5  deg . . . . . 

F-14A  Roll  Coefficient  Variation  with  a,  3 . 

F-14A  Yaw  Coefficient  Variation  with  a,  3... . 

F-14A  Pitch  Coefficient  Variation  with  a,  3 . 

F-14A  Longitudinal  Control  System. . . . . . 

F-14A  Lateral  Control  System . . . . 

F-14A  Directional  Control  System . . 

F-14A  Longitudinal  Static  Trim  Comparison . . . 

F-14A  Flight  234-1  Traces..... . 

F-14A  Model  Match  of  Flight  234-1 . 

F-14A  Flight  243  Traces . 

F-14A  Flight  236-1  Pull-Up . 

F-14A  Model  Simulation  of  Flight  236-1  Pull-Up . 

F-14A  Root  Migrations  with  AOA. . 

Ng^yn  Variation  with  AOA;  F-14A . . . 

F-14A  Eigenvectors;  otQ  -  20  deg,  30  “  1.5  deg . 

F-14A  Vector  Polygons  for  Dutch  Roll  Mode; 

oQ  -  20  deg,  30  ■  1*5  deg . . . 

N$d  Root  Migration  with  AOA;  F-14A.. . . 

N®stab  First-Order  Root  Migration  with  AOA;  F-14A.. 

F-14A  Static  Aerodynamic  Cross-Coupling 
Derivatives . . . . 

Migration  of  Roots  of  K( 6o)^coupl  Increasing  M 


I 


LIST  OF  ILLUSTRATIONS  (Continued) 


43.  Migration  of  Roots  of  Uncoupled  Characteristic 
Modes  as  Cross-Coupling  Increases.............. 


44.  n|  /A  Pole-Zero  Locations  Coupled  Versus  Uncoupled; 
F-W,  a  -  21  deg,  SD  *  5.5  deg . . 


45.  Closure  of  p/s  ■+•  to  Obtain  Time  Vectors, 

F-4J,  aQ  “  21  deg,  $Q  •  5.5  deg . 


46.  Vector  Polygons  for  Closed-Loop  tugj.  (Approximating  uu) 
F-4J;  aQ  «  21  deg,  •  5.5  deg . . . 


47.  Matrix  for  Approximation  of  Coupled  Numerator  Njfc  ^ . . .  75 

48.  System  Survey  of  Effect  of«£a  on  ■  0 . .  77 

49.  System  Survey  for  Effect. of  N«  on  NI 

(with  X'a  -  3.4251  1/sec2). . . .!$* .  78 

50.  Time  Vectors  for  Closed-Loop  Mode  1/Ta 

(Approximating  l/Tg^)  . . . . . .  81 


51.  Matrix  for  Approximation  of  Coupled  ^6atab  *•••••• 

52.  System  Survey  for  Effect  of  on  N5Bta^>  Na  ■  0. 

+  0 

53.  System  Survey  for  Effect. of  NI  on  N5 

(with  -  3.4251  1/sec2).... . . 


54.  System  Survey  for  Effect  of  Combined  Cross-Coupling 

“  "‘.Mb  . 


55.  Key  Parameters  for  Configuration  A. 


56.  Change  in  Cj  for  Configuration  B.. 

57.  Key  Parameters  for  Configuration  B. 


58.  Change  in  Cjg  for  Configuration  C. 


59.  Key  Parameters  for  Configuration  C. 


60.  Change  in  Cn^  for  Configuration  D., 

61.  Key  Parameters  for  Configuration  D. 


62.  Cn^  for  Configurations  Simulated . . . .  101 


ix 


LIST  OF  ILLUSTRATIONS  (Continued) 


Page 

63.  Weissman  Criterion  Predictions  for  the 

Simulation  Configurations . 102 

64.  Key  Parameters  for  Configuration  A2 . 104 

65.  Key  Parameters  for  Configuration  C2> . 105 

66.  Ur  Versus  Cdud  f°r  Augmented  and  Unaugmented 

Configurations  A  and  C.... . . . . .  106 

67.  Manned  Air  Combat  Simulator  I... . 110 

68.  Head-Up  Display  Used  in  Simulation....................... .  Ill 

69.  Lateral-Directional  Control  System....... . 112 

70.  Coefficients  as  Extended  for  |a|  <  180  deg 

and  1 3 1  <  90  deg......... . . . . . . . .  114 

71.  Cooper-Harper  Handling  Qualities  Rating  Scale .  122 

72.  Tentative  Departure  Rating  Scale...... . 124 

73.  Modified  Departure  Rating  Scale  (Mod  II) . 131 

74.  Rating  Comparisons  Between  Pilots. . 133 

75.  Configuration  A.  Departure  and  PSG  from  Wings  Level 

Tracking  Task;  Pilot  RC..... . . . .  139 

76.  Configuration  A?  Departure  and  PSG  from  Wings  Level 

Tracking  Task;  Pilot  RC . . .  149 

77.  Typical  F-4E  Roll  Departure/Spin  Entry . 150 

78.  Assumed  Pilot  Lateral  Loop  Closure;  $  *  6atk . 156 

79.  Survey  Plot  for  Roll  Loop  Closure;  aQ  ■  16  deg, 

0O  -  0 .  157 

80.  Survey  Plot  for  Roll  Loop  Closure;  aQ  -  23  deg, 

B0  -  0 .  159 

81.  Maximum  Pilot  Gain  for  Stable  $  +  ^stk  Closure; 

Configuration  A^  . . . 160 

82.  Roll  Loop  Closure  Gains  Employed  with 

Configuration  A^ . 160 


x 


LIST  OP  ILLUSTRATIONS  (Concluded) 


83.  Roll  Loop  Closure  Gains  Employed  with 

Configuration  A2 . . . . . . . .  162 

84.  Roll  Loop  Closure  Gains  Employed  with 

Configuration  C^ . . . .  162 

85.  Roll  Loop  Closure  Gains  Employed  with 

Configuration  C£ . . . . . .  163 

86.  Roll  Loop  Closure  Gains  Employed  with  Configuration  D .  163 

87.  Roll  Loop  Closure  Gains  Employed  with  Configuration  B .  164 

88.  Departure  Susceptibility  Rating  Versus  Lateral 

Closed-Loop  Divergence  Potential. . . . . .  167 

89.  Possible  Modifications  to  Welssman  Criterion .  168 

90.  Variation  of  Pilot  Ratings  with  Time  to  Double 
Amplitude  of  Aperiodic  Motion,  Landing 

Approach  Condition .  170 

91.  Variation  of  Pilot  Ratings  with  Equivalent  Aperiodic 

Inverse  Time  Constants . .  170 

92.  Migration  of  Closed-Loop  Pole  into  RHP  with 

Typical  Roll  Loop  Gains.. . . . . . . .  171 


xi 


LIST  OF  TABLES 


1.  Comparison  of  F-4E  and  F-4J  Veighta  an  tnertias 

2.  Comparison  of  F-14A  and  F-4J  Weights,  Inertias, 

end  Geometry . . . . . 


3.  Coupling  Numerator  Factors.. ... 

4.  N$  Numerator  Approximations, 


5.  Key  Maneuver-Limiting  Parameters. 


6.  Configuration  Matrix. 


7.  Los's  of  Control/Departure/Recovtry  debriefing  Guide. 


8.  Run  Matrix  for  Piloted  Simulation. 


9.  Number  of  Runs 


10.  Configuration  A]  Maneuver-Limiting  Factors . 

11.  Configuration  B  Maneuver-Limiting  Factors . 

12.  Configuration  Cj  Maneuver-Limiting  Factors . 

13.  Configuration  D  Maneuver-Limiting  Factors . 

14.  Configuration  A2  Maneuver-Limiting  Factors... . 

15.  Configuration  C2  Maneuver-Limiting  Factors. . 

16.  Comparison  of  Pilot  Ratings,  Initial  Versus  Reruns 

17.  Comparison  of  Predicted  and  Actual  Departure 

Susceptibility  Assessments . . . 


Page 

10 

46 

68 

80 

90 

107 

126 

128 

129 

137 

140 

142 

144 

147 

151 

154 

165 


xil 


DEFINITION  OF  AXIS  SYSTEMS 


Unless  otherwise  Indicated,  all  moment  coefficients  and  derivatives 
are  referenced  to  the  fuselage-centerline-oriented  body  system  of  axes. 
All  longitudinal  force  coefficients  and  derivatives  are  referenced  to 
the  flight  path  (or  wind)  system  of  axes.  Thene  axis  systems  are  de¬ 
fined  in  the  sketch  below. 
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unoDDcnoK 

k.  BACKGROUND 

For  modern  military  aircraft,  manageable  flight  at  high,  near-stall 
angles  of  attack  is  a  key  effectiveness  and  safety  factor.  The  safe 
angles  of  attack  and  sideslip  attainable  in  level  and  maneuvering  flight 
often  represent  the  important  differences  in  both  offensive/defensive 
capability  and  accident  histories  of  otherwise  equivalent  aircraft.  In 
the  past  generation  of  fighter  aircraft,  surprise  loss  of  control  depar¬ 
ture  generally  has  been  of  concern  as  a  key  maneuver-limiting  factor. 
The  more  recent  generation  has  emphsized  design  for  departure  resistance 
but  not  necessarily  prevention.  Many  operational  pilots  desire  the 
capability  to  depart  the  aircraft  on  command  as  a  last-ditch  defensive 
maneuver,  since  departure  is  generally  violent  with  unpredictable 
results  from  the  attacker's  viewpoint.  The  key  element  then  revolves 
around  departure  resistance  (or  susceptibility),  warning,  severity,  and 
recoverability.  Accordingly,  the  issue  of  improved  utility  and  safety 
at  high  angles  of  attack  has  been  re-recognized  as  an  area  of  fruitful 
research  and  development. 

The  recent  upsurge  in  interest  has  resulted  in  numerous  studies 
devoted  to  identifying  the  safety-related  maneuver-limiting  phenomena 
Inherent  in: 

•  Buffet 

•  Pitchup  or  "dig-in" 

•  Roll  reversal 

•  Wing  rock 

•  Nose  slice 

•  Rolling  departure  —  snap  roll 


1 


Some  of  these  phenomena  are  perturbations  about  a  steady  flight 
condition  (e.g.,  buffet,  wing  rock),  which  not  only  limit  tracking 
capability  (as  shown  in  Reference  1)  but  also  may  serve  as  warning  of 
approaching  “departure"  from  controlled  flight  (e.g.,  nose  slice,  roll¬ 
ing  departure).  The  more  interesting  and  mysterious  of  these  phenomena 
are  the  latter,  and  considerable  effort  has  been  expended  in  correlating 
the  angles  of  attack  at  which  such  behavior  occurs  with  various  aero¬ 
dynamic  stability  and  control  parameters  (e.g.,  References  2-4).  Such 
open-loop  "correlations,"  while  useful,  are  not  so  satisfying  or  in¬ 
structive  as  the  more  positive  identification  of  causal  factors  offered 
by  the  methods  developed  in  Reference  5. 

The  difficulty  is  in  identifying  causal  relationships  in  the  complex 
Interactive  pilot/vehicle  situation  during  uncoordinated  flight,  where 
aerodynamic  cross-coupling  can  compound  vehicle  dynamic  characteris¬ 
tics.  The  previous  (Reference  S)  analysis  and  pilot  simulation  showed 
the  static  aerodynamic  cross-coupling  coefficients  and  N£  to  contri¬ 
bute  to  closed-loop  nose-slice  departure  susceptibility  and  severity  in 
the  A-7  aircraft.  It  represented  an  initial  successful  and  promising 
new  attack  on  an  old  problem.  But  configuration  differences  among  air¬ 
craft  types,  models  and  even  loadings  have  been  observed  to  produce 
gross  differences  in  behavior  through  changed  aerodynamics  and  inertial 
characteristics.  Thus  it  is  desirable  that  the  same  methods  be  applied 
to  additional  high-performance  fighter  aircraft  having  widely  differing 
high  AOA  handling  characteristics  to  see  if  further  cause-effect  rela¬ 
tionships  can  be  identified. 

Accordingly,  the  stated  or  implied  goals  of  this  program  are  to: 

•  Identify  key  design  parameters  that  limit  high- 
angle-of-attack  maneuverability  for  contemporary 
high-performance  attack  and  fighter-type  aircraft. 

•  Postulate  fundamental  aerodynamic  and  control 
system  design  methodologies  that  will  alleviate 
the  limiting  conditions. 

•  Formulate  handling  qualities  requirements  for  high 
AOA  maneuvering  flight  to  be  incorporated  in  MIL- 
F-8785C,  the  military  flying  qualities  specifica¬ 
tion  for  piloted  airplanes. 
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B.  TECHNICAL  APPROACH 


The  research  encompassed  four  major  technical  areas.  The  first  was 
devoted  to  development  and  validation  of  aerodynamic  models  for  two 
fighter  aircraft  having  significantly  different  high  AOA  maneuver- 
limiting  characteristics.  The  two  aircraft  selected,  the  F-4J  and  the 
F-14A,  were  found  to  represent  almost  opposite  extremes  in  their  range 
of  departure  susceptibility  and  severity.  Attention  was  then  turned  to 
investigation  of  maneuver-limiting  characteristics  of  each  airframe  and 
identification  of  causal  relationships  based  on  quaai-linear*  analysis 
at  symmetric  and  asymmetric  ((3  i4  0)  flight  conditions.  The  causal 
parameters  were  varied  to  alter  the  dynamic  characteristics  of  each 
airframe  at  high  AOA.  Predicted  characteristics  were  checked  utilizing 
complete  slx-degree-of-f reedom  <6  DOF)  models  with  nonlinear 
aerodynamics,  and  any  differences  between  the  quasl-linear  frozen  point 
analytic  predictions  and  the  nonlinear  model  results  were  resolved. 

The  second  task  Involved  development  of  methods,  criteria,  and  an 
associated  pilot  rating  scale,  for  evaluation  of  handling  qualities  when 
approaching  controllability  limits.  Previous  analysis  and  simulation 
(Reference  5)  has  demonstrated  that  the  Cooper-Harper  handling  quality 
rating  scale  is  inapplicable  for  departure/recovery  flight  situations. 
Therefore,  a  new  pilot  rating  scale  is  required  specifically  for  high 
angles  of  attack. 

The  third  task,  pilot  simulation,  was  performed  to  evaluate  the 
influence  of  intentionally  varied  maneuver-limiting  characteristics 
(wing  rock,  nose  Bllce,  etc.)  under  different  normal  flying  situations 
such  as  training  flights  and  air  combat  tracking  tasks;  to  identify  key 


A 

It  is  recognized  that  these  linear  analysis  techniques  are  applied 
to  phenomena  which  may  be  nonlinear  in  nature.  However,  the  intent  of 
the  study  ia  to  investigate  in  an  understandable  manner  the  conditions 
which  can  lead  to  or  precipitate  rapid  changes  in  motion,  as  opposed  to 
analyzing  the  fully  developed  large  angle  phenomena.  Therefore,  piece- 
wise  linearity  is  appropriate  —  but  care  should  be  taken  in  attempting 
to  extrapolate  the  results. 
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flying  quality  parameters  in  terms  of  departure  resistance,  warning, 
severity,  and  recovery  for  possible  Inclusion  in  the  flying  qualities 
specification;  and  to  exercise  and  refine  as  necessary  the  new  pilot 
rating  scale . 

The  fourth  task  involved  assessment  and  correlation  of  results  of 
the  simulation  with  other  applicable  results  and  data;  and  formulation 
of  high  AOA  maneuvering  flight  generalized  design  guides  and  flying 
qualities  criteria  for  incorporation  into  MIL-F-8785C. 
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AIRCRAFT  HIGH  AOA  DTMAHIC  CHASACTSRISTICS 

Thla  section  presents  condensed  results  of  an  extensive  analysis  of 
the  high  AOA  dynamic  characteristics  of  the  two  example  aircraft.  The 
F-4J  is  treated  first,  followed  by  the  F-14A.  The  high  AOA  character¬ 
istics  of  each  aircraft  are  described  as  determined  in  flight  test  and 
operational  deployment.  The  analytical  models  are  briefly  summarized 
and  example  comparisons  shown  to  demonstrate  the  match  obtained  with  the 
actual  aircraft. 

Results  of  linear,  frozen-point  dynamic  analysis  at  a  series  of 
symmetric  and  asymmetric  (8  4  0)  flight  conditions  are  summarized.  This 
analysis  provides  insight  to  potential  causal  factors  behind  undesirable 
characteristics  at  high  AOA,  showing  examples  of  advantages  and  limita¬ 
tions  of  such  analysis  whan  applied  to  situations  in  which  aerodynamics 
are  obviously  quite  nonlinear.  In  particular,  this  section  Identified 
lateral-directional  static  aerodynamic  cross-coupling  as  a  key  factor  in 
determining  high  AOA  stability  and  controllability  characteristics. 

A.  F-4J  ANALYSIS 

The  F-4J  model  of  the  F-4  family  was  selected  because  a  previously 
validated  aerodynamic  model  (Reference  6)  was  in  hand.  However,  all 
versions  of  the  aircraft  with  an  unslattad  leading  edge  wing  are 
reported  to  exhibit  quite  similar  high  AOA  flying  characteristics: 
increasing  AOA  produces  buffet  of  increasing  severity,  increasingly 
adverse  yaw,  wing  rock,  mild  pitch-up,  stall,  and  finally  departure.  A 
representative  trimmed  lift  curve  with  onset  of  these  handling  phenomena 
is  shown  in  Figure  1  from  Reference  7.  Above  about  11  deg  AOA,  aileron 
adverse  yaw  becomes  significant  and  rudder  is  used  for  roll  control. 

Wing  rock  appears  primarily  as  a  rolling-sideslipping  motion.  It 
is  often  described  by  pilots  as  a  divergent  dutch  roll  which  is  easily 
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MODEL  F-4  AIRCRAFT 


*  1  tilM 


Figure  1.  Basic  F-4  Trimmed  Lift  Curve  (from  Ref.  7) 


aggravated  or  caused  by  pilot  lateral  stick  inputs.  The  angle  of  attack 
for  onset  and  the  severity  actually  depend  upon  several  facto..  , ,  includ¬ 
ing  c.g.  location,  roll  inertia  or  loading  (clean  vs.  external  stores), 
gear  and  flap  settings,  etc.  It  is  generally  more  pronounced  at  light 
weight,  with  low  roll  Inertia,  and  at  forward  c.g.  The  oscillations 
generally  require  10  sec  or  more  of  sustained  high  AOA  to  develop. 

Amplitudes  can  be  large  enough  (±30  deg  $  and  ±10  deg  ?)  to  interfere 
with  tracking  performance,  or  to  be  highly  uncomfortable  when  escorting 
considerably  slower  aircraft. 

The  mild  pitch-up  is  alleviated  by  a  pitch  rate  damper  but  does 

result  in  some  lightening  of  stick  forces.  It  has  been  described  by 
some  pilots  as  a  "dig-in"  during  slow  turns.  Possibly  it  is  more  of  a 
stall  warning  cue  than  a  maneuver-limiting  problem. 

Departure  is  described  variously  as  nose  slice,  roll  departure,  or 
both.  Its  onset  also  appears  to  depend  upon  c.g.  location  and  aircraft 

inertia  variations.  Although  shown  in  Figure  1  to  occur  at  considerably 

higher  AOA  than  wing  rock,  departure  is  not  necessarily  preceded  by  wing 
rock  —  especially  if  AOA  is  rapidly  increased. 
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1.  Analytic  Model 


The  F-4J  aerodynamic  model  (Reference  6)  was  derived  from  three 
excellent  but  separate  data  sources  (References  8-10).  It  was  assembled 
to  model  as  closely  as  possible  the  F-4J  sub-  and  transonlcally  for  a 
moving-base  simulation  used  to  train  U.S.  Navy  pilots  In  air  combat 
maneuvering.  It  therefore  had  to  be  a  very  good  representation  of  the 
complete  aircraft  in  high  AOA  maneuvering  up  through  stall.  The  instruc¬ 
tor  pilots  pronounced  the  simulation  "credible  and  realistic  of  actual 
F-4J  handling  characteristics." 

The  history  of  development  of  the  complete  aerodynamic  model  and  the 
data  package  is  presented  in  Part  III,  Appendix  I.  In  brief,  the  aero¬ 
dynamic  coefficients  represent  the  following  flight  regime: 

0  <  o  <  +110  deg 
-30  <  8  <  +30  deg 

For  this  analysis  the  look-up  table  data  were  restricted  to  low  Mach 
(<0.4)  and  s  single  altitude  (h  ■  15,000).  Weight  and  inertia  charac¬ 
teristics  representative  of  partially  full  internal  fuel  tanks,  an  empty 
600  gal  centerline  tank,  and  missile  pylons  on  wing  stations  2  and  8 
have  been  assumed. 

Conventional  6  DOF  equations  plus  kinematic  terms  are  employed  with 
the  moment  equations  in  body  centerline  axes,  the  force  equations  in 
wind  axes,  and  aircraft  orientation  angles  in  standard  Euler  axes.  The 
nonlinear  equations  of  motion,  aerodynamic  force  and  moment  equations, 
auxiliary  equations,  etc.,  are  detailed  in  Appendix  I  (Tart  III). 

A  block  diagram  of  the  basic  flight  control  and  augmentation  mech¬ 
anization  is  presented  in  Figure  2.  Aircraft  control  is  exerted  through 
horizontal  sfcabilator  (58tab)»  rudder  (6r),  aileron  (6a),  and  spoiler 

(6 _ )  surfaces.  The  aileron  and  spoiler  systems  are  Interconnected 

■P 

such  that  lateral  deflection  of  tha  control  stick  produces  downward 
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deflection  of  one  aileron  and  upward  deflection  of  the  opposite  spoiler. 
Stability  augmentation  is  provided  in  all  three  axes,  but  it  is  common 
practice  for  the  pilot  to  turn  the  roll  damper  off  during  high  AOA 
maneuvering.  The  complex  artificial  feel  system  was  not  modeled  in  this 
analysis  and  simulation. 

2.  Model  Validation 

As  a  part  of  the  checkout  and  acceptance  test,  the  simulation  based 
on  the  Reference  6  aerodynamic  data  was  flown  through  various  offensive 
and  defensive  air  combat  maneuvers,  stalls,  and  departures  by  Navy  in¬ 
structor  pilots,  who  indicated  that  it  adequately  represented  the  F-4J 
handling  and  performance.  This  provided  the  first  gross  validation  of 
the  aerodynamic  model. 

Actual  flight  test  traces  of  F-4J  high  AOA  dynamic  responses  are  not 
available  to  compare  against  those  of  the  mathematical  models.  The 
available  high  AOA  flight  traces  are  from  an  F-4E  stall/post-stall 
flight  test  (Reference  11).  However,  Reference  12  indicates  that  all 
hard-wing  models  of  the  F-4  have  approximately  the  same  stall/departure 
characteristics.  The  F-4E  has  a  nose  section  approximately  5  ft  longer 
than  the  F-4J;  this  could  cause  some  difference  in  side  force  due  to 
asymmetric  vortex  shedding  at  high  AOA  (Reference  13).  In  addition  the 
flight  test  vehicle  was  equipped  with  a  spin  chute,  had  a  reinforced  aft 
fuselage  structure  to  handle  the  spin  chute  loads,  and  is  presumed  to 
have  offsetting  ballast  in  the  nose.  Thus,  its  pitch  and  yaw  inertia 
characteristics  are  quite  different  from  the  F-4J,  as  shown  in  Table  1. 
This  difference  is  assumed  to  have  minor  Influence  on  the  sequence  of 
stall/departure  characteristics  with  increasing  AOA.  The  larger  yaw 
inertia  and  longer  nose  (greater  vortex  shedding)  of  the  F-4E  could  be 
offsetting  departure  susceptibility  factors.  But  the  inertia  differ¬ 
ences  should  produce  measurable  differences  in  dynamic  response  parame¬ 
ters  such  as  dutch  roll  frequency,  response  to  control  Inputs,  and 
angular  rates  associated  with  departure.  For  example,  scaling  parame¬ 
ters  by  the  ratio  of  inertias  gives  ratios  of  uigp  and  of  1.11. 
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TABLE 

1.  COMPARISON  OF  F-4E  AND 
WEIGHTS  AND  INERTIAS 

F-4J 

F-4E  (Reference  11) 

F-4J 

<£.  stores 

Empty  tank 

Empty  tank 

Wing  stores 

Pylons  1,  2,  8,  9 

Pylons  2,  I 

W  (lb) 

40,000 

37,000 

c.g.  (X  c) 

28.1 

29.3 

Ix  (slug-ft2) 

27,500 

23,850 

Iy  (8lUg-ft2) 

157,000 

127,400 

Iz  (slug-ft2) 

180,600 

146,000 

Ixz  (slug-ft2)  5,500 

2,210 

Example  comparisons  between  the  F-4E  test  results  and  the  nonlinear 
6  DOF  F-4J  math  model  are  shown  In  Figures  3-5.  Figure  3  shows  the  F-4E 
at  a  nearly  steady  AOA  averaging  about  23  deg.  The  rudder  trace  indi¬ 
cates  that  the  yaw  damper  is  on.  The  wing  rock  starts  as  a  divergent 
dutch  roll  but  then  appears  to  limit  at  approximately  40-50  deg  roll 
amplitude  with  a  $/8  ratio  between  2  and  3.  Also  the  divergence  grows 
linearly,  instead  of  exponentially  as  it  would  at  constant  c  •  During 
the  large,  constant-amplitude  4>  and  3  oscillation,  the  AOA  trace  shows  a 
small  oscillation  at  twice  the  frequency  of  the  lateral  oscillation. 

Figure  4  is  the  F-4J  simulation  trimmed  to  aQ  ■  23  deg  and  disturbed 
by  a  lateral  stick  pulse.  Here  the  dutch  roll  also  diverges  almost 
linearly  for  about  three  cycles  and  then  becomes  limited  in  amplitude. 
The  4/8  ratio  is  about  3  during  the  limit  cycle.  An  oscillation  devel¬ 
ops  in  the  AOA  trace,  again  with  a  frequency  double  that  of  the  lateral 
oscillation.  The  characteristics  are  remarkably  similar  except  for  the 
period  of  the  oscillations.  The  higher  frequency  of  the  F-4E  oscilla¬ 
tions  is  most  likely  due  to  a  higher  dynamic  pressure.  In  Figure  3  the 
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Figure  5.  Departure  During  Left  Turn 


aircraft  was  pulled  up  to  a  *  23  deg  within  about  10  sec  and  then  held 
at  this  AOA  for  about  13  sec.  The  initial  Mach  number  was  not  Indi¬ 
cated.  In  Figure  4  the  model  was  trimmed  at  an  Initial  speed  corre¬ 
sponding  to  aQ  -  23  deg.  If  all  other  factors  were  exactly  the  same, 
the  difference  in  frequency  shown  here  would  be  accounted  for  with  only 
a  30  percent  difference  in  speed. 

The  computer-generated  images  of  Figure  5  show  a  comparison  of  the 

F-4E  aircraft  and  the  F-4J  model  departures  from  a  wind-up  turn  to  the 

left.  The  departure  onset  is  quite  similar.  Both  depart  with  a  nose 
slice  away  from  the  turn  accompanied  by  a  roll  oscillation.  The  F-4E 
ends  up  with  more  of  a  rolling  departure,  while  the  F-4J  ends  with  more 
of  a  yaw  departure  (nose  slice).  This  difference  in  final  motion  is 
probably  influenced  greatly  by  the  differences  in  inertias.  The  F-4E 
has  24  percent  higher  yaw  inertia  and  7.3  percent  higher  ratio  of  yaw  to 

roll  inertia.  Thus  it  could  be  expected  to  be  more  resistant  to  yaw 

departure  and  prone  to  roll  departure. 
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Since  our  main  interest  is  in  modeling  and  establishing  causal  rela¬ 
tionships  leading  to  high  AOA  maneuver-limiting  factors  rather  than 
modeling  post-stall  gyration  and  spin,  the  foregoing  match  of  the  non¬ 
linear  wing  rock  behavior  and  nose  slice  onset  is  considered  to  be  one 
validation  of  the  aerodynamic  model.  Appendix  I  of  Part  III  and  por¬ 
tions  of  the  analysis  and  simulation  to  be  presented  later  offer  addi¬ 
tional  validation. 

One  key  aspect  is  that  the  aerodynamic  model  demonstrates  both  limit 
cycle  wing  rock  and  nose  slice  departure  without  introducing  artificial 
hysteresis  in  aerodynamic  coefficients  (Reference  14),  drastic  changes 
in  roll  damping  (Reference  10),  or  contrived  elevator  control  asymmetry 
(Reference  15). 

3.  Dynamic  Analysis 

In  high  AOA  maneuvering  flight  it  is  not  uncommon  for  the  aircraft 
to  be  in  asymmetric  (6  j  0 )  flight  —  either  intentionally  or  uninten¬ 
tionally.  Previous  studies  (e.g.,  References  5,  16-18)  have  shown  that 
static  aerodynamic  cross-coupling  due  to  sideslip  can  have  significant 
influence  on  aircraft  stability  and  dynamic  characteristics  and  possibly 
the  departure  characteristics.  For  example,  Figures  6-8  reflect  the 
pitch,  roll,  and  yaw  aerodynamic  moments,  respectively,  for  our  F-4J 
model  as  a  function  of  a  and  8.  For  simplicity  the  aerodynamic  moments 
were  assumed  symmetric  with  sideslip.  Over  the  region  15  <  a  <  25  deg 
pitching  moment  is  significantly  Influenced  by  sideslip,  and  rolling  and 
yawing  moments  due  to  sideslip  vary  greatly  with  AOA.  This  is  the  same 
AOA  region  in  which  the  various  maneuver-limiting  characteristics  of  the 
F-4  are  exhibited.  Thus,  in  the  analysis  to  follow,  attention  will  be 
devoted  to  examining  the  possible  influence  of  static  aerodynamic  cross- 
coupling  on  high  AOA  maneuver-limiting  factors.  Aircraft  open-loop 
stability  characteristics  will  be  analyzed  first,  and  then  cloBed-loop 
piloted  control  aspects  will  be  determined. 
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Figure  6.  F-4J  Pitch  Coefficient  Variation 
with  a,  $ 
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Figure  8.  F-4J  Yew  Coefficient  Variation 
with  a,  8 


a.  Open-Loop  Dynamic  Characteristics 

Figure  9  presents  a  locus  of  longitudinal  short-period  an  1  lateral 
dutch  roll  root  aigrations  with  AOA  in  wings-level,  symmetric,  1-g 
flight  at  an  altitude  of  15,000  ft.  These  dynamic  parameters  represent 
symmetric,  linearised,  fixed-operating-point  conditions.  With 
Increasing  AOA  the  dutch  roll  mode  is  seen  to  progress  from  a  lightly 
damped  stable  mode  to  oscillatory  divergence  and,  finally,  to  a  pair  of 
aperiodic  divergences  while  the  longitudinal  short  period  remains 
relatively  unchanged. 

Also  shown  in  the  figure  is  the  AOA  at  which  the  dynamic  stability 
parameter  Cng^yn  becomes  zero.  At  lower  AOA  the  parameter  is  positive, 
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Figure  9.  F-4J  Lateral-Longitudinal  Short-Period  Root 
Migration  with  a;  0O  “  0 

Indicating  directional  dynamic  stability.  At  high  AOA  the  parameter  la 
negative,  indicating  directional  dynamic  instability.  The  AOA  regions 
below  and  above  Cn3jyn  ■  0  are  identified  as  wing-rock  and  nose-slice 
regions,  which  correlate  quite  well  with  the  AOA  rangeB  Identified  with 
these  characteristics  in  the  F-4E  flight  test,  Figure  10  (from  Refer¬ 
ence  19 ) . 

The  influence  of  non-zero  sideslip  at  trim  on  the  dutch  roll  and 
short-period  frequencies  and  damping  1b  reflected  in  Figure  11.  The 
6  DOF  aircraft  was  trimmed  to  l.S  and  5.5  deg  sideslip  and  linearized 
transfer  function  parameters  obtained.  These  indicate  that  the  dutch 
roll  is  destabilized  and  the  short-period  stabilized  by  steady  sideslip. 
Thus  there  is  an  apparent  interchange  of  damping  between  the  two  moden 
while  the  overall  system  damping  remains  essentially  constant. 

Unfortunately,  comparison  of  dynamic  characteristics  predicted  by 
linear  frozen-point  analysis  in  Figures  9  and  11  with  the  actual  flight 
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Symbol 

Gross  Weight 
(lb) 

Load  Factor 
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Figure  10.  F-4E  Stall  Approach  Character  1st  ice 
(From  Ref.  19) 
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Figure  11.  F-4J  6  DOF  Linearized  Equation;  Lateral- 
Longitudinal  Root  Migration  with  a  and  3 


traces  of  Figure  3  or  nonlinear  6  DOF  model  traces  of  Figure  4  shows 
several  significant  differences.  The  wing  rock  linear  analysis  predicts 
a  dutch  roll  divergence  of  increasing  amplitude  as  a,  3,  or  both  ate 
Increased,  and  a  longitudinal  short-period  with  a  frequency  quite  close 
to  the  dutch  roll  frequency,  but  heavily  damped.  The  traceB  of 
Figures  3  and  4  show  the  dutch  roll  divergence  bounded  in  amplitude  and 
a  low-damped  longitudinal  oscillation  at  about  twice  the  dutch  roll 
frequency. 

Figure  12  shows  a  similar  set  of  nonlinear  time  response  character¬ 
istics  for  the  open-loop  6  DOF  nonlinear  model  trimmed  to  aQ  -  21  deg 
and  3Q  ■  0  deg  and  excited  with  an  aileron  pulse.  From  Figure  12  it  is 
apparent  that  the  longitudinal  oscillation  results  from  rectifying  the 
lateral  oscillation  and  the  sharp  peaks  of  pitching  acceleration 
coincide  with  sideslip  passing  through  zero.  Thus,  the  nonlinear-type 
behavior  appears  to  be  caused  by  lateral-longitudinal  coupling  due  to 
sideslip. 

A  nine-by-nine  matrix  of  aerodynamic  and  kinematic  terms  for 
coupled,  non-symmetric  flight  obtained  from  the  partial  derivative 
expansion  of  the  nine  equations  of  motion  is  presented  in  Figure  13. 
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Figure  13.  Matrix  for  Coupled  Equations  with  0O  t  0 


The  major  derivatives  and  kinematic  terms  are  identified  in  literal  form 
and,  for  comparison  with  Figure  12,  have  been  evaluated  at  aQ  «  21  deg 
aad  30  ■  l.S  deg.  The  principal  cross-coupling  at  this  small  sideslip 
angle  Is  seen  to  be  due  to  Mg,  «C  £  and  N£.  The  Influences  of  these 

derivatives  on  actual  vehicle  dynamic  response  can  best  be  seen  through 

the  use  of  eigenvectors  and  vector  polygons. 

Figure  14  contains  the  eigenvectors  for  the  dutch  roll  and  short- 
period  modes.  These  show  both  modes  to  be  dominated  by  roll,  side¬ 
slip,  and  pitch,  with  the  difference  between  the  modes  being  primarily 

motion  phasing  and  damping.  For  example,  in  the  dutch  roll,  q  and  8  are 
nearly  180  deg  out  of  phase,  while  in  the  short -period  they  are  nearly 
in  phase.  However,  the  short-period  is  so  well  damped  (Cap  *  0.6)  that 
It  would  Influence  only  the  first  second  or  two  of  the  time  responses 
shown  in  Figure  12.  Therefore,  the  observed  time  responses  must  be  due 
entirely  to  the  dutch  roll. 

Vector  polygons  for  the  four  principal  motion  equations  are  pre¬ 
sented  in  Figure  15.  Each  polygon  represents  a  force  or  moment  equa¬ 
tion.  The  polygon's  sides  show  the  contributions  of  each  motion  varia¬ 

ble  to  the  total  acceleration.  Vector  lengths  are  a  product  of  the 
appropriate  stability  derivatives  and  eigenvectors  evaluated  at  the 

dutch  roll  frequency.  The  predicted  motion  is  mostly  rolling  (p),  with 
•  « 

sideslip  (8),  and  pitching  (q)  about  the  same  magnitude,  and  with  com- 

•  as 

paratively  little  yawing  (r).  The  polygons  for  a  and  u  are  so  small 

that  they  can  be  neglected.  The  vector  polygons  provide  the  first 

insight  to  the  causae  for  the  change  in  modal  response  characteris¬ 
tics.  The  p  vector  polygon  shows  phasing  of  the  ta a  vector  opposite 
to  that  for  the  <£pp  vector.  ThuB  the  aerodynamic  cross-coupling 
coefficient,  in  effect  opposes  the  normal  aerodynamic  damping 

term,  tp,  and  therefore  increases  the  rolling  tendency  of  the  aircraft. 
Since  the  aerodynamic  cross-coupling  coefficients  Increase  in  magnitude 
as  the  trim  sideslip  increases,  one  can  readily  visualize  that  at  higher 
sideslip  angles,  ta  can  actually  become  greater  than  £p>  Then  one 
could  expect  the  model  to  exhibit  a  roll  divergence  for  large  a  and  8 
perturbations . 
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Figurt'  1A.  F-AJ  Eigenvectors;  aD  -  21  deg,  30  1,5  d®8 
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Figure  15.  F-4J  Force  Vector  Folygons; 
aQ  ■  21  deg,  $0  -  1.5  deg 


The  3  vector  polygon  indicates  that  this  vehicle  motion  derives  pri¬ 
marily  from  the  kinematic  coupling  of  roll  rate  about  the  body  center- 
line  axle  at  the  elevated  AOA. 

In  the  r  vector  polygon  the  p  vector  la  mlsaing  because  the  aero¬ 
dynamic  coefficient  Np  la  negligible  at  this  flight  condition.  The 
vector  N£a  has  a  relatively  large  component  in  phase  with  N£r  so  that 
the  aerodynamic  croas-coupllng  from  longitudinal  into  lateral  augments 
the  natural  yaw  damping  and  thus  contributes  to  reducing  body-axis 
yawing  motion.  All  of  this  contributes  to  the  dominance  of  roll  and 
sideslip  in  the  Figure  12  wing  rock. 

e 

Finally  the  q  vector  polygon  shows  that  the  MgB  vector  is  of  oppo¬ 
site  phase  and  is  larger  in  magnitude  than  the  damping  sector  derived 
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from  the  aerodynamic  derivatives  Mq  +  M£«  Thus  the  sideslip-induced 
pitching  moment  overpowers  the  normal  aerodynamic  damping  and  forces  the 
double  frequency  oscillation  observed  in  the  time  traces.  Mg  changes 
sign  with  0  (Figure  6).  Since  pitch  inertia  is  considerably  larger  than 
roll  inertia  for  this  aircraft,  a  transfer  of  energy  from  the  unstable 
lateral  mode  into  the  stable  longitudinal  mode  may  be  a  contributing 
factor  in  bounding  the  lateral-directional  mode  divergence.  The  influ¬ 
ence  of  reduci.1  •  is  shown  in  Figure  16.  In  this  case,  is  small 
but  not  zero.  Comparing  these  traces  with  those  of  Figure  12  demon¬ 
strates  the  strong  influence  Mg  has  in  bounding  the  lateral  divergence. 

It  is  concluded  for  our  model  that  the  roll/sideslip  oscillation 
starts  as  a  divergent  dutch  roll  with  small  sideslip  excursions.  As 
sideslip  ‘increases,  coupling  causes  a  nonlinear  behavior  which  bounds 
the  divergence.  Thus  F-4J  wing  rock,  per  se,  1b  a  nonlinear  phenomenon 
dominated  by  static  aerodynamic  lateral-longitudinal  cross-coupling. 
Further,  the  vector  polygons  show  that  except  for  H  £  the  damping 
derivatives  do  not  have  a  significant  influence  on  high  AOA  dynamic 
characteristics. 

This  analysis  has  also  shown  that  one  must  exercise  caution  in  the 
use  or  Interpretation  of  open-loop  transfer  function  parameters  even 
when  obtained  from  6  DOF  fixed  operating  point  conditions.  Such  param¬ 
eters  are  limited  to  relatively  small  sideslip  conditions  Buch  as  onset 
to  wing  rock  or  departure  and  hence  better  reflect  susceptibility  to 
these  maneuver-limiting  factors  rather  than  describing  the  resulting 
vehicle  motion. 

b .  Closed-Loop  Dynamic  Characteristics 

The  closed-loop  dynamic  characteristics  of  any  system  are  strongly 
Influenced  by  transfer  function  numerator  roots  (zeros)  in  or  near  the 
range  of  the  desired  bandwidth  of  control.  When  controlling  any  motion 
variable  with  any  controller,  the  numerator  zeros  of  that  transfer  func¬ 
tion  attract  denominator  (open-loop)  poles.  Hence  these  zeros  alter 
closed-loop  stability  and  modal  response  characteristics  of  the  vehicle 
in  proportion  to  the  tightness  (gain)  of  the  loop  closure.  Thus  certain 
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Figure  16.  Open-Loop  Wing  Rock  with  Cng  Small; 
00-21  deg,  30  ■  0  deg 


zeros  can  become  high  AOA  maneuver-limiting  factors  in  that  they  cause 
or  influence  closed-loop  roll  reversal,  wing  rock  (lateral  PIO),  and 
departure  (nose  slice  or  roll).  This  observation  applies  to  loops 
closed  both  by  the  pilot  and  by  automatic  feedback  systems.  The  key 
numerators  for  piloted  control  are  N^,  and  The  F-4J 

augmentation  systems  are  of  little  concern  due  to  limited  SAS  authori¬ 
ties  and  very  low  loop  gains  at  high  AOA  and  low  dynamic  pressure. 

Lateral  stick  control  of  roll.  Under  normal  low  AOA  flight 
situations,  the  3  DOF  equations  produce  a  transfer  function 

N? 

$(s)  -  "stk 


stk 


A^[s2  +  2c(j>ii)^s  + 


A  -  (s  +  l/Tg)(s  +  1/Tr)(b2  +  2cd^d“  +  <4) 

with  values  of  and  5^  in  the  vicinity  of  those  for  the  dutch  roll 
mode,  u)c;  and  With  this  loop  closure  the  zeros  attract  (modify) 

the  dutch  roll,  causing  the  closed-loop  root  locus  (with  increasing 
feedback  gain  K)  to  take  on  a  form  typified  by  the  sketch  below. 
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As  AOA  Is  increased,  generally  decreases  and  can  become  negative 
in  value.  When  this  occurs,  1*  of  the  form 

■  A+(* +^r)(*  +  t^) 

where  one  root  ie  positive  and  one  le  negative.  With  the  high  AOA  shift 
in  open-loop  denominator  roots,  the  root  locus  often  takes  on  the  form 
shown  below. 


"fy>2  TR  Ts  T*| 


In  this  case  the  spiral  root,  -1/Tg,  is  driven  toward  the  right 
half-plane  (RHP)  aero  and  a  first-order  instability  (divergence)  re¬ 
sults.  The  rate  of  divergence  depends  upon  how  far  the  zero  lies  in  the 
RHP  and  how  tightly  the  loop  is  closed. 

Migration  of  the  F-4J  open-loop  zeros  with  AOA  at  zero  sideslip 

is  reflected  by  the  solid  lines  in  Figure  17,  which  Indicate  an  almost 

2 

linear  decrease  in  w,*,  (at  constant  damping)  with  increase  in  AOA  up  to 
•  2 

a  -  18  deg.  Above  this  angle,  is  negative  and  the  RHP  root  location 
is  quite  AOA-sensitive.  Also  noted  on  the  figure  are  the  regions  asso¬ 
ciated  with  wing  rock  and  nose  slice  from  Figure  10.  Wing  rock  occurs 
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in  the  region  where  is  small;  nose  slice  occurs  in  the  region  where 
1/T^  is  large  and  negative. 

2 

Actually,  negative  ^  reflects  steady-state  roll  opposite  to  that 
commanded  by  stick  deflection.  This  is  generally  due  to  large  adverse 
aileron  yaw,  since  (after  Reference  20) 


where  -  86fl  +  11.66ep.  At  high  AOA,  Ng  is  small,  is  large 

negative,  and  N{  ^  determines  the  sign  of  Adverse  (negative)  Ngstk 
then  tends  to  produce  negative  and  the  RHP  zero.  In  physical  terms, 
adverse  yaw  produces  sideslip  and,  in  turn,  roll  reversal  due  to 
It  thus  becomes  apparent  that  the  nose  slice  region  identified  in  Fig¬ 
ure  17  may  be  related  to  strong  adverse  aileron  yaw  and  the  magnitude  of 
the  RHP  zero. 

2 

The  dashed  lines  of  Figure  17  show  the  values  of  and  2c^tify  (or 

1/T^,  l/T^)  for  b0  ■  5.5  deg.  Whereas  increasing  AOA  at  8  ■  0  causes 

u>4  to  decrease  steadily  and  finally  become  negative,  when  8^0  the 
*  2 

aerodynamic  cross-coupling  causes  to  stay  positive  raid  close  to  the 
dutch  roll  mode  but  the  damping  ratio,  becomes  negative  (see  sketch 
below) . 


Lateral  -  Directional 
Root  Locations  for : 

a0  =  2i° 

0o  =  55° 


WSH 
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The  promixity  of  and  means  low  residue  for  the  dutch  roll  mode 
In  roll  control  with  stick  when  {5  f  0.  Thus  sideslip  should  reduce  the 
roll  content  of  the  unstable  dutch  roll  mode;  this  is  consistent  with 
the  observed  bounding  of  the  roll  oscillation  divergence  in  Figures  4 
and  16. 

Rudder  control  of  roll.  When  the  region  of  adverse  aileron  yaw  Is 
reached,  roll  control  is  maintained  via  rudder.  This  produces  yaw  and 
sideslip-induced  roll  in  the  same  direction  and  thus  much  more  rapid 
turning  than  just  yawing  moment  and  side  force  would  produce.  The 
rudder-to-roll  numerator,  ,  is  of  the  form 


One  root  Is  always  positive  (RHP)  as  shown  in  Figure  18.  The  RHP  zero 
results  from  adverse  roll  due  to  It  is  apparent  from  Figure  18 

that  the  zeros  of  n|  are  not  greatly  influenced  by  either  AOA  or  side¬ 
slip. 

Again  this  control  structure  portends  closed-loop  instability  since 
some  pole  must  be  driven  toward  the  RHP  zero.  However,  these  zeros  lie 
at  high  frequencies  relative  to  the  low-frequency  bandwidth  normally 
achieved  with  rudder  control  via  the  pilot's  legs.  Rudder-to-roll  con¬ 
trol  handling  qualities  are  normally  dominated  by  the  low-frequency 
spiral  and  roll  subsidence  modes,  and  at  high  AOA  these  remain  stable 
and  well  behaved.  While  this  type  of  control  is  "unnatural,"  it  does 
not  result  in  roll  stability  problems  at  low  gain  closures. 


On  the  F-4J,  I/TQ2  decreases  with  increasing  AOA  and  finally  couples 
with  1/Tq^  to  form  a  complex  pair  which  almost  cancels  the  longitudinal 
phugoid  mode.  Thus  these  zeros  remain  well  behaved.  However,  when 
there  is  sideslip  so  that  static  aerodynamic  cross-coupling  becomes 
significant,  the  numerator  contains  additional  zeros  related  to  lateral- 
directional  modes. 

It  was  discovered  (Reference  5)  that  such  coupling  in  the  A-7  air¬ 
craft  produced  an  RHP  zero,  1/Tq^,  and  subsequent  pilot  control  of  pitch 
attitude  via  the  stabilator  caused  a  closed-loop  directional  divergence 
which  has  the  appearance  of  a  nose  slice.  A  similar  situation  exists 
in  the  F-4J  aircraft.  Figure  19  presents  pole  and  zero  locations  for 
aQ  -  23  deg  and  »  5.S  deg.  The  coupled  6  DOF  equations  (Figure  19a) 
produce  two  real  zeros  (one  RHP,  one  LHP),  one  complex  pair  which  nearly 
cancels  the  phugoid,  and  one  higher-frequency  complex  pair. 

The  source  of  the  low-frequency  zeros  can  be  determined  by  elimin¬ 
ating  the  velocity  equation,  which  should  reduce  the  phugoid  mode  to 
a  first-order  pole  at  the  origin  and  eliminate  the  speed  numerator  zero, 
l/Te^.  It  may  be  observed  In  Figure  19b  that  the  only  change  is  in  the 
phugoid  dipole  pair  which  became  first-order  as  expected. 

Next,  as  in  Reference  5,  the  cross-coupling  terms  and  N£  are  set 
to  zero  in  the  6  DOF  equations  to  identify  the  lateral-directional 
modes.  This  is  shown  in  Figure  19c  where  the  pole-zero  cancellations 
show  the  two  real  zeros  to  be  related  to  the  roll  subsidence,  1/Tr,  and 
spiral,  1/T#,  modes  and  the  higher-f requeucy  complex  zeros  to  be  related 
to  the  dutch  roll. 

Figure  20  is  a  system  survey  for  closure  of  the  pitch  attitude  loop 
A 

for  the  augmented  (SAS  on)  airframe  at  this  same  trim  condition.  The 
lower  figure  is  a  root  Iocub  and  the  upper  figure  is  a  Bode-siggle 
(Reference  20).  The  root  locus  shows  the  migration  of  roots  from  the 


The  slight  difference  In  airframe  dynamic  parametere  between 
Figures  20  and  11  indicates  the  very  minor  Influence  of  the  SASs  at  high 
AOA  and  low  dynamic  pressure  conditions. 
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open-loop  poles  toward  the  open-loop  zeros  as  gain  Is  Increased .  It 
provides  Information  concerning  closed-loop  frequency  and  damping  as  a 
function  of  gain  but  does  not  provide  Information  concerning  sensitivity 
to  gain  change.  The  Bode-siggle  presents  the  open-loop  transfer  func¬ 
tion  amplitude  and  phase  as  functions  of  frequency  where  the  amplitude 
is  evaluated  for  s  *  jw  (Bode)  and  s  ■  o  (Biggie).  The  -f»  symbols  iden¬ 
tify  migration  of  complex  LHP  closed-loop  pole  frequency  with  increasing 
gain  (Bode-root  locus,  Reference  20)  and  the  0  symbols  similar  migra¬ 
tion  of  closed-loop  piles  in  the  RHP.  Closed-loop  first-order  roots 
lie  at  the  intersection  of  the  gain  closure  lines  and  the  o  amplitude 
curves.  Thus  these  survey  plots  allow  identification  of  closed-loop 
roots  at  a  specific  gain  •-  ilue,  and  sensitivity  of  root  value  to  changes 
in  gain. 

In  order  to  provide  any  pitch  attitude  control  whatsoever,  the  loop 
gain  must  he  sufficiently  high  that  the  closure  line  crosses  the  Bode 
plot  below  the  -20  dB  low-frequency  asymptote,  i.e.,  with  a  gain  20  dB 
or  greater.  Generally  the  crossover  will  be  at  or  above  the  short- 

period  frequency  -  again  in  the  vicinity  of  -30  to  -40  dB.  The  root 

locus  plot  has  tick  marks  showing  closed-loop  root  locstions  for  gains 
of  20  and  32  dB.  This  shove  that  the  root  moving  toward  l/Tg^  has 
already  reached  the  RHP,  and  divergence,  with  the  gain  at  20  dB.  By 
32  dB  the  first-order  divergence  has  a  time  constant  of  4  sec.  Thua  any 
attempted  control  of  pitch  attitude  results  in  a  divergence.  This 
divergence  is  due  to  lateral-longitudinal  coupling  and  is  evidenced  in  a 
lateral  motion. 

Figure  21  shows  values  for  the  real  zeros  1/Tg.j  and  1/TR  at  @0  - 

5.5  deg  and  aQ  from  21  to  25  deg.  These  indicate  that  at  non-zero 

sideslip  and  A0A  increasing  above  20  deg  I/Tq^  rapidly  moves  to  negative 

values  which  portend  closed-loop  first-order  divergence  time  constants 
as  low  as  1  sec.  Again  this  coincides  with  the  nose  slice  region 
observed  in  flight  test. 
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Figure  21.  N$  ^  First-Order  Root  Migration* 
with  AOA;  F-4J 

Coupling  numerator* .  when  the  eircreft  ia  in  asymmetric, 
sideslipping  flight,  any  control  action  in  the  longitudinal  axis 
influences  the  lateral-directional  numerator  and  vice  versa.  This  can 
be  illustrated  as  follows: 


The  effective  transfer  function  for  each  control  loop  becomes 


»8 

N«atab 

__ 


+ 


stab°stk 


stab°stk 


where  N«atab*.tk  is  the  numerator  term  due  to  cross-coupling  effects.* 
Both  the  0  and  $  numerators  are  modified  by  the  same  coupling  numerator; 
the  extent  of  the  influence  depends  upon  how  tightly  the  other  loop  is 
being  closed.  This  Influence  can  be  determined  by  treating  the 
effective  numerator  as  a  closed-loop  system  and  identifying  the 
migration  of  numerator  roots  as  the  loop  gain  (Y^  or  Ye)  is  varied , 
e.g.,  the  roots  of  the  6  numerator  ars  obtained  from 
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The  uncoupled  numerator  roots  become  the  poles  of  this  system  and 
the  coupling  numerator  roots  are  the  zeros .  The  root  migration  is 
controlled  by  the  pilot  gain,  Y^,  For  the  Figure  20  example  case  of 
aQ  “  23  deg,  0O  *  5.5  deg,  the  transfer  function  is 
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The  reader  is  referred  to  page  163  et  seq.  of  Reference  20  for  a 
thorough  development  and  treatiae  of  multiloop  control  systems  analysis. 
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where  first-order  factors  (s  +  a)  are  represented  by  (a)  and  second- 
order  factors  [s  +  2cmns  +  u£]  by  [5;  unJ .  Assuming  to  be  a  simple 

positive  gain,  this  produces  a  root  locus  of  the  form  shown  in  the 
sketch  below.  Thus  pilot  control  of  the  roll  attitude  loop  via  aileron 
and  spoiler  causes  the  RHP  zero  of  the  pitch  numerator  to  move  further 
into  the  RHP  because  the  root  of  the  coupling  numerator  lies  beyond 

i/t03. 


Similarly,  pilot  control  of  pitch  attitude  via  Btabilator  alters  the 
roll  numerator  as 
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which,  for  Yg  assumed  a  simple  nsgative  gain,  produces  a  root  locus  of 
the  form  sketched  below.  This  results  in  the  RHP  complex  pair  being 
driven  to  the  real  axis  with  one  root  eventually  moving  further  into  the 
RHP  while  the  other  moves  into  the  LHP. 


<  I  5  v  •,  1  ti  >  I  JW  1 . 
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Thus  in  each  instance  the  influence  of  the  coupling  numerator  is 
definitely  not  beneficial;  the  coupling  can  further  degrade  closed-loop 
control  in  both  axes. 

Since  F-4  pilots  are  normally  instructed  to  use  rudder  to  roll  at 
high  AO  A,  the  Influence  of  this  control  technique  may  be  observed  from 
the  transfer  functions  (again  at  aQ  ■  23  deg,  (5  ■  5.5  deg): 


Use  of  rudder  to  roll  modifies  the  pitch  numerator  in  much  the  same 
manner  as  does  use  of  aileron  and  spoiler.  This  can  be  aeen  by 
comparing  values  of  the  coupling  numerator  terms.  On  the  other  hand, 
pitch  control  has  little  Influence  on  the  rudder  roll  numerator  because 
nearly  all  numerator  and  denominator  roots  cancel.  Thus  it  appears  that 
this  control  strategy  produces  the  most  consistent  closed-loop  lateral- 
directional  dynamic  properties  at  high  AOA  and  non-zero  sideslip,  which 
also  might  be  a  reason  pilots  do  not  mind  using  this  technique. 

Summary.  This  enalysls  has  shown  several  areas  in  which  control- 
numerator  dynamic  parameters  can  become  maneuver-limiting  factors  for 
the  F-4J.  These  limits  are  generally  associated  with  zeros  migrating  to 
the  RHP  at  high  AOA  in  asymmetric  (3  4  0)  flight. 

The  RHP  zero  of  is  due  to  adverse  "aileron"  yaw  (N^^  nega¬ 
tive)  and  results  in  roll  reversal.  The  attendant  sideslip  causes  the 
4  zeros  to  shift  toward  the  dutch  roll  poles.  This  proximity  reduces 
the  roll  contribution  in  the  dutch  roll  mode  and  produces  a  nonlinear 
bounding  of  dutch  roll  (or  wing  rock)  oscillation  amplitude. 

One  RHP  zero  exists  in  N$  whether  the  aircraft  is  in  asymmetric 

°r 

flight  or  not.  However,  that  root  is  at  a  frequency  sufficiently  above 
the  region  in  which  rudder  control  is  generally  exerted  that  it  does  not 
pose  a  closed-loop  stability  problem. 
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An  RHP  zero  can  occur  in  N68ta^  due  to  aerodynamic  cross-coupling 
with  non-zero  sideslip*  This  zero  results  in  a  first-order  divergence 
mode  if  stabilator  control  of  pitch  attitude  le  attempted. 

Asymmetric  flight  also  produces  lateral-longitudinal  control  cou¬ 
pling  that  further  aggravates  these  undesirable  RHP  zero  locations  of 
the  key  €  and  $  numerators. 

Thus  the  aerodynamic  cross-coupling  associated  with  sideslip  tends 
to  dominate  the  F-4J  high  AOA  flying  qualities.  Values  of  the  cross- 
coupling  derivatives  evaluated  at  B0  ■  5.5  deg  and  a  range  of  AOA  are 
shown  in  Figure  22.  Comparing  these  with  the  root  migrations  in 
Figures  9  and  11  demonstrates  their  significance.  The  regions  where 
large  shifts  in  dynamic  parametars  occur  coincide  with  large  magnitude 
of  cross-coupling  derivatives.  Also,  it  should  be  noted  that  the 
magnitude  of  the  cross-coupling  derivatives  is  directly  proportional  to 
the  magnitude  of  sideslip. 

The  foregoing  analytic  results  must  be  viewed  with  some  caution 
since  it  was  shown  that  linearized,  fixed-operating-point  dynamic 
parameters  do  not  necessarily  reflect  the  actual  open-loop  dynamics  of 
the  nonlinear  airframe.  On  the  other  hand,  the  apparent  correlation 
between  bounaing  of  dutch  roll  divergence  and  movement  of  zeros  with 


Figure  22.  F-4J  Static  Aerodynamic  Cross-Coupling 
Derivatives;  0n  *  5.5  deg 
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sideslip,  and  Che  AOA  region  for  nose  slice  observed  In  actual  flight 

with  the  region  where  the  RHP  zeros  of  8 tab  an<*  N^stk  are  *-ar8®  18 
quite  encouraging . 


B.  F-14A  ANALYSIS 

The  F-14A  configuration  does  not  have  wing  leading-edge  elats.  This 
aircraft  has  variable  wing  sweep;  at  M  <  0.5  the  wings  are  in  the  full 
forward  (A  ■  22  deg)  position.  The  unaugmented  airframe  is  classified 
as  departure  resistant,  but  it  can  be  departed  (References  21  and  22). 
It  remains  laterally  stable  and  longitudinally  controllable  at  full  aft 
stick  (AOA  approximately  50  deg  with  wings  forward).  It  has  a  mild 
dutch  roll  instability  starting  at  about  15  deg  AOA.  The  oscillations 
are  aggravated  if  the  pilot  attempts  to  oppose  the  motion  with  either 
lateral  stick  or  rudder.  The  dutch  roll  mode  becomes  stable  again  at 
about  23-25  deg  AOA.  Roll  reversal  due  to  adverse  yaw  from  the  aileron 
and  differential  horizontal  stabilizer  also  starts  at  about  18  deg 
AOA.  Prolonged  deflection  of  the  lateral  stick  at  a  >  deg  will 
result  in  a  departure,  characterized  by  a  rapid  Increase  in  “adverse'’ 
yew  rate,  which  may  develop  into  a  series  of  uncommanded  rolls  with  the 
appearance  of  snap  rolls.  At  low  speed,  departure  can  occur  at  AOA  as 
low  as  21  deg.  The  departure  characteristics  are  a  function  of  wing 
sweep,  c.g.  location,  and  Mach  number.  The  airplane  will  spin  and  has  a 
non- recoverable  flat  spli  mode. 

1.  Analytic  Modal 

The  F-14A  aerodynamic  model  was  derived  from  several  sources  (Refer¬ 
ences  23-27).  The  data  base  was  taken  from  the  NASA/Langley  Differen¬ 
tial  Maneuvering  Simulator  (DMS)  report  of  Reference  23.  This  fixed- 
base  piloted  simulation  was  considered  by  Grumman  test  pilots  to  provide 
a  valid  representation  of  high  AOA  dynamic  characteristics  observed  in 
early  flight  tests.  However,  time  traces  obtained  from  thor-o  same  aero¬ 
dynamic  data  did  not  match  available  flight  test  records.  Therefore 
the  data  were  modified  on  the  basis  of  information  obt fined  from  other 
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references!  discussions  with  F-14  aerodynamlciats  at  Grumman,  and 
finally,  as  necessary  to  match  available  flight  test  data. 

The  resulting  aerodynamic  model,  together  with  supportive  documents'- 
tlon  for  each  aerodynamic  coefficient,  is  detailed  in  Part  III,  Appendix 
II.  The  data  caver  the  ranges 

0  <  a  <  55  deg;  -20  <  6  <  +20  deg 


at  low  Mach  (0.2).  At  this  low  speed  the  variable-sweep  wing  is 
nominally  in  its  full  forward  (A  ■  22  deg)  position,  so  wing  sweep 
influence  is  not  Included. 

Based  upon  the  F-4J  analysis,  the  key  high  AOA  aerodynamic  moments 
are  Cy,(a,g),  Cn(u»0)  end  Cm(a,0).  The  roll  moment  (Figure  2.3)  Is  seen 
to  remain  essentially  constant  for  a  >  15  deg  and  0  <  10  deg,  as  opposed 
to  the  F-4J's  significant  decrease  in  Gfc(a,0)  fora  >  15  deg  (Figure  7). 
The  yaw  moment  (Figure  24)  ia  similar  to  the  F-4J's  (Figure  8)  for 
a  <  17  deg,  but  for  a  >  17  deg  both  the  slopes  with  respect  to  a  and  the 
maximum  values  with  respect  to  0  era  approximately  half  those  of  the 
F-4J.  The  combination  of  less  negative  Cn^  and  larger  negative  Cig  pro¬ 
duces  positive  values  of  the  dynamic  stability  parameter  0ngdyn  through¬ 
out  the  55  deg  AOA  range. 

The  pitch  moment  with  sideslip  (Figure  25)  is  also  similar  to  that 
of  the  I-4J  ovar  the  range  10  <  ci  <  25  deg,  where  it  Is  negative.  How¬ 
ever,  above  a  ■  35  deg,  becomes  positive. 

Weight,  c.g.,  and  moments  and  product  of  inertia  for  the  clean  (no 
external  stores)  aircraft  were  provided  by  Grumman.  These  are  shown  in 
Table  2  with  F-4J  values  for  comparison.  All  inertias  are  In  the 
fuselage  centerline  reference  axis  system.  Aircraft’  dimensions  are  also 
presented  for  comparison. 
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Figure  25.  F-14A  Pitch  Coefficient  Variation 

with  a,  fl 


TABLE  2 


COMPARISON  OF  F-I4A  AND  F-4J  WEIGHTS, 
INERTIAS,  AND  DIMENSIONS 


Parameter 

F-14A 

F-4J 

W  (lb) 

46,950 

37,000 

c.g.  (X  c) 

16.0 

29.3 

Ix  (slug-ft2) 

58,950 

23,850 

Iy  (slug-ft2) 

225,600 

127,400 

Iz  (slug-ft2) 

285,000 

146,000 

Ixz  (slug-ft2) 

-3030 

2210 

A  (deg) 

22 

45 

CM 

4-» 

■N— ' 

to 

565 

530 

b  (ft) 

64.1 

38.67 

C  (ft) 

9.8 

16.04 

The  Inertia  differences  between  the  vehicles  are  large.  Although 
the  F-14A  Is  some  27  percent  heavier  than  the  F-4J  with  comparable 
fuel  loading,  the  inertias  range  from  77  to  147  percent  higher  and  the 
product  of  inertia  is  of  opposite  sign.  This  has  considerable  Influence 
on  the  ratios  of  aerodynamic  to  kinematic  terms  in  the  model  and  greatly 
influences  coupling  effects,  as  will  be  shown  later. 

Aircraft  control  is  exerted  through  a  horizontal  stabilizer  (6Btab), 
twin  rudders  (6r),  spoilers  and  differential  horizontal  (6D)  sur¬ 
face  deflections.  Simplified  block  diagrams  of  the  flight  control  and 
augmentation  systems  are  shown  in  Figures  26-28.  The  longitudinal  sys¬ 
tem  contains  a  feel  system  (spring  cartridge,  variable  gearing  ratio, 
bobweighta),  a  limited-authority  pitch  rate  SAS,  and  a  surface  actuator. 
However,  the  analysis  and  simulation  are  based  on  a  simple  fixed  stick- 
to-surfaoe  gain  and  pitch  rate  SAS. 
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Figure  26.  F-14A  Longitudinal  Control  System 


Figure  27.  F-14A  Lateral  Control  System 
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Figure  28.  F-14A  Directional  Control  System 


The  lateral  system  (Figure  27)  has  a  roll  rate  command  augmentation 
system  (CAS)  in  parallel  with  a  mechanical  stick-surface  link.  At  high 
AOA  the  spoilers  are  ineffective  and  all  roll  control  is  obtained 
through  differential  horizontal,  <$D.  The  roll  rate  CAS  operates  only 
through  the  differential  horizontal,  and  the  electrical  path  gains  are 
varied  with  AOA  such  that  the  CAS  command  and  feedback  fades  to  zero 
when  20  <  a  <  31  deg  and,  simultaneously,  the  electric  stick-surface 
path  fades  in  to  cancel  all  but  ±2  deg  of  mechanical  path  surface 
command.  Thua,  for  a  >  31  deg  there  is  no  roll  rate  damping  augmenta¬ 
tion  and  the  pilot  can  only  command  ±2  deg  of  differential  horizontal 
surface. 

The  rudder  system  (Figure  28)  contains  a  fixed  gain  mechanical  link 
from  pedals  to  surface  actuators,  a  conventional  yaw  damper  with  yaw 
rate  and  lateral  acceleration  feedback,  and  a  lateral  stlck-to-rudder 
crossfeed  (SRI).  At  a  >  10  deg  the  yaw  SAS  is  cut  out  and  the  SRI  faded 
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in.  The  SRI  signal,  passed  through  a  first-order  lag  filter  (break 
frequency  8  rad/sec),  reaches  full  gain  at  a  >  20  deg.  Above  a  ■  20  deg 
full  lateral  stick  deflection  produces  19  deg  of  rudder  deflection. 

Thus  for  the  AOA  region  of  Interest,  a  >  20  deg,  the  airframe  is 
essentially  unaugmented  in  the  lateral-directional  axes  and  the  control 
stick  effectively  moves  the  stabllator  for  pitching  and  twin  rudders  for 
rolling. 

2.  Model  Validation 

F-1AA  flight  test  data  covering  the  appropriate  AOA  region  were 
provided  by  the  aircraft  manufacturer  (Reference  28).  Unfortunately, 
the  data  were  extracted  from  tests  of  three  aircraft  having  slightly 
different  weights,  inertias,  flight  control  systems,  and  flight  test 
instrumentation.  Since  the  flight  test  goals  involved  flight  control 
and  leading  edge  flap  system  development  rather  than  identification  of 
basic  airframe  dynamics,  most  runs  were  contaminated  by  SAS  effects, 
pilot  control  inputs,  or  configuration  influences.  The  model  validation 
therefore  was  based  on  general  response  m«ffthirvg  of  a  few  short  time 
"windows"  of  bare-airframa  response  or  rel«M  'uly  uncontaminated  traces, 
plus  observation  that  the  sequence  of  dynamic  characteristic  changes 
with  increasing  AOA  described  previously  was  indeed  achieved. 

Figure  29  presents  results  of  a  static  longitudinal  validation  check 
based  upon  a  1-g  stall  in  which  AOA  is  increased  steadily  from  10  deg  to 
almost  40  deg  with  a  ramp  horizontal  stabilizer  input.  The  circled 
points  reflect  the  excellent  agreement  for  horizontal  stabilizer  inputs 
required  to  achieve  the  noted  AOA  trim  condition  at  35,000  ft  altitude 
with  the  aerodynamic  model. 

Figure  30  is  a  short  time-slice  of  the  divergent  dutch  rolx  at 
approximately  15  dsg  AOA.  Motion  is  excited  by  a  lateral  stick  doublet 
which,  through  the  SRI,  produces  both  differential  horizontal  and  rudder 
doublets.  Both  the  yaw  and  roll  rate  SASs  are  then  turned  off,  and  the 
subsequent  oscillation  represents  the  "bare"  airframe  in  the  lateral- 
directional  sense.  As  with  the  F-4J,  the  envelope  of  the  oscillation 
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Figure  29.  F-14A  Longitudinal  Static  Trim  Comparison 
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seems  to  grow  linearly  rather  than  exponentially,  which  Indicates 
nonlinear  damping.  AIbo,  the  AOA  trace  shows  an  oscillation  at  twice 
the  dutch  roll  frequency.  The  motion  on  the  trace  indicates  the 

pitch  SAS  was  on  and  responding  to  the  longitudinal  oscillation. 

Figure  31  presents  time  traces  for  the  F-14  aerodynamic  model 
trimmed  at  15  deg  AOA  and  excited  by  aileron  and  rudder  doublets  similar 
to  Figure  30.  This  also  shows  an  almost  linear,  rather  than 
exponential,  growth  in  dutch  roll  amplitude.  The  frequency  and  cycles 
to  double  amplitude  are  the  same  as  in  Figure  30.  The  small-amplitude 
AOA  oscillation  also  is  precisely  twice  the  frequency  of  the  dutch  roll. 

Figure  32  presents  another  short  time-slice  in  which  AOA  is 

maintained  at  18  ±  2  deg.  In  this  angle  of  attack  region  the  SAS  gain 

schedule  shows  that  the  yaw  SAS  should  be  off  and  the  roll  SAS  on. 

However,  comparison  of  the  yaw  rate  and  rudder  position  traces  indicates 
that  the  rudder  is  correlated  with  yaw  rate  at  approximately  the 
magnitude  and  phase  angle  that  would  result  from  the  yaw  damper.  Since 
the  rudder  pedal  and  lateral  stick  are  both  constant  end  very  small,  the 
aircraft  response  is  essentially  that  of  the  controls-fixed  augmented 
airframe.  This  dutch  roll  oscillation  is  essentially  the  same, 
indicating  that  the  yaw  SAS  is  quite  ineffective  at  this  AOA  and 

airspeed. 

Other  comparisons  which  can  be  made  among  Figures  30,  31,  and  32  are 


summarised  below: 

Parameter 

F-14A 

Model 

Figure  32 

Figure  30 

Figure  31 

Period  of  oscillation  (Bee) 

4.5 

4.5 

4.5 

Cycles  to  double 

2 

2 

2 

$/8  (deg/deg) 

2 

* 

2 

p/r 

25 

** 

22.5 

a  (deg) 

18  ±  2 

15 

15 

A  A  A 

3  trace  not  available.  r  trace  contaminated. 
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Taken  together,  the  foregoing  is  considered  an  acceptable  match  between 
actual  and  simulated  airframe  dynamics  in  the  region  of  15  deg  AOA  where 
the  dutch  roll  divergence  and  wing  rock  are  most  pronounced. 

Figure  33  shows  a  pullup  initiated  from  wings-level  flight  in  which 
AOA  is  steadily  Increased  from  13  to  38  deg  at  about  0.75  deg/sec.  A 
lateral  stick  oscillation  is  introduced  at  about  15  deg  AOA.  The  stick 
is  then  centered  and  the  SASs  turned  off  at  about  18  deg  AOA.  A  small 
rudder  ramp  is  then  introduced,  possibly  to  counter  the  Blow  roll-off  in 
bank  angle.  Simulation  traces  for  a  matching  pullup  with  similar  stick 
and  rudder  inputs  are  shown  in  Figure  34.  Both  Figures  33  and  34  show 
a  divergent  dutch  roll.  However,  the  amplitude  appears  to  increase 
linearly  rather  than  exponentially,  again  reflecting  nonlinear  damping. 
Both  sets  of  traces  show  the  period  of  oscillation  to  increase  from 
about  4.5  sec  to  over  5.5  sec  as  AOA  approaches  30  deg.  Again  there  is 
good  agreement  between  the  actual  flight  and  simulation  responses. 

3.  Dynamic  Analysis 

The  principal  shortcomings  of  the  F-14A  bare-airframe  flight  char¬ 
acteristics  at  high  AOA  are  wing  rock  (unstable  dutch  roll)  and  roll 
reversal.  Otherwise  the  airplane  is  quite  docile.  It  is  departure- 
resistant,  but  can  be  departed  if  lateral  control  deflection  1b  sus¬ 
tained  without  large  accompanying  rudder.  Thus,  the  following  discus¬ 
sion  will  be  devoted  to  identifying  potential  factors  underlying  these 
key  differences  between  the  F-14A  and  F-4J. 

a.  Open-Loop  Dynamic  Characteristics 

A  Burvey  plot  of  root  migration  with  angle  of  attack  for  the  range 
15  <  a  <  30  deg  and  8  ■  0  is  shown  in  Figure  35.  This  indicates  excel¬ 
lent  lateral-directional  agreement  with  the  previously  reported  high  AOA 
flight  characteristics  and  with  the  flight  traces  of  Figures  30,  31,  and 
33.  The  dutch  roll  mode  becomes  divergent  at  just  under  15  deg  AOA  and 
stays  unstable  until  approximately  22-23  deg  AOA.  This  is  the  region  of 
reported  wing  rock.  Above  23  deg  AOA  the  dutch  roll  is  stable  and  well 


Figure  35.  F-14A  Root  Migrations  with  AOA 


damped  but  the  coupled  roll-subsidence  and  spiral  becomes  oscillatory. 
The  AOA  at  which  Ng^^  ■  0  is  not  shown  on  this  plot  because  it  remains 
positive  (as  shown  in  Figure  36). 

The  longitudinal  short-period  mode  exhibits  relatively  low  frequency 
and  high  damping  at  15  deg  AOA.  Frequency  then  increases  and  damping 
decreases  with  increasing  AOA  until  the  dutch  roll  and  short-period  are 
about  equal  at  30  deg  AOA.  Note  that  this  is  not  in  agreement  with  the 
flight  traces  for  15  deg  AOA  (Figure  30),  which  showed  a  short-period 
frequency  twice  that  of  the  dutch  roll. 

Referring  back  to  Figure  25,  15  deg  AOA  is  where  Cmg  is  a  maximum. 
Thus,  cross-coupling  through  Mg  may  again  create  pitch  frequency  cou¬ 
pling  and  contribute  to  the  nonlinear  dutch  roll  damping,  as  was  the 
case  with  the  F-4J.  However  for  3^0  the  shift  in  the  Figure  35  root 
locations  is  so  small  as  to  be  negligible,  and  this  is  not  the  same  as 
with  the  F-4J. 
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Figure  36.  Variation  with  AOA;  F-14A 

Eigenvectors  for  the  coupled  dutch  roll  and  $hort-*period  modes 
of  the  6  DOF  airframe  at  aQ  *  20  deg  and  60  ■  1.5  deg  are  shown  in 
Figure  37.  There  is  essentially  no  longitudinal  coupling  into  the  dutch 
roll  mode  and  only  alight  lateral  coupling  (mainly  rolling  motion)  into 
the  short-period  mode.  These  eigenvectors  are  noticeably  different  from 
those  for  the  F-4J  (see  Figure  14)  at  a  similar  flight  condition. 

Vector  polygons  for  the  dutch  roll  mode  at  this  flight  condition 
are  shown  in  Figure  38.  The  p  equation  indicates  that  £.  g  and  3f '  com¬ 
pletely  dominate  rolling  motion.  All  damping  derives  from  £  there 
is  no  cross-coupling  influence  from  All  terms  in  the  i  equation 

have  been  doubled  to  make  this  polygon  legible.  The  relative  size 

o  * 

of  the  p  and  r  polygons  Indicates  that  in  this  AOA  region  static  and 
dynamic  yawing  force  and  moment  terms  are  negligible  in  the  dutch  roll 
mode.  However,  the  yawing  moment  that  does  exist  derives  primarily  from 
Np.  The  4  equation  has  significant  magnitude  compared  to  the  £  equa¬ 
tion.  The  key  is  the  dominance  of  the  M^lBl  term  which,  as  with  the 


59 


F-4J,  indicates  that  coupling  from  the  lateral  into  the  longitudinal 
produces  the  frequency  doubling  effect. 

Based  upon  the  above  argument,  it  ia  concluded  that  the  F-14A  wing 
rock  1b  due  to  a  divergent  dutch  roll  which  derives  from  very  low 
damping  in  roll,  £ p,  rather  than  lateral-longitudinal  cross-coupling 
(since  the  open-loop  lateral-directional  dynamic  characteristics  are 
relatively  insensitive  to  sideslip).  The  principal  static  aerodynamic 
cross-coupling  is  a  transfer  of  energy  from  rolling-sideslippiug  motion 
into  pitch  motion  and  this,  as  with  the  F-4J,  results  in  frequency 
doubling  into  pitch  and  may  contribute  to  amplitude  bounding  of  the  wing 
rock* 

b.  Closed-Loop  Dynamic  Characteristics 

Lateral  stick  control  of  roll.  Figure  39  presents  the  key 
parameters  for  the  numerator  NJ^.  The  lower  plot  indicates  that  the 
movement  of  one  zero  into  the  RHP  with  increasing  AOA  is  not  as  rapid 
for  this  airframe  as  for  the  F-4J  (see  Figure  17).  However,  departure 
can  occur  where  one  of  the  real  roots  is  large  negative.  This  plot  also 
demonstrates  that  the  roll  numerator  zeros  are  quite  insensitive  to 
sideslip  when  compared  to  the  F-4J. 

Rudder  control  of  roll.  Same  characteristics  as  for  the  F-4J. 

Stabilator  control  of  pitch.  The  potential  for  divergence  lies 

6 

with  the  RHP  zero,  I/Tq^,  in  Figure  40  is  a  plot  of  l/Tjj^  at 

80  «  4  deg  over  a  range  of  AOA.  Note  that  the  parameter  scale  factor  is 
an  order  of  magnitude  less  than  the  one  used  for  the  F-4J  in  Figure  21, 
and  the  largest  value  (approximately  -0.01)  is  two  orders  of  magnitude 
less  than  for  the  F-4J  (equal  to  -1.0).  Thus,  pilot  control  of  pitch 
attitude  should  have  no  noticeable  coupling  into  the  lateral-directional 
modes  for  the  F-14A. 
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Figure  40.  ^5ata|j  First-Order  Root  Migration  with  ADA,  F-14A 


Coupling  numerator.  The  lead  coefficient  (gain  term)  for  the  F-14A 
coupling  numerator,  *8»  ^or  Practical  purposes,  zero. 
Thus  closure  of  either  loop  (pitch  or  roll)  has  essentially  no  influence 
on  the  other. 

Summary .  The  F-14A  exhibits  somewhat  better  open-loop  static  and 
dynamic  stability  than  does  the  F-4J.  It  also  has  considerably  less 
lateral-longitudinal  static  cross-coupling  in  asymmetric  flight, 
and  hence  little  change  or  degradation  in  dynamic  parameters  with 
sideslip.  The  static  aerodynamic  cross-coupling  derivative  values  at 
80  ■  4  deg  and  15  <  a  <  30  deg  are  presented  in  Figure  41.  Comparison 


Figure  41.  F-14A  Static  Aerodynamic  Cross-Coupling 

Derivatives 
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of  these  values  with  those  for  the  F-4J  (Figure  22)  shows  a  drastic 
difference  between  the  two  aircrafts 

The  most  significant  maneuver-limiting  factor  for  the  F-14A  appears 
to  be  roll  reversal  associated  with  movement  of  one  root  of  the  roll 
numerator,  into  the  RHP  as  AOA  is  increased.  The  F-14A  exhibits  a 

more  gradual  shift  of  this  root  with  AOA  than  does  the  F-4J. 
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SECTION  III 


AERODTNAHIC  CROSS-CODFLING  INFLUENCES 


The  previous  section  showed  a  significant  migration  of  denominator 
and  numerator  roots  for  the  F-4J  due  to  aerodynamic  cross-coupling  at 
30  *  5.5  deg.  In  this  section  literal  analytic  expressions  are  devel¬ 
oped  to  examine  in  greater  detail  the  influences  of  specific  cross¬ 
coupling  coefficients  on  the  airframe  characteristic  (open-loop)  modes 
and  the  key  piloted  control  numerators  for  roll  due  to  lateral  stick  and 
pitch  due  to  stabllator  deflections. 

With  lateral-longitudinal  cross-coupling,  at  least  5  DOF  equations 
of  motion  (e.g.,  Figure  13)  must  be  used  for  analysis.  As  a  result,  the 
wealth  of  existing  knowledge  on  the  dynamics  of  vehicles  in  three  (or 
fewer)  degrees  of  freedom  is  not  adequate,  and  qualitative  insight  into 
the  motions  of  the  aircraft  is  sacrificed. 

The  Important  cross-coupling  terms  in  the  equations  of  motion  iden¬ 
tified  empirically  by  Porter  and  Loomis  (Reference  29)  and  Johnston  and 
Hogge  (Reference  5)  are  both  aerodynamic  and  kinematic.  Effects  of 
cross-coupling  on  the  characteristic  aircraft  modes  were  evaluated  by 
Hamel  (Reference  16)  by  applying  servo  analysis  techniques  to  a  5  DOF 
aircraft  mathematical  model. 

The  approach  applied  next  will  follow  and  expand  upon  that  utilized 
by  Hamel. 

A.  CHARACTERISTICS  MODES 

Analysis  of  an  uncoupled  (3  =  0)  airframe  is  relatively  straight¬ 
forward.  Approximations  to  the  classical  3  DOF  lateral  or  longitudinal 
expressions  can  be  developed  for  a  wide  range  of  airframe  and  flight 
conditions  (e.g.,  Reference  20). 

However,  sideslip  characteristically  introduces  cross-coupling,  both 
kinematic  and  aerodynamic.  Then  the  only  clear  simplification  one  can 
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apply  Is  the  assumption  of  constant  forward  velocity,  thereby  removing 
the  velocity  equation  and  associated  terms  (see  Figure  13).  This  can 
be  justified  since  forward  velocity  changes  comparatively  slowly.  No 
significant  coupling  occurs  due  to  velocity,  and  its  effect  is  primarily 
on  the  low-frequency  characteristics  (i.e.,  the  phugold  mode).  The 
three  kinematic  equations  of  Figure  13  can  be  resolved  and  substituted 
in  the  Z,  M,  and  Y  equations  to  reduce  the  number  of  equations  to  five. 

Hamel  takes  f ive-degree-of-freedom  equations  of  motion  "using  body 
axes  which  have  the  x-axis  initially  aligned  with  the  flight  path"; 
from  the  context  we  interpret  that  for  initial  sideslip,  the  x-axis  la 
aligned  with  the  projection  of  the  flight  path  vector  onto  the  plane  of 
symmetry.  In  expanding  the  characteristic  determinant  he  separates  out 
the  effects  of  the  aerodynamic  coupling  terms  N„  and  Mg  (neglect¬ 
ing  any  other  coupling  between  the  longitudinal  set  and  the  lateral- 
directional  set)  thus: 

A  "  AlongAlat  +  K<MAcoupl 
The  uncoupled  (0O  £  0)  denominator  terms  are 

Along  •  a2  +  2cspwaps  +  «2p 

hut  -  (s  +  1/Tb)(s  +  l/TR)(s2  +  2;du>ds  +  uj) 

Hamel  defines  the  classical  longitudinal  (short  period)  and  lateral- 
directional  (spiral,  roll  subsidence,  dutch  roll)  by  approximate  factors 
(Reference  20).  The  coupling  component,  K(0o)Acou^,  can  be  manipulated 
to  represent  a  sum  of  two  terms: 

K(eo)ACOupl  “  *5»o»(«  +  1/T0)(a2  +  2c0u>0s  +  u>2) 

+  M0N'Bo(s  +  1/TiXs  +  1/T2) 

The  approximate  factors  for  these  coupling  terms  are  reproduced  (in  the 
notation  of  this  report)  in  Table  3. 
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TABLE  3.  COUPLING  NUMERATOR  FACTORS 


1/T0  -  -Mq 

2;0u>0  -  -yp  -  n;  + 

(i)2  -  Ng  +  YgNj  -  (N5/Xi)(tg  +  Ypi;) 

1/Ti  -  -tj,  +  (  jC^/Na)(Np  -  g/U0) 

i/x2  -  ti8/[u0(*£/nsKh;  -  £;)] 


The  roots  of  the  coupling  component  are  evaluated  by  applying  servo 
analysis  techniques  for  obtaining  the  roots  of  a  characteristic  equa¬ 
tion,  setting  it  into  the  form 

1  +  Mg(Na/*aXs  +  1/Tl)(»  +  _  0 

s(s  +  1/T0)[s2  +  2CpU0s  +  u)£] 

Figure  42  is  a  root  locus  plot  for  this  equation.  In  this  Instance 
iQ  >  1  and  <  0.  The  "closed-loop"  roots  for  this  case  are 

KCMAcoupl  "  <£aM8  +  -W(s  +  .47) 

X  [s2  +  2(-.48)(1.87)s  +  (1.87)2J 


This  same  technique  can  be  applied  on  a  larger  scale  to  find  the 
roots  for  the  coupled  charactf  '.atic  modes,  i.e.,  the  roots  of  the  rela¬ 
tionship 


B° )4c°Up| 
^long^lat 
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Figure  42.  Migration  of  Roots  of  K(B0)ACOUpi  with 
Increasing  Mg  (<*0  -  21  deg,  B0  ■  5.5  deg) 

Figure  43  is  a  root  locus  plot  for  this  expression  showing  the  migration 
of  ta8p  and  iod  with  increasing  t  Note  that  l/Tg  and  l/Tg  are 
essentially  cancelled  by  zeros  obtained  from  the  approximate  factors  and 
the  preceding  "closure." 

The  actual  6  DOF  roots  of  the  F-4J  aQ  *  21  deg,  ■  5.5  deg  case 
are  shown  by  the  dashed  lines  and  &  symbols  in  Figure  43.  There  is 
fair  agreement  between  the  simplified  and  exact  loci  for  the  dutch  roll 
but  not  for  the  other  modes.  The  problem,  it  appears,  lies  with  the 
location  of  the  first-order  coupling  zeroB.  If  these  were  to  lie  closer 
to  the  origin,  the  uisp  and  Wj  loci  would  be  rotated  counterclockwise  and 
1/Ta  and  1/TR  would  be  driven  to  the  right.  It  is  not  surprising  that 
the  approximate  factors  employed  here  do  not  precisely  predict  the  F-4J 
root  shift  since  they  were  developed  in  Reference  16  for  an  airframe 
configuration  considerably  different  from  the  F-4J.  However,  they  do 
demonstrate  the  trends  and  show  specifically  that: 
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A  Actual  6  DOF  Roots 


Figure  43.  Migration  of  Roots  of  Characteristic 
Modes  as  Cross-Coupling  Increases; 
aQ  -  21  deg,  0O  -  5.5  deg 


•  The  location  of  the  roots  of  K(B0)Ac0upi  (Fig¬ 
ure  42)  and  hence  the  direct  ton  of  the  root  shift 
in  Figure  43,  is  controlled  by”  the  sign  and  mag¬ 
nitude  of  Mp(Ng/jfa). 

•  T^e  magnitude  of  root  shift  in  Figure  43  is  con¬ 
trolled  by  the  sign  and  magnitude  of  £  30. 

B.  N$otk  NUMERATOR 

The  first  step  in  identifying  the  Influences  of  cross-coupling  terms 
on  the  key  numerator  for  lateral  piloting  control,  N$  require  elimi¬ 
nation  of  all  insignificant  terms.  This  is  accomplished  by  a  series 
of  simplifying  steps.  The  6  DOF  characteristic  poles  and  the  zeros 
for  the  numerator  ^  ere  shown  in  Figure  44,  again  for  the  F-4J  at 
«0  "  21  deg,  8q  ■  5.5  deg.  Figure  44a  represents  a  case  for  which  the 
cross-coupling  terms  are  zero.  At  this  high  AOA  the  uncoupled  numerator 
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Figure  44,  N$  stlc/4  Pole-Zero  Locutions,  Coupled  vs.  Uncoupled, 
F-4J,  aQ  -  21  deg,  ■  5.5  deg 


consists  of  two  real  roots,  1/T^  and  l/T^*  non-minimum  phase  and 
minimum  phase,  respectively;  and  two  complex  pairs,  and  u>£p,  which 
exactly  cancel  the  uncoupled  longitudinal  poles.  With  croBs-coupling 
(Figure  44b)  the  first-order  zeros  couple  into  a  complex  pair  located 
near  the  divergent  dutch  roll  mode.  In  addition,  the  numerator  counter¬ 
part  to  map,  4p,  has  moved  apart  from  the  short-period  poles.  The 
effects  of  coupling  on  the  phugold  numerator  term,  u>pi  are  negllble  (see 
Figure  44b).  This  indicates  that  the  velocity  equation  of  Figure  13  may 
be  eliminated.  The  phugold  then  reduces  to  a  first-order  pole-zero  pair 
at  the  origin,  but  all  other  roots  are  affected  only  slightly  (Figure 
44c).  However,  even  with  the  Vy  terms  removed  from  Figure  13,  and 
recognizing  that  the  equations  for  the  Euler  angles  ^  and  8  contribute 
no  significant  kinematic  coupling  and  can  be  deleted,  a  rather  unwieldy 
aix-by-slx  matrix  still  remains. 

The  matrix  can  be  further  reduced  by  use  of  time  vectors  evaluated 
for  the  particular  mode  of  interest.  Inspection  of  Figure  44b  shows 
that  the  ^  zeros  of  greatest  concern  are  the  complex  pair  which  lie 
in  the  right  half-plane,  near  the  already  divergent  dutch  roll  poles. 
Therefore  further  analysis  of  n|  will  concentrate  on  this  complex 
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pair,  and  on  the  effects  of  cross-coupling  in  driving  them  off  the  real 
axis,  where  Figure  44a  places  the  corresponding  zeros  of  the  uncoupled 

response* 

For  denominator  roots  the  relative  magnitude  and  phase  of  each 
motion  variable  in  the  mode  corresponding  to  any  particular  real  pole 
or  pair  of  complex  poles  can  be  determined.  The  result  can  be  expressed 
numerically  as  an  eigenvector  or  graphically  as  phasors.  In  order  to 
determine  such  characteristics  for  transfer  function  numerator  roots,  a 
pole  is  driven  into  each  zero  to  be  examined.  The  poles  are  driven  Into 
the  zeros  by  closing  a  tight  feedback  loop,  as  shown  in  Figure  45.  In 
Figure  44c  the  real  poles  are  driven  into  the  right  half-plane  pair  of 
complex  zeros.  Actually  the  (p/sJSg^  transfer  function  is  utilized 
here  rather  than  <|>/6stlc  since  the  two  are  almost  identical  (compare 
Figures  44c  and  45). 

Figure  46  presents  the  resulting  phasors  for  the  closed-loop  roll- 
0 

spiral  mode,  oigg,  which  approximates  the  numerator  "mode"  w^.  Using 
these  phasors,  vector  polygons  have  been  plotted  for  this  mode  (see 
Figure  46)  from  the  equations  of  motion,  Figure  13.  Minor  vectors 


Figure  45.  Closure  of  p/s  ♦  6gt^  to  Obtain  Time  Vectors, 
F-4J ,  aQ  «  21  deg,  -  5.5  deg 
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correspond  to  negligible  terms;  small  polygons  indicate  little  contri¬ 
bution  to  the  modal  response.  Examination  of  these  vector  polygons 
leads  to  elimination  of  a  number  of  terms  from  the  equations  of  Fig¬ 
ure  13.  Specifically,  the  kinematic  equations  are  no  longer  needed,  and 
so  the  matrix  for  approximation  of  the  coupled  numerator  reduces 

A 

to  the  form  shown  in  Figure  47.  The  derivatives  Yy,  N$  and  Zw  are 
also  minimal  for  this  particular  airframe  and  flight  condition.  How¬ 
ever,  under  other  circumstances  they  can  be  quite  large,  so  they  have 
been  retained  to  obtain  a  more  general  expression  for  . 

8uiv 

The  determinant  of  this  matrix  may  be  developed  as  a  combination  of 
uncoupled  longitudinal  and  lateral-directional  terms  plus  cross-coupling 
terms,  as  follows: 

Nttk  :  *  «stk  {t82  '  +  Mq  +  M&)a  ‘  (M«  +  zwMq)3 


This  is  of  the  form 

"U  ;  coupled  +  **  +  "■ 

where  (^^stk)unC0upled  consist9  of  two  second-order  terms:  a  longitu¬ 
dinal  term  equivalent  to  the  conventional  short-period  (denominator) 
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Figure  47.  Matrix  for  Approximation  of  Coupled  Numerator  N$ 

approximation;  and  a  lateral-directional  term  which  is  the  3  DOF  N$  ^ 

(with  Y$  =0)  (Reference  20).  The  coupling  "gain"  terms  are: 

8  CK 


_u  t  -  N|58tk _ 

jn -  cob  «0 


and 


stk 


Kn  -  MgNJ  coa  o0 


Servo  analysis  techniques  can  now  be  utilized  to  examine  the  influ¬ 
ence  of  each  of  the  cross-coupling  derivatives  by  writing  the  numerator 
in  characteristic-equation  form,  i.e., 


uncoupled 


+  *x 


+  % 


0 


and  examining  the  migration  of  the  roots  as  each  of  the  "gains"  is 
varied. 


Setting  the  coupling  term  KN  to  zero,  the  characteristic  equation 
may  be  written  as: 


1 


+ 


^ _ 

uncoupled 


0 
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Figure  48  is  a  system  survey  plot  showing  the  effect  of  varying 

r  $ 

positive  *  a  on  the  uncoupled  numerator  Ng  The  upper  plot  is  a  root 
locus  and  the  lower  plot  a  Bode-siggie  (see  Section  II)  shoving  the 
influence  of  increasing  (negative)  gain.  These  indicate  that  the  com¬ 
plex  zero  associated  with  the  longitudinal  short  period,  ^gp,  is  driven 
to  higher  frequencies  parallel  to  the  ju»-axie;  and  the  conventional  roll 
zeros,  1/T<f>i  and  l/T^,  are  driven  apart.  The  influence  of  £  £  is  found 
to  be  negligible  in  this  case,  but  this  is  due  to  a  very  small  value  of 
the  ratio  (at  the  flight  condition  being  evaluated)  which 
is  a  factor  in  the  "gain"  term. 

The  prime  contributor  to  cross-coupling  influence  on  the  F-4J  roll 
numerator  is  %  ■  MgN£  cos  a0  as  shown  in  Figure  49,  which  is  a  survey 
for  the  equation 


where  primes  have  been  used  to  denote  that  Kg  has  already  been  ac¬ 
counted  for,  and  double  primes  indicate  that  both  and  KN  have  been 

accounted  for.  KN  produces  a  locus  opposite  in  direction  to  that  of 
K ,  since  both  are  pure  gains  but  have  opposite  sign.  The  root  loca¬ 
tion  1b  much  more  sensitive  to  changes  In  K^,  as  the  relative  Bode 
amplitudes  in  Figures  48  and  49  show. 

The  "gains"  K^,  and  KN  are  functions  of  N£,  Mg,  Ng  and 

*^stk*  *£o  an<*  No  have  been  stressed  in  the  foregoing,  it;  should 

be  obvious  that  the  influence  of  changes  in  magnitude  of  any  of  these 
derivatives  can  be  equally  represented  by  Figures  48  and  49.  Changes 
in  sign  of  any  of  the  derivatives  would  reverse  the  direction  of  root 
migration  in  either  case. 

The  results  shown  here  represent  the  usual  derivative  signs  that 
prevail  as  AOA  increases  toward  stall,  l.e.,  £  ^  positive,  N'  negative, 
and  Mp  negative.  However,  Mg  depends  upon  several  configuration  details 
and  can  be  either  positive  or  negative  (Reference  30).  For  the  above 
signs,  uncoordinated  (8  j1  0)  flight  resultB  in  reduced  dutch  roll 
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stability  (Figure  43  or  44)  and  a  corresponding  shift  of  the  roll 

numerator  zeros.  These  shifts  portend  increasing  closed-loop  roll 

control  problems  with  sideslip  (Figure  45). 

The  accuracy  of  the  5  DOF  N^stjc  approximation  is  evident  from 

Table  4:  the  coupled,  6  DOF  numerator  has  been  matched  almost  identi¬ 
cally.  With  this  validation,  of  the  approximation,  the  following  simple 
coupled  approximate  factors  were  developed  for  this  flight  condition  (it 
must  be  stressed  that  these  literal  approximate  factors  apply  to  this 

aircraft  and  flight  condition  only, and  may  not  be  accurate  for  other 
situations).  For  the  general  form 


N^atk  "  t®2  +  +  *4Ka2  +  24sp“ap8  +  “Jp) 

2^$  -  -(Nr  +  V  “  cos  <*0 

2  •  *  - 

ui^  -  cos  a0  +  /N£Mg  cos  a0 


2^p^p  “  -<Mq  +  Mi)  +  /N^Mg  cos  aQ 
<4p  "  ~Ma  +  zwMq  +  /NaM3  cos  aQ 

Table  4  includes  the  values  given  by  these  approximate  factors, 
indicating  that  they  also  are  accurate  and  therefore,  encompass  the  key 
derivatives. 

These  literal  factors  clearly  indicate  that  longitudinal-lateral 
cross-coupling  in  N$  ^  primarily  through  the  product  N^Mg.  If 

either  derivative  is  zero,  the  crosB-coupling  terms  disappear  and  the 
approximate  factors  become  similar  to  the  literal  expressions  for  un¬ 
coupled  lateral  dutrh  roll  and  longitudinal  short-period  (Reference  20). 
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TABLE  4 


N$6tlt  NUMERATOR  APPROXIMATIONS 


NUMERATOR 

«♦ 

— 

4  8  P 

(ii“ 

8p 

6  DOF  coupled 

-.531 

.909 

.778 

1.361 

5  DOF  (approximate) 
coupled 

. 

Cn 

LO 

U> 

.892 

.768 

1.370 

Based  on  literal 
approximate  factors 

-.574 

.940 

.667 

1.544 

In  addition,  for  non-zero  N£  and  Mg  the  aircraft  coupling  is  aensitive 
to  sideslip  since 


,  32N'  .  9sb 


so  that  Figure  49  may  be  viewed  as  a  locus  of  root  migration  with  side¬ 
slip  for  constant  Cnga* 

From  the  standpoint  of  vehicle  design,  any  airframe  changes  which 
modify  ilther  of  these  cross-coupling  derivatives  can  be  expected  to 
have  a  significant  influence  on  high  AOA  flying  qualities  and  departure 
characteristics. 

C.  N^tab  NUMERATOR 

Analysts  of  the  influence  of  aerodynamic  cross-coupling  on  the 

0 

numerator  for  longitudinal  control,  N6sta^»  U8ea  an  approach  identical 
to  that  taken  with  .  The  first  step  is  to  make  the  development 

tractable  by  eliminating  insignificant  terms  in  the  nine-by-nine  matrix. 

It  wao  shown  in  Subsection  ll.A.3.b  (Figure  19)  that  constraining 
velocity  does  little  to  change  the  zeros  influenced  by  cross-coupling, 
i.e.,  I/Tq^,  ^tR»  and  ui<j.  Additionally,  we  find  from  Figure  50  that  the 
equations  for  the  Euler  angles  and  6  can  again  be  eliminated.  This 
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Figure  50.  Time  Vectors  for  Closed-Loop  Mode  1/TS 
(Approximating  1/T03),  F-4J; 
a  «  21  deg,  B  -  5,5  deg 
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still  leaves  a  six-by-six  matrix  which  Is  too  large  and  complex  to  work 
with  comfortably. 

Applications  of  time  vectors,  as  was  done  with  N^atjc»  *s  accom¬ 
plished  by  closing  a  pitch  rate  loop  with  a  feedback  element  of  K/s 
(to  obtain  body-axis  9).  Modal  response  coefficients  are  then  calcu¬ 
lated  for  the  closed-loop  pole  driven  into  I/Tq^,  and  the  time  vectors 
for  each  equation  constructed  as  shown  in  Figure  50.  Except  f^r  the 
at  equation,  few  terms  are  small  enough  to  be  eliminated  by  inspection. 
Further, it  seems  unwise  to  restrict  the  general  applicability  of  the 
results  by  eliminating  the  a  equation.  The  main  benefit  of  the  time 
vectors  is  to  demonstrate  that  the  modal  response  for  a  root  at  1/Tq ^ 
contains  considerably  more  lateral  than  longitudinal  motion. 

Thus,  it  is  necessary  to  expand  the  complete  six-by-six  matrix  of 
Figure  51  in  literal  form,  evaluate  each  of  the  terms,  and  eliminate 
products  of  small  numbers.  The  resulting  polynomial  is  then  separated 
into  uncoupled  and  coupled  terms  to  form  the  approximation: 
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in  the 
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of  Reference  18, 
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Figure  51 .  Matrix  for  Approximation  of  Coupled  N5 

(4>0  Mo-0,  0  -  a,  80  Small)  8tab 


The  terms  not  containing  cross-coupling  derivatives  are  seen  to 
separate  into  a  firat-order  and  quartic  function  of  s.  The  first- 
order  is  the  conventional  literal  approximation  for  the  longitudinal 
zero  1/Tg2>  i.e.,  -Zw.  The  quartic  contains  all  lateral-directional 
derivatives  and,  with  the  exception  of  the  sin  aQ  terms,  is  the  same  as 
the  conventional  body-axis  uncoupled  lateral  characteristic  equation 
(Reference  20).  This  polynomial  thus  factors  into  the  parameters 
1/TR,  and  1/TS. 

The  first  cross-coupling  term  involves  as  the  "gain,"  a  free  s, 
and  a  quadratic  containing  only  yawing  derivatives.  The  second  cross¬ 
coupling  expression  haB  N£  as  a  "gain,"  a  free  s,  and  a  quadratic  con¬ 
taining  only  rolling  derivatives. 

The  kinematic  terms  Zp  and  Zr  were  negligible  for  hut  here 

dictate  the  influence  of  the  aerodynamic  cross-coupling.  As  a  result, 
coupling  ir.  the  0  *  ^8tab  '‘Operator  18  un  explicit  function  of  sideslip. 
It  is  interesting  to  note  that  Mg  does  not  appear  in  the  approximation; 
a  quadratic  coupling  term  involving  MBz68tub^T0  was  found  to  be  negli¬ 
gible. 


The  Influence  of  the  individual  cross-coupling  terms  can  again  be 
observed  by  servo  analysis  methods,  setting 


(»L.J 


stab  coupled 


- 


ata^  uncoupled 


(“lublc. 


+  (n5  ) 

^  °stab^ 


K 


and  performing  system  surveys.  The  survey  for 


(N§stab^ 

1  +  — g - =2 -  -  0 

^8tab^unc0Upiec| 


is  shown  in  Figure  52.  The  "poles"  1/T®  and  1 /Tq^  couple  to  create  a 
second-order  closed-loop  root,  Wqs,  which  is  on  the  ju-axis.  The  dutch 
roll  "pole"  is  moved  to  the  left  half-plane,  while  1/T^  is  moved  slight¬ 
ly  to  the  left. 

Figure  53  shows  the  additional  movement  when  the  cross-coupling  due 
to  Ng  is  included.  The  latter  moves  things  further  in  the  same  direc¬ 
tion.  However,  it  can  be  determined  that  the  first-order  coupled  zero 

in  the  left  half-plane  comes  from  the  roll  subsidence  mode  and  the  zero 

farthest  into  the  right  half-plane  (previously  called  1  /Tq^)  comes  from 
a  coupling  of  the  spiral  mode  and  the.  conventionale  1/Tq2»  Since  the 
time  vectors  indicated  this  mode  is  dominated  by  lateral-directional 
motion  it  is  now  labeled  (l/Tg)  .  The  prime  notation  in  Figures  52 
and  53  is  used  to  keep  track  of  the  number  of  "loops"  being  closed 

and  the  influence  of  each  "closure."  From  here  on  the  prime  notation 
will  be  dropped  and  the  lateral-longitudinal  coupled  zeros  identified 
as  w®,  1/tr*  1/Tg>  etc. 

Since  each  of  the  croBs-coupling  terms  has  similar  influence  on 

the  migration  of  zeros,  no  insight  la  lost  by  combining  the  two  coupling 
terras  via  the  concept  of  the  previously  noted  "dynamic"  terra  Then 
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The  system  survey  for  the  combined  coupling  influence  on  Nfistab  ia  shown 
in  Figure  54;  the  coupled  roots,  of  course,  are  the  same. 

Finally,  it  should  be  noted  that  aerodynamic  cross-coupling  in  the 
0 

numerator  ^6sta^  ifl  a  direct  function  of  sideslip  squared.  Since  is 
proportional  to  C^^(e0)  and  Na  is  proportional  to  cnpa(M.  the  "gain" 
term  is 

*«dyn0°  ~  (%0Bo  COs  a°  +  °VB°  sin  ao^ 

~  (c£g0  cos  ao  +  cnga  8i"  ao)Bo 

Thus  the  zero  1/t|  is  drastically  influenced  by  static  aerodynamic 
cross-coupling  and  magnitude  of  sideslip. 

D.  SUMMARY 

Literal  analytic  expressions  have  been  developed  to  show  the  influ¬ 
ence  of  static  aerodynamic  cross-coupling  derivatives  and  asymmetric 
flight  on  denominator  and  numerator  roots.  The  analysis  has  shown  that 
the  shift  in  denominator  roots  depends  upon  all  three  static  cross¬ 
coupling  derivatives,  NJ,  and  Mg.  The  direction  of  root  migration 

is  controlled  by  the  sign  and  magnitude  of  MgNa/  The  magnitude  of 

root  migration  is  determined  by  ^a80. 


^6stab^coupled 
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Figure  54.  System  Survey  for  Effect  of  Combined  Cross-Coupling 
°n  N§stab;  F-4J;  ac  -  21  deg,  $0  “  5.5  deg 
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For  the  F-4J  the.  principal  croBB-coupling  Influence  on  the  roll 
numerator,  is  due  to  the  term  N^Mg  cob  aQ.  The  literal  approxi¬ 
mate  factors  which  describe  the  coupled  numerator  roots  are 


25^ 

4 

*4p4p 


-(n;  +  Yv)  -  /N'MB  cob  n0 
Ng  cos  aQ  +  /NqMb  cos  oq 
“(Mq  +  M*)  +  cos  a0 


“sp  *  "Ma  +  ZyMq  +  ooa~a0 

ft 

The  pitch  numerator,  N68tajj»  Influenced  by  a  more  complicated 
combination  of  kinematic  and  aerodynamic  terms.  Existence  of  a  first- 
order  RHP  zero  depends  upon 


—  X  gN0  —  Yv(  tC  y  cos  &q  “  p  sin  ct0) 
- *£  cos  a0'  +  N’  sin  a;  -  <  0 


The  above  approximations  may  not  apply  to  configurations  significantly 
different  from  the  F-4J. 
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SECTION  IV 


CONFIGURATIONS  SELECTED  FOR  SIMULATION 


The  seven  aerodynamic  derivatives  which  have  been  shown  to  dominate 
open-  and  closed-loop  dynamic  characteristics  at  high  AOA  are  summarized 
in  Table  5.  These  are  comprised  of  six  static  stability  and  cross- 
coupling  derivatives  and  one  dynamic  damping  derivative.  Together  they 
limit  the  useful  flight  envelope  via  roll  reversal  and  wing  rock  tenden¬ 
cies,  and  they  limit  the  safe  flight  envelope  via  departure  suscepti¬ 
bility.  One  goal  of  the  piloted  simulation  was  to  demonstrate  this 
premise. 


TABLE  5 


KEY  MANEUVER-LIMITING  PARAMETERS 


Negative  Nsa  or  N$D 

Roll  reversal 

Positive  Mg 

Pitch  Up 

£b  *a.  /p  ' 

l  Wing  rock 

Ng,  N'a 

> 

\  Nose  slice 

HS  ) 

\  Roll  divergence 

This  section  presents  the  vehicle  aerodynamic  configurations  se¬ 
lected  for  the  piloted  simulation  and  the  predicted  high  AOA  dynamic 
characteristics. 

Four  "configurations"  were  obtained  by  altering  the  key  aerodynamic 
coefficients  of  a  baseline.  The  three  altered  configurations  incor¬ 
porated  changes  in,  respectively,  C^p(a),  Cjj(a,B),  and  Cn(a,fJ)  combined 
with  a  change  in  Cmg,  to  provide  specific  alteration  of  key  open-  and 
closed-loop  parameters  discussed  in  the  preceding  two  sections.  The 
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attempt  was  to  Influence  the  denominator  only,  the  numerator  only,  and 
the  combined  denominator  and  numerator,  thereby  to  provide  distinct 
variations  in  high  AOA  departure  susceptibility,  onset  warning,  and 
motion  severity.  The  influence  of  artificially  altering  aerodynamic 
coefficients  via  a  lateral-directional  augmentation  system  was  also 
investigated  with  two  of  the  configurations. 

A.  AERODYNAMIC  CONFIGURATIONS 

1 ■  Configuration  A 

The  baseline  case  is  iden  'led  as  Configuration  A.  This  F-4J- 
based  configuration  has  been  thoroughly  reviewed  in  the  preceding  sec¬ 
tions.  The  key  open-  and  closed-loop  dynamic  parameters  are  summarized 
in  Figure  55  for  comparison  with  the  other  configurations  discussed 
subsequently.  Dutch  roll  and  numerator  root  migrations  are  shown  for 
increasing  a  at  zero  0  (solid  lines)  and  at  5.5  deg  8  (dashed  lines). 
The  solid  triangle  (A)  indicates  the  AOA  above  which  Cn^  becomes 
negative  at  zero  0.  Recalling  the  discussion  in  Section  II  about  the 
shortcomings  of  linearized  frozen  point  analysis,  the  important  aspects 
to  note  here  are: 

The  low  dutch  roll  damping  at  AOA  below 
(wing  rock). 

2 

or  1  /T jj,  negative  at  AOA  above  18  deg  or  at 
small  8  (roll  reversal). 

1/T*.  and  1/Tg  both  large  and  negative  when  6^0 
(nose  slice  parameters). 

e  Dutch  roll  and  roll  numerator  roots  sensitive  to 
sideslip. 

The  plots  in  the  lower  left  corner  indicate  the  magnitude  of  the 
aerodynamic  coupling  terms  evaluated  at  5.5  deg  8.  In  the  AOA  range 
from  20  to  30  deg,  is  of  approximately  the  same  magnitude  as  3?|  but 
is  of  opposite  sign.  In  this  same  AOA  range,  is  of  about  the  same 
magnitude  as  Njj  and  of  the  same  sign.  For  this  relatively  mild  sideslip 
condition  (5.5  deg),  the  static  cross-coupling  derivatives  are  of  equal 
magnitude  to  the  basic  static  derivatives. 


e 
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2.  Configuration  B 


For  Configuration  B  the  roll  rate  damping  coefficient,  C^,  was  made 
large,  negative  and  essentially  invariant  over  the  range  10  <  a  <  35  deg 
(see  Cs  ,  Figure  56).  The  object  of  this  change  is  to  increase  the 

*  4 

dutch  roll  damping  and  hence  reduce  the  wing  rock  tendency.  Since  all 
other  parameters  remain  unchanged,  this  should  not  Alter  the  basic  roll 
reversal  tendency  or  the  nose  slice  departure  characteristics;  but  It 
should  eliminate  or  reduce  wing  rock  as  a  warning  of  impending  depar¬ 
ture.  However,  note  again  in  Figure  57  the  sensitivity  of  the  dutch 
roll  mode  to  sideslip.  At  8  -  5.5  deg  there  is  little  difference 
between  Configurations  A  (Figure  55)  and  B  (Figure  57),  which  indicates 
that  the  Influence  of  «Cp  is  small  with  respect  to  that  of  jCjj. 


-DDK 

Clp 

PER  RRB 

-can 

-□as 

-o.za 

-0.3E 

-□.3D 


Figure  56.  Change  in  Cg  for  Configuration  B 
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3.  Configuration  C 


For  Configuration  C  the  roll  due  to  sideslip  coefficient,  C$g,  was 
altered  (see  Figure  58)  to  approximate  that  of  the  F-14A  over  the  range 
15  <  a  <  45  deg,  taking  into  account  the  difference  in  roll  inertias  of 
the  two  aircraft,  i.e., 


'&F-14 


r*F-4 

1j«F-14 


Figure  59  presents  the  parameter  plots  for  Configuration  C.  The 
upper  left  figure  shows  that  the  change  Introduced  in  Cg^  and  Cga  has 
shifted  the  dutch  roll  mode  toward  the  left.  It  remains  in  the  left 
half-plane  for  8-0  and  only  bulges  into  the  right  half-plane  over  a 
small  a  range  at  8  -  5.5  deg.  A  major  difference  in  this  configuration 
occurs  in  the  roll  subsidence  and  spiral  modes,  which  in  the  previous 


Figure  58.  Change  in  for  Configuration  C 
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cases  have  remained  on  the  real  axis  and  therefore  have  not  been  shown. 
In  this  case,  1/TR  and  1/T8  couple  to  form  a  "lateral  phugoid”  mode, 
uf)R,  as  AOA  increases.  At  36  deg  AOA  the  frequency  of  this  mode  is 
approximately  twice  that  of  the  dutch  roll,  with  approximately  the  same 
damping.  Comparison  of  Figures  35,  39,  and  59  show  this  to  be  a  valid 
representation  of  F-14A  root  locations,  but  mode  Identity  may  be  differ¬ 
ent. 


Comparison  indicates  that  increasing  6  shifts  the  u>gR  locus  toward 
the  right  half-plane.  However,  both  the  dutch  roll  and  lateral  phugoid 
modes  are  now  relatively  unaffected  by  3. 

The  upper  right-hand  plot  for  is  essentially  the  same  for  Con¬ 
figuration  C  as  that  for  Configurations  A  and  B.  While  this  configura¬ 
tion  has  better  open- loop  stability,  the  closed-loop  characteristics 
should  be  little  different  from  the  previous  two  configurations  because 
the  roll  numerator  roots  are  essentially  unchanged. 


Absence  of  the  solid  triangle  symbol  indicates  that 
tive  throughout  the  AOA  range. 


is  posi- 


The  lower  left-hand  figure  shows  that  the  term  (which  derives 
from  the  variation  of  Cj^  with  « )  is  small  and  negative.  Since  this 
configuration  also  has  a  large  increase  in  (see  Figure  58),  the 
aerodynamic  cross-coupling  derivative  is  now  much  smaller  than  £  | 
throughout  the  angle-of-attack  range.  Comparing  the  values  of  £ „  for 
Configurations  A  and  C  with  the  differences  in  denominator  root  migra¬ 
tions  with  8  makes  it  apparent  that  /j  is  the  dominant  purameter  in 
shifting  the  Configuration  A  dutch  roll  mode  into  the  right  half-plane 
when  8^0.  The  lower  right-hand  plot  of  Figure  59  indicates  that  the 
change  in  and  £  jj  has  decreased  1/t|  at  AOA  above  20  deg. 
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4.  Configuration  D 


For  the  fourth  case  the  yaw  due  to  sideslip  coefficient,  Cn^ ,  was 
modified  to  approximate  that  of  the  F-14  as  shown  by  the  dashed  line  In 
Figure  60,  In  addition,  for  this  case  only,  the  sign  of  Cmg  was  re- 
versed  to  provide  positive  pitch  with  sideslip.  As  predicted  by  the 
approximate  factors  of  Section  III,  the  result  Is  an  airframe  quite 
Insensitive  to  sideslip  (Figure  61)  and,  most  pertinently,  having  numer¬ 
ator  roots  (l/T^  and  1/T^)  which  do  not  penetrate  as  far  into  the  right 
half-plane . 

Roll  reversal  is  delayed  until  a  >  20  deg  and 
tive  until  a  *  25  deg. 

aaaao 

o.aoaa 

Chit 

cch  sen 

a.aaiD 

□.no 

-□.□□in 

-a, Damn 

-0.DD3D 

-□.□m-a 

-□.□□so 

Figure  60.  Change  in  Cn^  for  Configuration  D 


Cn^  remains  posi- 


5.  Anticipated  High  AQA  Maneuver- 
Llaltlng  Characterlatica 


Review  of  the  preceding  four  configurations  should  go  a  long  way 
toward  sorting  out  the  relative  Importance  of  denominator  (open-loop) 
vs.  numerator  (closed-loop)  characteristics  on  high  AOA  departure  sus¬ 
ceptibility  and  severity.  Configurations  A,  B,  and  C  have  nearly  iden¬ 
tical  RHP  numerator  root  locations  and  therefore  should  have  similar 
closed-loop  divergence  tendencies,  although  the  motion  (l.e.,  roll  vs. 
yaw)  may  differ  somewhat  due  to  C$,pl  Cjg,  and  Cjj,a  changes.  On  the  other 
hand,  the  Configuration  D  numerator  roots  do  not  penetrate  so  far  into 
the  RHP  and  are  relatively  insensitive  to  sideslip.  Thus  one  might 
expect  this  vehicle  to  appear  less  susceptible  to  or  have  less  severe 
departure. 

Alternately,  if  departure  susceptibility  or  tendency  should  turn  out 
to  be  strongly  influenced  by  dutch  roll  (open-loop)  root  location,  then 
one  would  expect  Configurations  C  and  D  to  be  similar  Bince  both  are 
near  neutral  stability  and  relatively  insensitive  to  sideslip.  Any 
differences  between  Configurations  A  and  B  would  be  dependent  upon  side¬ 
slip. 

The  dynamic  stability  parameter,  Cna  ,  has  been  referred  to  pre- 

^dyn 

vlously.  The  AOA  at  which  this  parameter  becomes  negative  was  identi¬ 
fied  in  Figures  55,  57,  and  6i.  Figure  62  presents  plots  of 
Cn^  vs.  a  for  the  four  configurations.  BaBed  upon  the  criterion 

of  Cn  >  0.002  for  departure  resistance  (Reference  31),  then  only 

Pdyn 

Configuration  C  should  be  resistant,  Configurations  A  and  B  should  be 
quite  departure  prone,  and  Configuration  D  might  be  mildly  departure 
prone. 

Another  criterion,  Reference  32,  employs  both  open-  and  closed-loop 
parameters,  i.e., 


Cnat1  “  (cn*  cAr/cJU  ] 


cns  C08  « 


■  sln  a 


Figure  62. 


Cna  for  Configurations  Simulated 
»dyn 


For  a  number  of  years  it  has  been  recognized  (Reference  33)  that  unde¬ 
sirable  ratios  lead  to  roll  reversal  and  pilot-aggravated  wing 
rock  (PIO).  However,  In  Reference  32  Weissman  plots  Cn&  vfl’  LCDP  and 
identifies  four  regions  of  increasing  departure  and  spin  susceptibility 
and  severity  (Figure  63)  based  upon  experience  with  a  number  of  airframe 
configurations.  The  loci  of  our  unaugmented  airframe  parameter  values 
over  the  range  of  16  <  a  <  35  deg  are  shown  in  Figure  63;  circles  repre¬ 
sent  Configurations  A  and  B  (since  C*p  has  no  influence  on  either  param¬ 
eter),  squares  represent  Configuration  C,  and  triangles  Configuration 
D.  The  criteria  predict  high  departure/spin  susceptibility  with  strong 
rolling  departures  for  Configurations  A  and  B,  low  spin  susceptibility 
and  mild  rolling  departures  for  Configuration  C,  and  no  departures  for 
Configuration  D.  Thus,  all  regions  of  departure  susceptiblity  and 
severity  are  exercised  with  the  aerodynamics  selected. 

The  Figure  63  prediction  that  Configuration  D  should  be  less  suscep¬ 
tible  to  departure  is  in  agreement  with  the  first  assessment  above  bused 
on  numerator  root  locations. 
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B.  AUGMENTED  OORFIGUBAIIONS 


Lateral-directional  augmentation  was  synthesized  for  Configuration  A 
(the  most  departure  prone)  and  for  Configuration  C  (Cn^  always  posi¬ 
tive).  The  augmented  configurations  are  identified  as  A2  and  C2,  re¬ 
spectively,  while  the  unaugmented  configurations  will  be  referred  to  as 
A^  and  Cj .  The  augmentation  mechanization  is  representative  of  that 
used  in  current  high-performance  fighters  to  improve  both  open-  and 
closed-loop  handling  parameters,  i.e.,  washed-out  yaw  rate  and  lateral 
acceleration  feedbacks  to  rudder  to  increase  dutch  roll  damping;  roll 
rate  feedback  to  differential  stabilator  to  reduce  or  eliminate  wing 
rock  and  augment  dutch  roll  damping;  and  lateral  stick  to  rudder  cross- 
feed  (SRI)  to  reduce  or  eliminate  adverse  aileron  yaw  for  pilot  inputs, 
thereby  providing  a  positive  shift  in  the  closed-loop  parameters  LCDP, 
«?,  and  1/Tx. 


Undesirably,  the  roll  damper  also  produces  adverse  yaw  when  opposing 
uncommanded  roll;  it  may  increase  nose  slice  departure  tendencies  and 
pro-spin  surface  deflections.  A  block  diagram  for  the  augmentation 
mechanization  is  shown  in  Section  V, 


The  root  migration  surveys  for  Configuration  A2  are  shown  in  Fig¬ 
ure  64,  Comparison  with  Figure  55  shows  that  at  0  -  0  the  dutch  roll 
damping  is  increased  considerably;  it  does  not  go  unstable  until  AOA 

exceeds  that  for  Cn„  becoming  negative  ( A  symbol).  Also,  the  AOA  at 
?  &dyn 

which  becomes  negative  has  been  increased  approximately  5  dc g  and  the 

negative  values  of  both  1/T^  and  1/T®  have  been  reduced.  Thus,  within 
the  augmentation  system  authority  limits  and  near  zero  0,  Configura¬ 
tion  A2  stability  is  increased  and  departure  susceptibility  decreased  up 
to  about  23  deg  AOA.  However,  for  B  1*  0  the  augmentation  benefits  are 
overcome  and  the  dynamic  parameters  revert  back  to  air4?rurae-alone 
values. 
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Figure  64.  Key  Parameters  for  Conf iguration  A2 


Root  migration  survey  plots  for  Configuration  C2  are  presented  in 
Figure  65.  Comparison  with  Figure  57  again  shows  Improvement  in  dutch 
roll  stability.  However,  the  most  dramatic  change  is  in  ui|  which  now 
remains  positive  up  to  35  deg  AOA  (but  note  is  small  and  negative 

beyond  31  deg  AOA).  Thus  roll  reversal  problems  are  essentially  elimin¬ 
ated  and  numerator  root  sensitivity  to  sideslip  is  decreased  consider¬ 
ably  —  again  subject  to  augmentation  system  authority  limits. 

The  change  in  departure  and  spin  characteristics  cannot  be  pre¬ 
dicted  for  these  configurations  using  the  static  aerodynamic  parameters 
of  Figure  63  because  the  augmentation  system  influence  is  frequency 
dependent.  But  a  somewhat  similar  plot  can  be  constructed  using  the 
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Figure  654  Key  Parameters  for  Configuration  C2 
2  * 

parameters  and  C^d*  Such  a  plot  la  presented  in  Figure  66  for 

both  unaugmented  and  augmented  versions  of  Configurations  A  and  Ct  The 
2 

intercept  of  the  axis  with  trf^d  has  been  arbitrarily  shifted  to  en¬ 
hance  similarity  for  the  unaugmented  configurations  with  Figure  63. 
Note  that  the  root  plots  for  the  unaugmented  Configurations  Aj  and 
(solid  lines)  are  quite  similar  to  those  of  Figure  63.  The  plots 

(Figure  66)  for  Configurations  A2  and  C2  (dashed  lines)  are  moved  fur¬ 
ther  into  the  upper  right  quadrant,  which  is  identified  as  departure 
resistant  in  Figure  63.  Similar  interpretation  would  forecast  A2  to 
be  departure  resistant  up  to  about  24  deg  AOA  but  to  have  moderate  to 
severe  departure  tendencies  at  higher  AOA,  while  C2  should  have  no 
departure  tendencies. 


Despite  the  direct  relationship  between  Cn  and  it  is  move 


appropriate  to  plot  because  the  changes  in 

of  Cn  and  are  in  closer  agreement. 

8dyn 


dyn 

KXgt 


gn  (or  lack  thereof) 
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Figure  66.  Versus  for  Augmented  and  Unaugmented 

Configurations  A  and  C 


C.  SUMMARY  OF  CONFIGURATIONS  AND 

PREDICTED  HIGH  a  CHARACTERISTICS 

Table  6  summarizes  the  six-conf iguration  matrix  employed  in  the 
piloted  simulation.  For  A^  a  sequence  of  roll  reversal,  wing  rock,  nose 
slice,  and  finally  strong  rolling  departure  with  high  spin  suscepti¬ 
bility  is  predicted  with  increaeing  AOA.  For  A2  the  augmentation  system 
and  crossfeed  are  expected  to  minimize  or  eliminate  the  roll  reversal 
and  wing  rock  predeparture  warnings.  The  predicted  characteristics  are 
nose  slice  followed  by  moderate  rolling  departure.  For  B,  predicted 
characteristics  are  roll  reversal,  nose  slice,  and  strong  rolling  depar¬ 
ture.  Warning  In  the  form  of  wing  rock  should  not  be  present  due  to  the 
large  roll  damping  of  this  configuration.  Note  that  Configurations  A2 
and  B  allow  comparison  of  high  AOA  stall/departure  characteristics  with 
high  aerodynamic  roll  damping  (Configuration  B)  and  with  artificially 
augmented  roll  damping  (Configuration  A2). 

The  unagumented  Configuration  C^  is  predicted  to  exhibit  roll  rever¬ 
sal,  wing  rock,  and  mild  rolling  divergence  with  increasing  AOA  but  to 
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TABLE  6.  CONFIGURATION  MATRIX 


CONFIGU¬ 

RATION 

FCS 

AERODYNAMIC  VARIANT 

PREDICTED 

CHARACTERISTIC 

A1 

Basic 

Basic  F-4J 

Wing  rock 

Roll  reversal 

Nose  slice 

Roil  departure 

A2 

Augmented 

Nose  slice 

Roll  departure 

B 

Basic 

Increased  Co 

P 

Roll  reversal 

Noae  slice 

Roll  departure 

C1 

Basic 

Increased  Co.  ) 

p  [  15  <  a  <  45 
Decreased  C«  ) 

Wing  rock 

Roll  reversal 

Roll  departure 

c2 

Augmented 

Cna.  >  0 
None  ’’dyn 

LCDP  >  0 

D 

Basic 

Increased  Cno  ) 

P  J  0  >  15 

Decreased  Cn  ) 

“a 

Positive  C0g 

Wing  rock 

Roll  reversal 

Pitch  up 

have  no  nose  slice.  Augmentation  is  employed  with  C2  to  Improve  maneu- 
vering  control  and  to  determine  If  it  will  degrade  the  departure  resis¬ 
tance  of  this  configuration.  On  the  basis  of  C„fl  and  (or  LCDP), 

Pdyn  r 

no  departure  tendency  should  be  anticipated  for  C2.  Finally,  D  should 
exhibit  roll  reversal  and  wing  rock  warnings  but  little  lateral-direc¬ 
tional  departure  tendency.  However,  thu  change  in  sign  of  Cg^  to  pro¬ 
vide  positive  pitching  moment  with  sideslip  should  result  in  pitch-up, 
which  would  be  expected  to  aggravate  any  high  AOA  departure  tendency. 
All  other  configurations  have  negative  C^,  should  pitch  down  with 
increasing  sideslip,  and  hence  should  require  more  effort  of  the  pilot 
to  maintain  high  AOA  (longitudinal  stick  cue). 
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SECTION  V 


PILOTED  SIMULATION 


The  purposes  of  the  piloted  simulation  were  to  validate  the  influ¬ 
ences  of  the  key  aerodynamic  coefficients  in  determining  departure 
characteristics,  evaluate  the  Influences  of  varied  maneuver  limiting 
factors  on  high  AOA  maneuvering  control,  and  identify  potential  flying 
qualities  criteria.  While  aircraft  motion  is  highly  desirable  for  such 
validation  and  assessment,  departure  and  post-stall  gyration  (PSG) 
severity  obtained  with  the  6  DOF  analytic  models  indicated  that  a 
moving-base  simulation  would  offer  very  little  benefit.  Noee  slice 
motion  was  so  rapid  that  130  deg/sec  yaw  rate  could  be  reached  within 
two  seconds  after  departure  onset;  peak  angular  accelerations  as  large 
as 

p  ■  500  deg/sec^ 

q  “  200  deg/sec^ 

i  •*  150  deg/sec^ 

were  obtained.  Washout  requirements  to  prevent  hitting  displacement 
stops  would  have  to  be  so  rapid  as  to  negate  motion  benefits.  There¬ 
fore,  a  fixed-base  simulation  was  selected  having  an  unrestricted  out- 
of-cockpit  horizon  display  capable  of  such  rates  and  accelerations.  The 
aerodynamic  model  was  also  expanded  to  accommodate  Inverted  flight 
maneuvers. 

This  section  presents  an  overview  of  the  simulator,  including  dis¬ 
plays,  airframe,  and  flight  control  system  modeling,  cockpit,  and  force 
and  moment  equations. 
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A.  SIMULATOR 


The  simulation  was  performed  in  the  McDonnell  Aircraft  Company's 
fixed-base  Manned  Air  Combat  Simulator  dome,  a  20  ft  hemisphere  identi¬ 
fied  as  MACS-1.  Physical  aspects  of  the  simulation  are  summarized  in 
Figure  67.  The  horizon  and  target  are  projected  on  the  inside  of  the 
hemisphere.  The  cockpit  is  located  at  the  center  of  the  dome.  The  out- 
the-window,  head-up  (HUD)  and  head-down  (HDD)  displays,  and  cockpit 
layou;  are  as  indicated  in  Figure  67 .  Seat  cues  consisted  of  normal 
acceleration  and  buffet  motion  provided  through  an  inflatable  seat 
bladder.  The  frequency  and  amplitude  of  the  buffet  oscillation  varied 
as  the  aircraft  angle  of  attack  increased,  starting  at  about  14  deg. 

A  TV,  projection  of  a  gimballed  model  target  aircraft  provided  a 
maneuvering  tracking  task.  Target  motion  equations  provided  realistic 
maneuvers  while  allowing  the  aircraft  to  fly  well  beyond  the  departure 
region  for  the  test  aircraft.  The  target  was  controlled  by  the  computer 
operator  via  a  special  control  panel. 

B.  HEAD-UP  DISPLAY 

The  head-up  display  (HUD)  provided  Information  on  aircraft  attitude, 
heading,  airspeed,  altitude,  and  trim  conditions.  Figure  68  illustrates 
the  HUD  presentation  as  originally  installed.  Display  of  normal  accel¬ 
eration  provided  a  reference  for  trimming  longitudinal  control  to  1  g 
prior  to  the  start  of  a  run.  Trim  thrust  was  set  by  adjusting  throttles 
until  change  from  trim  thrust  registered  near  zero.  The  Flight  Path 
Angle  Ladder  displayed  both  flight  path  elevation  angle  and  roll  atti¬ 
tude.  Orientation  of  the  total  velocity  vector  was  provided  through 
the  Velocity  Vector  Indicator;  when  the  velocity  vector  was  outside  the 
HUD  field  of  view,  the  Indicator  remained  at  the  edge  of  the  HUD  and 
blinked. 

At  pilot  request,  the  HUD  was  simplified  late  in  the  simulation 
program.  Airspeed,  Heading,  and  Altitude  bar  displays  were  removed. 

The  fixed  reticle  tracking  gunsight  was  a  part  of  the  HUD. 
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Figure  68.  Head-Up  Display  (HUD)  Used  In  Simulation 


C.  FLIGHT  OOtmOL  SYSTEM 

The  manual  flight  control  mechanization  (atick,  pedals,  surface  de¬ 
flections  and  rates,  etc.)  represented  the  basic  F-4  aircraft  (e.g.,  see 
Figure  2). 

A  special  lateral-directional  stability  and  command  augmentation 
system  (SCAS)  was  mechanized  as  shown  in  Figure  69.  The  SCAS  feedback 
gains  are  production  F-4  values.  The  stlck-to-rudder  crossfeed  gain 
and  shaping  were  optimized  for  10  <  a  <  20  deg  and  faded  to  zero  at 
a  <  10  deg  to  prevent  adversely  influencing  the  low  AOA  handling  quali¬ 
ties. 


D.  EQUATIONS  OF  MOTIO* 

Slightly  modified  versions  of  the  standard  McDonnell  Aircraft  Com¬ 
pany  piloted  simulation  equations  of  motion  were  used.  These  are  essen¬ 
tially  the  same  as  those  presented  in  Appendix  I,  Part  III,  i.e.,  moment 
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equations  in  body  centerline  axes,  force  equations  in  wind  axes,  and 
aircraft  orientation  in  standard  Euler  axes.  The  only  change  required 
was  to  expand  the  small  angle  a  and  $  assumptions  and  to  eliminate 

the  tan  0  and  cos"'1’  0  expressions  of  the  <1  equation  which  gave  discon¬ 
tinuity  problems  as  0  approached  90  deg.  This  required  substituting 

the  ^equation  for  4,  i.e., 

4  -  q  -  tan  0(p  cos  a  +  r  sin  a)  +  Zy/CmVij  cos  0) 

in  ■  qVj  cos  0  -  VT  sin  0(p  cos  a  +  r  sin  a)  +  2^/m 

and  using  the  simple  expressions  a  ■  sin-*  (w/V^)  and  0  ■  sin-*  (v/VT) 
for  entry. into  the  aerodynamic  look-up  tables. 

E.  AERODYNAMIC  MODELS 

Since  the  major  interest  in  the  simulation  was  the  analysis  of 

departure  onset,  particular  attention  was  given  to  realistically  model¬ 
ing  aerodynamic  data  for  the  range  of  0  <  a  <  45  deg  and  0  <  ±30  deg. 
The  basic  aerodynamic  coefficients  were  selected  with  the  goal  of  pro¬ 
viding  a  reasonably  accurate  dynamic  model  of  the  F-4J  over  this  a  and  0 
range.  To  allow  greater  aircraft  motion  freedom  typical  of  severe 
departures,  post-stall  gyrations,  and  spins  without  encountering  unreal¬ 
istic  discontinuities  in  the  data,  certain  approximations  were  made  for 
a  >  45  deg,  a  <  0  deg,  and  0  >  ±30  deg,  i.e., 

•  For  a  >  45  deg  the  coefficients  were  faired  to 
110  deg,  then  extended  to  180  deg  by  either  main¬ 
taining  the  value  at  110  deg  or  fading  the  coef¬ 
ficient  to  zero  at  a  ■  180  deg. 

•  The  coefficients  were  assumed  to  be  similar  In 
value  for  positive  and  negative  A0A  so  that  all 
derivatives  were  taken  to  be  even  functions  of  A0A 
(i.e.,  mirror  images  about  a  -  0)  and  the  coeffi¬ 
cients  and  Cm  to  be  odd  functions  of  AOA. 

•  Variation  of  the  coefficlonts  with  0  was  assumed 
to  be  linear  for  0  >  30  deg. 

The  resulting  aerodynamic  coefficients  are  shown  in  Figure  70. 
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Figure  70,  (Concluded) 


It  is  recognized  that  the  assumptions  made  above  could  result  in 
considerable  error  for  extreme  angles  of  attack  and  sideslip;  however, 
the  use  of  full-ranged  data  allowed  the  simulation  to  be  operated 
throughout  a  violent  departure  without  demanding  an  immediate  return  to 
initial  conditions  when  lower  aerodynamic  data  limits  were  reached# 

The  aerodynamic  data  were  stored  in  the  digital  computer  as  look-up 
tables#  Total  forces  and  moments  were  computed  digitally  and  fed  to  the 
analog  systems  of  the  simulator  to  drive  the  cockpit  displays. 

Physical  dimensions  of  the  airplane  were  based  upon  the  P-4J.  Iner¬ 
tias  and  weights  are  summarized  in  Section  II# 

F.  DATA  RECORDED 

Parametric  data  were  recorded  via  three  eight-channel  Brush  record¬ 
ers.  All  Euler  angles,  body  axis  rates,  and  control  deflections  were 
recorded,  as  well  as  accelerations,  altitude,  velocity,  and  thrust  set¬ 
tings. 

For  a  limited  number  of  runs  the  total  body  axis  accelerations 
(p»  q.  ?)  were  recorded  and  inertial  and  aerodynamic  components  of  4 
were  also  recorded  separately. 

G.  TASKS  AMD  MANEUVERS 

Piloting  tasks  were  selected  which  exercised  the  open-  and  closed 
loop  departure  parameters  of  Interest.  These  were  dlvidd  into  two 
phases.  The  first  consisted  of  familiarization  maneuvers  typical  of 
"feeling  out"  stall/departure  motions,  warning,  etc.,  of  a  new  airframe. 
These  are  modeled  after  the  suggested  test  matrix  of  Reference  34  plus 
an  additional  aggravated  input  maneuver  for  determining  departure  sus¬ 
ceptibility,  proposed  in  Reference  33.  The  second  phase  consisted  of 
tracking  tasks  and  maneuvers  as  suggested  in  Reference  36. 
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1.  Handling  Familiarization  and  Aaaemeut 


Five  basic  familiarization  (F)  maneuvers  were  employed: 

a  Straight-ahead  stall  (Fj):  1  g  ,  fixed  thrust , 

longitudinal  stick  pull  to  produce  d  *  1  deg/sec 
until  stall  departure. 

•  Constant  attitude  stall  (F?):  1  g,  wings  level, 

holding  constant  pitch  attitude  while  slowly  re¬ 
ducing  thrust  until  stall/departure. 

•  dank-to-bauk  turns  (F3):  Constant  altitude, 

±60  deg  bank-to-bank  turns  at  increasing  AOA  until 
departure;  rudder  coordination  optional. 

•  Wind-up  turn  (F^):  Constant  thrust,  wind-up  turn 
to  stall  departure. 

•  Full-stick-deflection  maneuver  (F,-)i  From  60  deg 
bank,  rapidly  ramp  full  aft  longitudinal  stick 
followed  by  full  lateral  stick  out  of  the  bank. 

After  8  sec  neutralize  controls. 

With  the  exception  of  maneuver  F5  there  was  no  restriction  on  the 
type  or  magnitude  of  controls  used  by  the  pilots.  In  fact,  pilots  were 
advised  to  investigate  various  control  techniques  Just  as  they  would  in 
a  stall  evaluation  of  a  new  aircraft.  Maneuver  F^  (from  Reference  35) 
uas  evaluated  as  an  indicator  of  departure  susceptibility. 

H.  Precision  Tracking  Evaluation 

The  target  aircraft  <a  TV  picture  of  a  model  airplane,  projected 
onto  the  dome  at  the  proper  aspect  and  size)  was  programmed  to  perform 
three  basic  maneuvers  during  a  slow  pull-up.  Target  climb  rate  was 
fixed  to  produce  4*1  deg/eec  to  guide  the  subject  pilots  into  repeat¬ 
able  tracking  runs  of  approximately  30  sec  minimum  duration  before 
departure.  The  tasks  started  with  the  subject  aircraft  trimmed  at 
u  *  12  dig,  essentially  in  trail.  The  pilot  attempted  to  stay  in  the 
target's  6  o'clock  position  as  if  trying  to  overextend  a  firing 
solution.  The  tracking  (T)  maneuvers  were: 
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•  Straight-ahead  pullup  (T^):  Constant-thrust, 

straight-ahead  pullup  keeping  pipper  on  target 
until  aircraft  no  longer  controllable. 

•  Climbing  roll-reversals  (T2):  Track  target  air¬ 
craft  through  series  of  climbing  roll  rever¬ 
sals  ($  *  ±40  deg)  until  aircraft  no  longer  con¬ 
trollable  .  Constant ' thrust . 

•  Wind-up  turn  (T<j)s  Starting  wings-level,  track 
target  into  count ant-thrust  60  deg  bank  wind-up 
turn  until  aircraft'  no  longer  controllable. 

3.  Departure  Recovery 

Recovery  techniques  employed  were  at  the  discretion  of  the  pilots, 
based  upon  their  experiences  in  fighter  aircraft.  Recoveries  were  ini¬ 
tiated  either: 

•  Immediately  upon  positive  indication  of  uncom¬ 
manded  motion  (R}),  or 

•  After  about  a  3  sec  delay  (Rj)* 

Recoveries  were  also  attempted  in  the  SCAS-on  configurations  (Cases  A2 
and  C2)  by  turning  the  SCAS  off  upon  departure  (R3). 

0.  FLYING  QUALITY  RATING  SCALES 

Prevloua  stall/departure  simulations  for  flying  quality  assessment 
(e.g.,  Reference  5)  have  shown  the  familiar  Cooper-Harper  scale  (Fig- 
w-e  71)  to  be  Inadequate  because  departure  ia  centered  on  the  first  step 
i  v  the  decision  tree,  i.e.,  controllability  instead  of  performance.  It 
was  preordained  for  this  simulation  that  control  would  be  lost  (CH  ■  10) 
because  the  configurations  were  specifically  selected  to  provide  differ¬ 
ing  departure  warning,  severity,  and  recovery  in  order  to  identify  gra¬ 
dations  in  these  attributes.  Therefore,  one  requirement  of  this  program 
was  to  develop  and  validate  a  morn  appropriate  stall/departure  flying 
quality  rating  scale. 

The  Reference  37  experimental  program  succeeded  In  developing  rating 
scales  for  a  similar  loss-of-control  and  recovary  situation  induced 
by  flying  into  a  very  strong  wake  vortex.  This  experience  indicated 
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that  each  attribute  must  be  evaluated  Individually  and  then  the  overall 
hazard  assessed.  Figure  72  is  a  tentative  departure  rating  scale  adapt¬ 
ed  from  the  wake  vortex  experience.  Here  the  stall/departure  attributes 
are  divided  into  four  open-ended  scales.  The  first  is  departure  onset 
warning  —  clarity  or  lack  of  warning.  Next,  the  attitude  and  attitude 
rate  excursions  are  evaluated.  Recovery  characteristics  are  separately 
evaluated  (whether  the  recovery  was  rapid  and  achieved  by  simple  control 
application  or  whether  it  degraded  down  to  the  point  of  no  recovery). 
Finally,  an  overall  assessment  of  hazard  1b  rated. 

In  anticipation  that  overall  hazard  assessment  (acceptability)  might 
vary  with  aircraft  primary  mission,  provision  is  made  for  indicating 
scalar  differences  between  training  and  operational  squadron  usage.  For 
example,  an  operational  aircraft  might  have  very  clear  warning  of  depar¬ 
ture  onset  and  extreme  excursions  if  the  pilot  persisted,  and  be  diffi¬ 
cult  to  recover,  but  be  given  a  relatively  good  hazard  rating,  i.e.,  a 
2,  on  the  basis  of  the  very  clear  and  consistent  warning.  Similarly,  a 
training  aircraft  might  have  little  warning  of  departure  onset  and  rela¬ 
tively  large  attitude  excursions,  but  recover  very  easily  and  quickly  by 
merely  releasing  the  controls,  and  again  might  have  a  relatively  good 
overall  hazard  rating,  such  as  a  2,  on  the  basis  of  recovery. 

As  an  adjunct  to  the  rating  scales,  the  pilot  is  also  asked  to  pro¬ 
vide  qualitative  information  as  to  the  nature  of  the  aircraft  motions 
observed. 

Although  the  rating  scale  of  Figure  72  was  conceived  to  cover  both 
open-  and  cloaed-loup  situations  and  specifically  is  not  spin  oriented, 
it  was  employed  in  a  separate  departure/spin  susceptibility  piloted 
simulation  (Reference  38)  to  obtain  an  early  assessment  of  its  useful¬ 
ness.  Results  were  somewhat  inconclusive  since  the  Reference  38  task 
involved  preprogrammed,  open-loop  pilot  inputs  to  assess  the  departure/ 
spin  susceptibility  of  an  airframe  having  variations  in  surface  author¬ 
ity  and  actuator  rate  limiting  characteristics.  Nevertheless,  valuable 
comments  were  obtained  (Reference  39)  along  with  a  recommended  expansion 
and  restructuring  of  the  pilots'  qualitative  assessment  to  provide 
greater  flexibility  and  provoke  additional  commentary  to  identify  what 
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DEPARTURE  ATTITUDE  AND  ATTITUDE  RECOVERY  OVERALL 

WARNING  RATE  EXCURSIONS  FROM  DEPARTURE  ASSES MENT  OF 

DEPARTURE 

CHARACTERISTICS 


the  pilot  is  rating.  Accordingly,  a  loas-of-control/departure/recovery 
debriefing  guide  was  prepared,  Table  7,  to  augment  the  rating  chart. 
This  guide  addresses  the  four  topics  of  the  Figure  72  rating  scale  plus 
an  additional  section  devoted  to  post  stall  gyrations  (PSG)  characteris¬ 
tics  and  severity.  Provision  is  unde  under  each  topic  for  the  numerical 
rating  (STI)  from  the  Figure  72  scale.  The  pilot  is  requested  to  answer 
those  questions  appropriate  for  the  maneuver  flown.  Most  questions  per¬ 
tain  to  characteristics  of  the  vehicle  motion  and  are  constructed  to 
encourage  simple  yes/no,  multiple  choice,  or  short  written  answers. 
Cooper-Harper  (CH)  ratings  are  also  requested  for  the  two  tasks  which 
Involve  pilot  effort  to  exert  control  on  the  vehicle  (l.e.,  prevent 
departure  and  accomplish  recovery).  Finally,  it  should  be  noted  that 
the  questionnaire  was  prepared  for  use  in  flight  as  well  as  simulation. 
Thus,  some  questions  are  not  appropriate  for  fixed-base  simulation 
(e.g.,  assessment  of  vehicle  acceleration). 
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TABLE  7.  LOSS  OF  CONTROL/DEPARTURE/RECOVERY 
DEBRIEFING  GUIDE1. 


Iph  ni1  (Vnt  r.P1  /n-mirf  IV'iriOf  t'V  1'nh'n  I 

ianawur  ttiuau  qua&tiona  uuc  aro  uppropr mu  tot  thu  maneuver 

flown)  Clrcla  intWara  vharavar  poaalbl*. 


a)  la  warning  clmir/unelanr/nonaKiatant?  - 

b)  Khat  U  (ha  nalur*  of  tlia  warning?  (Aircr.Ut  motion*, 

vibrations,  Inatrinu.'it  i  mil  to  t  loni ,  vigual  cutb,  motion  cuu»,  control 
ijiUffl  (Aal  auua)  . 

c)  la  warning  maxkad  Hy  bom*  othar  aircraft  ciiatactoriat ic ?_ 

<i)  la  cha  winning  auch  that  it  Allow*  tha  aircraft  to  bw 

flown  eloior  to  the  Unit  of  th*  anvalopa?  __ _ — 

a)  Data  cho  warning  Intiifara  with  tit*  prut.iiy  us.c  or 

raduca  Altaian  of factivanoai?  . _ _ 

0  Oo  tha  aircraft  motion*/ vibration*  ctmatiuu  a  mimiih 

haiard  worthy  of  a  flight  raatilction?  _ _ _ 

g)  ia  tlia  turning  to  inadaguau  nut  .i 1 light  ro»t notion 
would  ba  nacuaaary  to  pravant  loaa  of  control!-  _  m _ 

II.  Pnnarturo.  (Sou  NIL-S-8 J691  definition  anil  d  tat  union) 

4)  Dlu  dopurtur*  occur?  _ 

bj  Wliat  won  th«  aircraft  .tiotiotu'/  (Wing  rock,  nunc  yluo, 
pitch  up,  rolling  daparsura,  divergent  oacllUtlnna) 

4}  What  war  tha  itvarlty  of  tha  aircraft  motion?  . 

Mild  «  Mild  Acceleration  tu** 

Incurnadiat*  -  Rapid  motion  dr  acceleration  In  one  or 

more  a*ia, 

Savore-  Vary  rapid  notion  or  acceleration  In  I'mi  or  note 

AXlBt 

1)  How  lor;*  war*  tha  chattgoa  in  airuNil  aLUtuJu? 

2)  How  (*at  ware  tha  rataa?  _  .  _ _ 

J)  How  large  war  a  tha  acceleration*?  _ _ _ _ 

d)  Ware  tha  aircraft  motion*  ati>J  attltudua  mLupn.i.iblu L . 
a)  Wara  tha  altera?*  motion*  dlaor lain lug  and/or  del' illt.it 

l)  Ac eo riling  to  MlL-S-fl3&9i  definition*!  huw  would  tlio  rot  »- 

unca/auaeapcibiUty  to  departure  ba  doucrlbud  for  ,EU  1 *  manuuvur? _ 

g)  Should  pilot  action  praminc  or  da  Lay  clo|inrcure?  (t(  mo 
pilot  action  takon  ao  ataco)  .  ,  _ _ _ 

1)  Uliuc  pilot  actiona  wui*  takun? _ _ 

2)  Wliat  uara/would  ba  tin*  demand*  on  tha  pilot  to 

pravant  JapartuM?  (Kafur  to  Coopar-Harpar  Scale)  _  .  _ . _ 

3)  Old  pilot  actiona  aggravate  departure? _ 

1  ".i)*  Jhat  ware  ^he^treraf caecum*?  - - - - — 

b)  lili at  wai  tha  aavarity  of  tlia  a  (rural  t  motion? _ 

Mild  -  Mild  anceiarotlon*  ami  rata* 

tntarmenUta  -  Rapid  motion  or  aucelarotlim  in  ona  or 

mor*  anla* 

Savara  -  Vary  rapid  motion  or  need  I  aralloti  in  .mu  or 

mor*  oalit 

1)  War*  tha  change*  in  aircraft  attitude  inrm?  _ 

2)  War*  tlia  rote*  feat?  _ _ _ 

1)  War*  tha  acuelaratlona  large? _ 

e)  Wuru  tlia  aircraft  noiion*  dUorltfrtt  mg  and/or  dibiii-  „ 

taking?  _  ^ 

IV  R  eco  wry 

•li  How  rapid  uu*  tha  recovery? 

Immediate 

Slow  •<  Altar  a  pariod  of  lima,  short  oiiougn  cu  iiruvcnt 
doubt  concerning  avantual  recovery 

fiucaiaivo  -  Aftar  a  puriod  of  time,  long  an.wi,;li  to  p  to¬ 
il  uc  a  aurioua  doubta  cuocartting  avancu.il  tocovirry 

b)  Wli.it  raciovary  contiala  wuro  used?  —  - 

c)  Warn  thu  racovary  contrnj  .nplic.H lunui 

Simpla  *  Ona,  two,  or  thrao  actiona  tli.it  do  nut  ruium* 
pilot  pract.'ca  to  h«  affactlva, 

NiU.brui  -  Cururoi  application*  that  would  hu  cont idered 
nutnal  pilot  procudura* 

Complli-.itad  -  Mora  th-in  thru*  .tutluna  und/ur  actinna 
ro«|uirn  uonaidcraiil.,  pilot  praeLtea  tu  bu  aituctivo, 

Uim.il  oral  •  Control  dpi«l  leal  lona  tliat  irr  nn|r|uu  tn  Mil* 
out-iif-eoni ml  roe o vary  and/or  w.iuid  nit  ba  Lunaidotvd  * r ' .i I  pilut 
prui’odura. 

Ag”  r.i'Jiit  1  .in  -  to  eh -i  iil»»>f  mi  ennd  1 1  mn, 
dJ  What  war*  the  diinai.d*  cm  tha  pilvt  to  irromiinh 
rocovary?  (Rofar  tu  Cuopor-Karpor  Snla  whuri  .tppror r i.n n  fnr  _ 

Wi.rc*  das  rlbino  Canamlt  mi  pllni,  Jloili'v  '.-..rda  a*  «'r  i.irv 
cumidurlnu  that  thii  uik  bucan  with  an  not -m-rancr >1  air.  r.ill,) 

V  lUiarl  F.irh‘-«r  _  /r^.riar 


SECTION  VI 


HIGH  AOA  FLYING  QUALITIES  ASSESSMENTS 


Effects  of  the  maneuver-limiting  factors  on  high-angle-of-attack 
flying  qualities  were  assessed  by  three  pilots  utilising  the  prescribed 
set  of  non-tracking  and  tracking  maneuvers.  The  two  principal  pilots 
were  highly  qualified  U.S.  Air  Force  fighter  test  pilots,  graduates  of 
the  Air  Force  Test  Pilot  School  at  Edwards  AFB,  with  a  broad  range  of 
practical  experience.  The  first  pilot,  RC,  has  flown  a  varied  selection 
of  aircraft,  from  the  F-4  to  the  B-52;  the  second  pilot,  JF,  was  experi¬ 
enced  in  operational  fighter  and  attack  aircraft,  including  the  A-9 
stall/spin  evaluation,  and  is  a  former  instructor  at  the  Test  Pilot 
School  responsible  for  stall/spin  curricula.  The  third  pilot,  RH,  was  a 
handling  quality  engineer/general  aviation  pilot  who  also  accomplished 
all  of  the  Initial  simulation  debugging,  procedural  checkout,  primary 
pilot  brieflngB  and  debriefings,  etc. 

The  combination  of  airframe  configurations,  familiarization  and 
tracking  maneuvers,  and  recoveries  resulted  in  the  98-cell  matrix  of 
Table  8.  This  matrix  was  flown  first  as  a  sequence  of  runs  progressing 
horizontally  across  the  various  maneuvers  for  a  given  airframe.  After 
the  complete  matrix  had  been  accomplished  once,  a  second  series  was  per¬ 
formed  progressing  vertically  through  the  airframe  configurations  for 
each  tracking  task. 

A  total  of  1088  formal  runs  were  conducted.  These  were  fairly 
evenly  divided  between  the  various  configurations  and  pilots,  as  indi¬ 
cated  by  Table  9. 

One  of  the  first  taskB  of  the  simulation  was  to  validate  or  modify 
the  ratings  scale  and  debriefing  questionnaire  to  the  satisfaction  of 
the  evaluation  pilots.  The  final  versions  were  used  throughout  the 
remainder  of  the  simulation.  Accordingly,  this  section  starts  with  the 
rating  acale  assessment  and  then  presents  the  airframe  configuration 
assessments. 
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TASK/KECOVERY  TEQJNIQtJE 


TABLE  9.  NUMBER  OF  RUNS 


CONFIGU¬ 

RATION 

PILOT 

TOTAL 

RC 

JF 

RH 

A1 

59 

79 

79 

217 

a2 

68 

48 

79 

195 

B 

61 

94 

58 

213 

C1 

44 

77 

40 

161 

c2 

53 

62 

58 

173 

D 

35 

62 

32 

129 

TOTAL 

320 

422 

346 

_ 

1088 

A.  RATINGS  SCALE  AMD  QUESTIONNAIRE  ASSESSMENT 

The  debriefing  questionnaire  was  answered  at  the  completion  of  the 
familiarization  runs  and  again  at  the  completion  of  the  first  set  of 
tracking  runs.  During  this  period  modifications  were  worked  out  for  the 
rating  scale.  The  modified  scale  was  then  used  for  the  final  Bet  of 
tracking  runs. 

1 .  Questionnaire 

The  debriefing  questionnaire  proved  essential  to  obtaining  a  consis¬ 
tent  set  of  information  for  interpretation  of  results.  However,  the 
simple  yes  or  no  answers  requested  did  not  suffice.  The  tasks  and 
resulting  aircraft  motions  were  so  complex  as  to  require  considerable 
explanation.  For  example,  more  than  a  single  mode  of  departure  or  spin 
was  observed  in  some  cases.  Departure  warning  and  severity  also  varied 
with  rate  of  stall  onset  (&)  and  pilot  usage  of  controls.  Thus  accom¬ 
plishment  of  the  questionnaire  often  became  quite  tlme-coneuming,  with 
several  pages  of  accompanying  explanation.  This  was  augmented  by  tape 
recordings  made,  on-line,  of  pilot  commentary  while  performing  the  simu¬ 
lation  runs.  As  a  consequence,  pilot  responses  were  so  lengthy  —  and 
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in  so  many  cases  repetitive  -  that  it  is  not  practical  to  incorporate 

them  in  this  report. 

However,  the  basic  conclusion  reached  by  both  pilots  and  experi- 
mentors  was  that  such  a  questionnaire  is  highly  desirable  for  stall/ 
departure/recovery  investigation.  The  only  change  in  format  deemed 
necessary  is  provision  of  more  space  to  accommodate  lengthy  answers. 

2.  Rating  Scales 

As  in  the  previous  simulation  (Reference  5),  the  pilots  felt  that 
the  Cooper-Harper  rating  scale  was  inappropriate  for  stail/depat ture/ 
recovery  evaluation  where  aircraft  characteristics  rather  than  task 
performance  are  being  assessed.  On  the  other  hand,  the  overall  concept 
of  the  Figure  72  departure  ratings  scale  was  appreciated,  although  somu 
alterations  were  required  before  this  scale  met  with  full  approval.  In 
particular,  the  pilots  did  not  like  the  open-ended  aspects  of  the  scales 
(the  scale  extends  beyond  the  first  and  last  descriptors).  This  clearly 
gave  problems,  and  they  preferred  that  the  scale  be  bounded  with  a 
simple  descriptor  at  each  end.  Second,  the  time  required  to  recover  and 
the  control  application  necessary  to  effect  recovery  needed  to  be  sepa¬ 
rated.  Recovery  controls  also  deserved  a  further  breakdown,  with  one 
rating  for  required  control  complexity  and  a  second  rating  for  control 
application  timing.  Finally,  it  was  decided  that  an  assessment  of 
hazard  was  difficult  to  make  because  this  changes  with  situations  such 
as  low  versus  high  altitude.  The  pilots  preferred  to  provide  an  overall 
assessment  baaed  upon  the  acceptability  of  the  flying  characteristics  at 
high  angles  of  attack.  The  end  point  descriptors  selected  by  the  pilots 
were,  in  their  jargon,  "Sierra  Hotel,"  meaning  they  would  like  to  fly 
the  aircraft  at  high  AOA  and  "Delta  Sierra,"  meaning  unacceptable  flying 
characteristics.  The  resulting  rating  scale  is  shown  in  Figure  73. 
Both  test  pilots  were  happy  with  this  final  form  and  considered  that  it 
covered  all  key  factors,  minimized  the  descriptor  conflicts,  and  did  not 
constrain  the  pilot's  quantitative  evaluation. 

The  rating  scale  of  Figure  73  was  employed  in  the  final  set  of 
tracking  tasks.  Resulting  characteristic  assessments  for  the  six 
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Figure  73.  Modified  Departure  Rating  Scale  (Mod  II) 


vehicle  coni i gurations  are  summarized  In  Figure  74.  The  numerical 
ratings  shown  for  each  pilot  are  the  averages  of  Beparate  ratings  given 
for  each  of  the  three  tracking  tasks  (Tj_  through  T3).  Considering  the 
small  statistical  run  sampling  and  the  fact  that  all  configurations  were 
intended  to  have  basic  handling  deficiencies,  the  results  are  considered 
to  be  very  good.  It  is  to  be  expected  that  all  ratings  would  tend 
toward  the  poor  (high  numerical  value)  end  of  the  scales  and  they  gener¬ 
ally  are.  Most  importantly,  they  also  reflect  the  specific  dynamic  and 
handling  variations  intended  to  be  exhibited  between  configurations  and, 
in  most  instances,  the  ratings  of  RC  and  JF  demonstrate  agreement  in 
variational  differences  if  not  in  absolute  levels.  The  ratings  of  pilot 
RH  tend  to  vacillate  between  those  o t  RC  and  JF  and  also  (except  for 
recovery  control  timing)  tend  to  agi.ua  with  the  poorer  rating.  This 
downrating  probably  reflects  pilot  RH'g  relative  inexperience  in  depar¬ 
ture/  opin/recovery  and  the  violent  PSGs  which  high-performance  fighters 
often  exhibit.  Thus  the  ratings  of  RC  and  JF  are  considered  the  more 
significant. 

The  separation  between  RC  and  JF  ratings  bn  warning  and  motion 
severity  is  consistent  with  differences  in  piloting  techniques.  Gener¬ 
ally  JF  is  about  one  rating  point  harsher  then  KC  because  of  a  more 
rapid  pull  into  departure  and  spin.  This  will  be  discussed  in  more 
detail  later. 

Recovery  control  timing  assessment  shows  very  good  agreement  between 
the  pilots  except  for  Configuration  C2.  This  primarily  reflects  timing 
of  control  release.  The  differential  is  about  one  rating  point  or  less 
except  Configuration  C£  where  RC  had  difficulty  with  the  two  different 
spin  modes  the  aircraft  could  have.  Timing  of  control  release  deter¬ 
mined  whether  an  oscilletory  (recoverable)  or  flat  (unrecoverable)  spin 
was  obtained. 

Time  to  recover  and  overall  rating  show  almost  Identical  trends  and 
thus  reflect  predominant  concern  for  loss  of  altitude  in  any  departure/ 
spin  situation.  Configuration  C^  had  an  oscillatory,  recoverable 
spin.  Configuration  D  did  not  spin  but  could  end  up  with  appreciable 
altitude  loss  if  the  pilot  kept  fighting  the  "wallowing"  departure  mode 
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(described  in  Subsection  II. B).  All  other  configurations  had  flat  spin 
modes  that  were  difficult  if  not  impossible  to  stop.  The  significant 
disapproval  of  the  Configuration  D  pitch-up  is  also  evident  for  pilot  RC 
on  both  of  these  ratings. 

The  rating  scale  for  recovery  control  complexity  may  need  additional 
refinement.  Pilot  JF  rated  all  configurations  at  1  because  he  merely 
released  controls  to  initiate  recovery  and.  by  definition,  this  is  as 
simple  as  things  can  be.  Pilot  RC,  with  less  severe  departures,  was 
able  to  discern  some  influence  of  the  Btick  and  rudder  in  effecting 
recovery  and  this  was  also  reflected  in  his  being  more  critical  of 
recovery  control  timing. 

However,  aspects  pertaining  to  recovery  (control  usage  and  recovery 
time)  for  configurations  ocher  than  D  are  not  necessarily  significant 
because  the  simulation  cannot  be  considered  valid  for  spin  and  spin 
recovery.  The  motions  were  generally  so  violent  and  the  out-the- 

windscreen  sky/earth  display  so  featureless  that  the  pilots  had  to 
resort  to  the  all-attitude  ball  for  recovery.  This  is  not  completely 
unnatural  for  the  military  pilots,  but  did  make  the  task  more  diffi¬ 
cult.  Additionally,  the  aerodynamic  coefficients  employed  In  the 
simulation  do  not  represent  steady  spin  conditions.  Probably  the  most 
that  can  be  attributed  to  recovery  evaluations  is  that  they  exercised 
the  rating  scales. 

In  summary,  the  departure  rating  scales  of  Figure  73  were  accepted 
and  supported  by  the  pilots  and,  based  upon  this  small  statistical 

sample,  the  numerical  ratings  given  reflect  characteristics  "designed 
into"  the  vehicle  configurations.  The  spread  in  ratings  between  the  Air 
Force  teat  pilots  (generally  about  O.S  on  a  5-point  scale)  is  consistent 
with  spreads  normally  expected  in  Cooper-Harper  ratings  (about  1.0  on  a 
10-polnt  scale).  Greater  spreads  generally  reflect  sensitivity  to 
piloting  technique;  that  sensitivity  itself  may  be  an  attribute  of  the 
scales. 

Since  vehicle  attributes  rather  than  a  task  performance  level 

are  being  rated,  use  of  the  rating  scales  must  be  accompanied  by  a 
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qualitative  assessment/description.  The  qualitative  information  should 
include,  as  a  minimum,  the  following: 

•  Warning 

-  Type 

-  Clarity 

-  Margin 

•  Departure 

«  Resistance  (susceptibility) 

-  Type 

-  Severity 

Ability  of  pilot  to  delay  or  prevent 
•  Control  action  taken 

e  Demands  on  the  pilot 

•  Post-Departure  Motion 

-  Type  of  aircraft  motion 

-  Severity 

•  Recovery 

-  Rapidity 

-  Recovery  controls 

-  Demands 

«  Ability  to  recognise 

a  Ability  to  perform  necessary  control  action 

B.  COHFIGORATIOH  ASSESSMENTS 

As  noted  previously,  the  goals  were  to  validate  the  influence  of  the 
key  aerodynamic  coefficients  in  determining  departure  characteristics; 
evaluate  the  influence  of  varied  maneuver-limiting  factors  on  high  AOA 
maneuvering  control;  and  identify  potential  flying  quality  criteria 
in  terms  of  departure  resistance,  warning,  and  severity  for  possible 
application  to  the  flying  quality  specification.  In  this  subsection, 
predicted  vs.  observed  characteristics  aie  reviewed  and  demonstrate 
that  the  above  goals  have  been  achieved,  although  therw  were  seme  sur¬ 
prises  in  the  process.  The  unaugmented  airframe  configurations  will  be 
assessed  firot,  and  then  the  augmented  cases. 
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It  is  necessary  to  rely  heavily  on  qualitative  statements  contained 
in  the  debriefing  questionnaire  and  commentary  taped  on-line  as  the 
pilots  talked  themselves  through  the  runs.  The  questionnaires  were 
filled  in  only  during  the  initial  progression  through  the  test  matrix. 
At  that  time  the  pilots  were  "feeling  out”  the  various  configurations 
and  produced  voluminous  commentary  on  each,  covering  the  primary  obser¬ 
vations  plus  variations  due  to  different  control  applications,  tech¬ 
niques,  rates  of  onset,  tracking  maneuvers,  etc.  The  assessments 
reported  in  the  following  paragraphs  have  been  distilled  to  the  simplest 
possible  factors  consistent  with  the  goals  of  the  simulation. 

Only  assessments  of  the  fighter  test  pilots  are  reported  because  RC 
and  JF  are  most  experienced  in  stall/departure/spln  testing  and  were 
more  observant  of  differences  between  configurations.  It  became 
apparent  early  in  the  data  analysis  that  the  two  test  pilots  were 
employing  almost  opposite  approaches  to  stall/departure.  RC  employed  a 
cautious,  slow  increase  in  AOA  and  generally  was  able  to  detect  subtle 
changes  in  vehicle  stability  or  response  characteristics.  JF  used  an 
aggressive,  rapid  AOA  Increase  which  did  not  give  time  to  detect  such 
warnings.  As  a  consequence,  he  generally  pulled  to  a  higher  AOA  with 
leas  speed  bleedoff  and  obtained  more  severe  post-stall  gyrations 
(PSG).  Thus,  differences  in  pilot  closed-loop  control  techniques  pro¬ 
duced  a  confounding  influence  which  must  be  recognized  before  delving 
into  other  details)  departure  warning  and  severity  are  a  function  of 
pull-up  rate  and  lateral-directional  controls  application  immediately 
preceding  departure. 

1 .  Unangawnted  Airframe 

a.  Configuration  A} 


Predicted  versus  observed  pre-  and  post-departure  characteristics 
are  shown  in  Table  10.  The  warning  Pilot  RC  observed  with  increasing 
AOA  is  consistent  with  the  graceful  (gradual)  degradation  suggested 
by  the  open-loop  dutch  roll  and  roll  numerator  root  migrations  of 
Figure  55.  It  appears  the  negative  Cm^  contribution  was  observed  as  a 
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TABLE  10.  CONFIGURATION  A^  MANEUVER-LIMITING  FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning: 

Wing  rock 

(16  <  a  <  22  deg) 

Wing  rock  low  ampli¬ 
tude,  high  frequency 

None 

Roll  revereal 

Cot  >  IB  deg) 

Roll  reversal 

(a  >  20  deg)i  g-break 

(20-22  deg  a) 

None;  g-break 

Wing  rock  large  ampli¬ 
tude,  low  frequency 

Rating: 

2.8 

4.8 

Departure: 

Noee  slice 
(a  >  22  deg) 

. 

Nose  slice  followed  by 
roll 

Nose  slice  and  roll 

Poet  stall  gyrations 

Large  yaw,  pitch,  roll 
oscillation  and  spin 

Violent  ooclllatlons 
and  spin 

Severity  Rating: 

2.2 

4.6 

Departure /Spin 
Susceptibility 
(Weissman  criterion): 

o  <  21  deg:  mild  or 
no  departure 

• 

Release  o  -  20  deg; 
yaw  initially  slow, 
builds  rapidly 

21  <  a  <  24  deg: 
moderate  rolling 
departure 

Release  a  ■  23  deg; 
abrupt  nose  slice, 
builds  exponentially 

a  >  24  deg:  strong 
rolling  departure, 
high  spin  suscepti¬ 
bility 

Rating: 

Resistant 

Extremely  Susceptible 

Comments : 

Maneuvers  tend  to  mask 
warning;  not  masked  if 
pay  attention 

Most  violent  post 
stall  gyration  oscil¬ 
lations  of  any  con¬ 
figuration 
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"g-break."  The  nose  slice  departure  is  as  predicted,  but  the  onset  is 

noticed  prior  to  Cnft  becoming  negative.  Thus  onset  might  be  influ- 
pdyn 

enced  by  pilot  control,  and  thus  associated  with  negative  l/T^.  Pilot 
RC  considered  the  configuration  departure  resistant  (R)  because  of  ample 
warning  and  his  tendency  to  initiate  recovery  at  a  *  20  deg  where  the 
Welssman  criterion  predicts  mild  or  no  departure.  Pilot  RC  therefore 
gave  fairly  good  high  AOA  ratings. 

Pilot  JF  eaw  no  warning  except  the  "g-break."  He  obtained  an  abrupt 
nose  slice  and  violent  poet-stall  gyration  (PSG).  The  tendency  to  flat 
spin  may  be  due  to  the  large  negative  Cn^  and  small  C$B  at  hia  recovery 
initiation  point  and  the  further  destabilizing  influence  of  sideslip  on 
the  denominator  and  numerator  roots,  all  of  which  portend  high  body-axis 
yaw  rate  and  little  roll.  Inertia  effects  would  tend  to  raise  the  nose 

and  align  the  yuw  rotational  axis  with  the  velocity  vector - a  flat 

spin. 

This  configuration  also  exhibited  violent  oscillatory  spin  charac¬ 
teristics  which  could  be  reduced  if  a  large  nose-down  pitch  rate  devel¬ 
oped  as  the  nose  was  slicing.  This  could  be  produced  via  forward  stick 
coinciding  with  the  "g-break."  The  resulting  rq  inertia  crosa-ccmpling 
produced  a  large  roll  acceleration  and  quite  wild  gyrations  (s.g., 
Figure  75).  Pitot  JF  noted  this  configuration  to  have  the  most  violent 
PSG  characteristlca  of  any  configuration  and  rated  it  extremely  suscep¬ 
tible  (ES)  to  departure. 

While  the  Weiasman  criterion  correctly  predicted  the  spin  suscep¬ 
tibility,  the  initial  departure  was  in  yaw  rather  than  roll.  The 
maneuver-limiting  factors  are  the  strong  nose  slice  and  spin  character¬ 
istics.  The  causal  factors  are  combined  open-  and  closed-loop  (roll) 
instabilities  accompanied  by  signifeant  negative  C^. 

b.  Configuration  B 

As  predicted,  this  configuration  is  similar  to  Configuration  A,  but 
with  less  departure  warning  (see  Table  li).  The  low-frequency  yaw/roll 
oscillation  that  Pilot  RC  observed  above  21  deg  AOA  is  consistent  with 
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TABLE  11.  CONFIGURATION  B  MANEUVER-LIMITING  FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning: 

Roll  ravaraal 
(a  >  18  dag) 

Resists  roll  command; 
g-break 

Low  frequency  yaw/roll 
21  deg  a 

No  warning 

Rating: 

4.2 

4.8 

Departure: 

Noaa  alien 
(a  >  22  dag) 

Nose  slice 

Nose  slice 

Poat  atall  gyrationa 

Violent  roll  reversal 
or  spin 

Severe  and  violent 
oscillation 

Severity  Rating: 

2.7 

4.5 

Departure/Spin 

Susceptibility 
(Welsaman  Criterion) 

a  <  21  deg:  mild  or 
no  departure 

Release  a  *  21  dag; 
yaw  followed  by  roll 
and  pitch  down 

21  <  a  <  24  deg: 
moderate  rolling 
departure 

Releaae  o  *  24  deg 

a  >  24  deg:  strong 
noaa  slice,  high  apin 
susceptibility 

Nose  drift  followed  by 
roll;  prone  to  flat 
spin 

Rating: 

Susceptible 

Extremely  Susceptible 

Comments: 

Cannot  discern  onset 
of  departure;  insidi¬ 
ous 

Low  acceleration  in 
noao  slice 
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open-loop  dutch  roll  dynamics  (see  Figure  57).  The  lower  frequency 
resulted  in  less  acceleration  Into  the  nose  slice.  Coupled  with  roll 
reversal ,  this  apparently  masks  the  departure  onset.  Although  he  Ini¬ 
tiated  recovery  at  about  the  same  AOA  and  thought  the  departure  severity 
(pilot  rating  ■  2.7)  Is  about  the  same  as  Configuration  Aj,  he  consid¬ 
ered  Configuration  B  to  be  departure  susceptible  (S)  because  of  the  poor 
warning  (PR  ■  4.2).  Pilot  JF  could  see  little  difference  between  A^  and 
B  except  the  increased  tendency  to  flat  spin  due  to  the  combined  high 
aerodynamic  roll  rate  damping  and  strong  directional  divergence.  He 
again  rated  it  extremely  susceptible  (ES). 

This  configuration  exhibited  all  anticipated  flying  characteris¬ 
tics.  The  Welssman  criterion  correctly  identified  departure  end  spin 
severity  but  not  type.  Again,  the  maneuver-limiting  factors  are  the 
severe  nose  slice  and  flat  spin  tendencies*  The  causal  factors  are  the 
open-  and  closed-loop  directional  instabilities  aggravated  by  high 
negative  C^* 

c.  Configuration 


Configuration  Cj  (Table  12)  with  positive  dynamic  stability  (due  to 
large  negative  C^)  had  a  much  greater  tendency  to  develop  roll  oscilla¬ 
tions.  Both  pilots  noted  that  any  lateral  control  input  excited  wing 
rock  and  produced  a  tendency  to  lateral  PIO.  Since  this  sensitivity 
persisted  over  a  wide  AOA  range,  it  was  not  considered  a  warning  but 
actually  a  masking  of  departure  onset.  Both  pilots  considered  the 
warning  to  be  poor  (PR  >  4). 

This  configuration  is  stable  open-loop  end  driven  unstable  and  into 
departure  by  closing  the  attitude  loops.  Departure  characteristics  are 
less  consistent  (predictable)  because  the  nature  of  departure  depends 
upon  pilot  control  activity*  It  exhibited  at  least  two  and  possibly 
thrse  different  types  of  departure  onset,  which  further  aggravated  the 
poor  warning.  The  lateral  phugoid  mode  could  contribute  the  additional 
departure  mode  (see  Figure  59).  One  pilot  detected  more  tendency  to  yaw 
and  the  other  more  tendency  to  roll.  However,  they  generally  were  in 
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TABLE  12.  CONFIGURATION  C1  MANEUVER-LIMITING  FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning: 

Wing  rock 

(16  <  a  <  26  dag) 

Wing  rock:  any 
disturbance 

Wing  rock:  tendency 
to  low  frequency 
pilot-induced 
uec.lTlatlan 

Roll  reversal 
(a  >  18  deg) 

Low  frequency,  high 
amplitude  oscillation 
masks  departure 

Rating: 

4.3 

4.0 

Departure: 

Roll  departure 

Post-stall  gyrations 

Nose  slice 

Roll 

Two  distinct  types: 

a)  roll,  little  yawi 

b)  yaw,  followed  by 
roll 

Severity  Rating: 

1.6 

3.7 

Departure/Spin 

Susceptibility: 

(Weissman  Criterion) 

a  <  21  deg:  no  depar¬ 
ture 

a  >  21  deg:  mild 
rolling  departure,  low 
spin  susceptibility 

Release  a  *  32.5  deg, 
nose  slice  and  roll 

Release  a  *  25.5  deg, 
milder  departure, 
oscillatory  spin 

Rating: 

Susceptible 

Susceptible 

Comments: 

Takes  tine  to  see  yaw 
ie  not  going  to  stop 

Not  prone  to  flat  spin 
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closer  agreement  for  this  configuration  concerning  warning,  nature,  and 
severity  of  and  susceptibility  to  departure,  probably  due  to  the  posi¬ 
tive  dynamic  stsbility  which  prolonged  the  time  period  the  airplane 
could  be  maintained  under  semi-control.  Thus,  both  pilots  had  suffi¬ 
cient  time  to  observe  its  high  AOA  characteristics  fully  before  closed- 
loop  instability  predominated. 

Both  pilots  considered  the  departure  and  FSG  severity  to  be  somewhat 
lees  than  for  either  A^  or  B.  The  configuration  was  spin-prone,  but 
spins  were  generally  highly  ocillatory  with  some  possibility  of  recov¬ 
ery.  It  was  not  prone  to  flat  spins.  Departure  and  spin  characteris¬ 
tics  and  susceptibility  are  in  agreement  with  the  Weissman  criterion. 

The  initial  or  basic  maneuver-limiting  factor  1b  the  strong  wing 
rock  tendency  which  extends  from  relatively  low  AOA  through  departure 
and  is  aggravated  by  closed-loop  instability  (roll  PIO).  Departure  is  a 
secondary  factor.  Warning  may  also  be  adversely  Influenced  by  the  low- 
frequency,  large-amplitude,  lateral  phugoid  oscillation. 


d.  Configuration  D 


This  configuration  (Table  13)  had  less  directional  Instability  but 
did  have  sideslip-induced  pitch-up.  The  expected  wing-rock  tendency 
from  a  lightly  damped  dutch  roll  was  not  observed  by  either  pilot.  Both 
noted  departure  onset  as  a  slow,  persistent  yaw  (noBe  slice)  accompanied 
by  mild  roll  reversal  and  pitch-up.  The  onset  was  sufficiently  mild 
that  departure  was  semi-controllable  throughout  the  achievable  AOA 
range.  A.  poor  departure  warning  rating  (PR  >  3.5)  resulted. 


Open-  and  closed-loop  Instability  parameters  (Wj  and  u)2  or  Cn^ 
and  LCDP)  are  small  for  this  configuration  and  lead  to  sloppy  control 
and  low-frequency  wallowing  which  masks  departure.  At  the  same  time, 
the  wallowing  does  not  generate  sufficiently  rapid  motion  to  excite 
Inertia  croee-coupllng  and  PSG.  All  pilots  tended  to  continue  fighting 
to  maintain  control  well  past  full  stall,  incurring  excessive  eltitude 
loas.  However,  if  controls  were  released  at  any  time  the  aircraft  would 
immediately  go  into  a  noae-low  spiral  and  recover  by  itself. 
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TABLE  13.  CONFIGURATION  D  MANEUVER-LIMITING  FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning] 

Wing  rock 

(16  <  o  <  26  deg) 

Roll  reversal 
(a  >  20  deg) 

Yaw,  roll  reversal, 
pitch  up 

Sloppy  control,  large 

6  and  $  overshoots 
mask  departure 

Large  definite  yaw 
followed  by  reversal 
and  pitch  up 

Rating: 

4.0 

3.5 

Departure: 

Pitch  up 

Nose  slice/pitch  up 

Nosa  slice/pitch  up 

Poet  atall  gyration 

None,  never  spin 

Mild,  no  spin 

Departure /Spin 
Susceptibility 
(Weissman  Criterion): 

All  a:  no  departure 

Release  a  *  27.6  deg, 
slow  nose  slice,  pitch 
up,  some  roll 

Release  a  *  27.6  deg; 
slow  consistent  nose 
slice,  roll,  departure 
nose  high 

Rating: 

Resistant 

Resistant 

Comments: 

Do  not  like  pitch  up 
at  departure;  cannot 
get  nose  down 

Always  recoverable  If 
releaae  controls 
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Departure  warning  was  considered  poor  (e.g.,  4.0)  by  RC  because  of 
the  wallowing,  but  at  the  sane  tine  the  configuration  was  considered 

departure  resistant  (R).  Pilot  JF  assessed  It  as  resistant  primarily 

because  it  would  not  spin.  Pilot  RC's  principal  objection  to  this 
configuration  was  the  positive  C^.  Even  though  it  eliminated  the 

severe  pitch-down  in  nose  slice  and  the  resulting  Inertia-coupled  PSG 

exhibited  by  other  configurations,  he  did  not  like  the  nose  rising,  no 
matter  how  mildly,  when  in  or  near  departure  conditions. 

The  maneuver-limiting  factor  for  this  configuration  is  the  wallow¬ 
ing  -  a  general  insensitivity  to  or  sloppiness  of  response  to  pilot 

control  inputs  prior  to  the  departure  AOA.  The  key  causal  factors  are 
considered  to  be  the  relatively  small  (due  to  small  %  and  C&g) 
which  produces  very  low  DC  gain  in  the  roll  frequency  (Bode)  domain  and 
positive  which  apparently  contributed  to  overcontrol  in  pitch. 

e.  Summary 

The  intended  high  AOA  departure  characteristics  were  obtained  in  the 

four  unaugmented  airframe  configurations.  The  results  demonstrated  the 

influences  that  changes  in  the  six  key  aerodynamic  coefficients  cun  have 

on  departure  warning,  severity,  and  susceptibility.  The  .susceptibility 

of  Configuration  C^  to  departure  was  somewhat  revealing:  the  open-loop 

parameter  Cn.  ,  by  Itself,  is  not  a  sufficient  indicator  of  departure 
sdyn 

resistance.  However,  when  combined  with  the  closed-loop  parameter  LCDP 
(or  w^),  the  two  are  quite  accurate  indicators  of  departure  suscepti¬ 
bility  and  severity  as  advanced  by  Weissman. 

Departure  severity  also  appears  influenced  by  the  static  cross¬ 
coupling  derivatives  <£&,  Nq,  and  Mg.  Severe  departure  and  PSG  were 
obtained  for  the  two  configurations  :f.n  which  combinations  of  these 
derivatives  shifted  lateral-directional  denominator  and  numerator  roots 
into  the  right  half-plane  of  the  root  locus  plot  with  i  4  0  (Figures  55 
and  57).  Lesa  revere  departure  was  obtained  for  combinations  which 
primarily  shifted  numerator  roots  info  the  right  half-plane  (Figure  59), 
and  the  mildest  departure  was  obtained  for  combinations  which  produced 
very  little  shift  of  lateral  denominator  and  numerator  roots  with 
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sideslip  (Figure  61).  The  literal  approximate  factors  developed  in 
Section  III  indicate  that  the  small  shifting  is  due  to  the  combined 
Influence  of  small  N£  and  positive  Mg.  However,  the  nose  rise  asso¬ 
ciated  with  large  positive  Mg  was  disconcerting  to  one  pilot. 

The  simulation  results  also  demonstrated  that  wing  rock  or  lateral 
PIO  is  not  a  satisfactory  warning  of  impending  departure  if  ''.here  is 
sufficient  elevator  power  to  pull  through  the  warning  region  rapidly  or 
if  the  AOA  range  in  which  it  occurs  is  too  wide. 

2.  Augmented  Airframe 


a.  Configuration  A2 


The  flight  control  augmentation  was  configured  to  provide  a  high 
level  of  roll  rate  damping,  and  to  improve  roll  maneuvering  response  via 
roll  rate  command  (pc)  and  reduced  adverse  aileron  yaw  (SRI).  This 
augmentation  was  also  expected  to  eliminate  all  cues  of  impending  depar¬ 
ture.  The  departure  was  anticipated  to  be  a  strong  nose  slice  due  to 
increasingly  negative  Cng  as  AOA  increased  and,  once  sideslip  started  to 
build,  it  was  expected  that  the  characteristics  would  be  very  similar  to 
those  of  Configuration 

Results  summarized  in  Table  14  for  Pilot  JF  are  very  much  as  antici¬ 
pated.  However,  Pilot  RC  detected  a  pre-nose-slice  resistance  to  roll 
commands  (at  AOA  above  approximately  21  deg)  when  attempting  to  track 
the  maneuvering  target.  This  resistance  was  observed  only  during  the 
initial  tracking  sets  and  resulted  in  the  pilot's  neutralizing  all  con¬ 
trols  before  the  directional  instability  became  severe  (rating  1.8). 
The  aecond  set  of  tracking  runs  were  done  mure  hurriedly,  the  warning 
was  not  obtained  (rating  4.3),  end  more  severe  departures  ensued  (rating 
3.7).  The  susceptibility  assessment  as  resistant  (R)  by  Pilot  RC  was 
given  only  on  the  initial  set  of  rune. 

The  maneuver-limiting  factor  for  this  configuration  is  ths  nose 
■lice.  But.  since  there  la  no  warning,  en  AOA  limit  wculd  have  to  be 
set  to  provide  a  safety  margin  against  lnadvsrtsnt  departure.  The  key 
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TABLE  14.  CONFIGURATION  A2  MANEUVER-LIMITING  FACTORS 


PREDICTED 

OBSERVED 

1 

RC 

JF 

Warning : 

None 

Resists  above 

21  deg  a 

Cluar  warning  with 
target 

None 

Rating: 

2.8 

(4.3  on  reruns) 

4.8 

Departure : 

Nose  slice 

Nose  slice 

Nose  Blice 

Roll  and  severe  PSG 

Severity  Rating 

1.8 

(3.7  on  reruns) 

4.5 

Departure /Spin 
Susceptibility 

o  <  23  deg:  no  depar¬ 
ture 

a  >  23  deg:  nose 
slice  and  revert  to 
Configuration  A| 
characteristics 

Release  a  *  25  deg; 
smooth,  slow  yaw 
followed  by  roll 

Release  a  A  32  deg; 
exponential  yaw  diver¬ 
gence  tendency  to  flat 
spin 

Retina: 

Resistant 

Extremely  Susceptible 

Comments : 

Depart  without  pilot 
aware;  only  cue  is 
resistance  to  $c 
against  target 

Easier  to  fly 

Rate  of  yaw  dependent 
upon  4  or  g 
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causal  factors  are  the  strong  static  directional  divergence  and  the 
elimination  of  departure  warning  by  the  augmentation  system. 

Time  traces  of  a  Configuration  A2  departure  from  a  winga-level  pull- 
up  are  shown  in  Figure  76.  The  lateral  Btick  trace  reflects  about  the 
same  activity  prior  to  departure  as  in  Figure  75  but  far  less  sideslip 
and  roll  oscillation.  After  pilot  neutralization  of  controls,  the  aug¬ 
mentation  loops  continue  to  deflect  the  lateral  and  directional  control 
surfaces  and  apparently  aid  recovery  (again  compare  with  Figure  75). 
Interestingly,  the  actual  departure  and  PSG  motions  of  Figure  76  are 
remarkably  similar  to  Reference  38  flight  traces  for  the  F-4E  with 
augmentation  systems  on  (Figure  77).  Prior  to  departure  the  traces  of 
the  two  figures  are  not  comparable  because  the  stick  to  rudder  intercon¬ 
nect  in  the  simulation  reduced  sideslip  excitation  and  provided  better 
roll  attitude  control. 


b.  Configuration  C2 


Based  on  the  positive  Cn^  and  Improved  (or  LCDP),  this  aug¬ 
mented  configuration  was  expected  to  have  the  best  high  AOA  flight 
characteristics  and  little  or  no  departure  tendency.  Xt  did  allow 
both  pilots  to  consistently  maintain  control  to  higher  AOA  than  did  any 
other  configuration,  though  abrupt  nose  slice  departures  were  common 
(Table  15).  In  this  case  the  roll  damper  countered  the  dynamic 
stability  contribution  and  allowed  the  large  negative  Cn^  otatlc  insta¬ 
bility  to  predominate.  The  aircraft  responses  by  themselves  provided  no 
warning  of  Impending  departure,  although  Pilot  RC  was  able  to  extract 
soma  low-frequency  PIO  cues  with  respect  to  the  target  aircraft.  He 
rated  the  configuration  doparture  resistant  (R)  because  of  this  cue  and 
the  mild  (to  him)  nose  olice  motion.  Without  the  target  he  rated  the 
configuration  departure  eueceptible  (S) ,  as  did  Pilot  JF. 


This  configuration  was  also  spin-prone  and,  similar  to  C^,  exhibited 
two  distinct  spin  modes  as  noted  in  Table  15.  It  appeared  to  be  more 
prone  to  spin  flat  than  was  the  unaugmented  configuration  (C^),  because 
the  roll  rate  CAS  automatically  opposed  any  uncoromandad  roll. 


TABLE  15.  CONFIGURATION  C2  MANEUVER-LIMITING  FACTORS 


PfJSDICTED 

OBSERVED 

RC 

JF 

Warning: 

None 

Low  frequency  pilot- 
induced  oscillation; 
g-break 

None,  very  steady 

Rating: 

3.2 

4.8 

Departure: 

None 

Nose  slice  followed  by 
roll  and  pitch  down 

Abrupt  nose  slice 

Severity  Rating: 

1.5 

3.8 

Departure/Spin 

Susceptibility: 

Reslitant 

Release  a  *  31  deg; 
slow  nose  slice; 
almost  always  spin 

Release  o  *  30  deg; 
abrupt  nose  slice;  two 
Bpin  modes: 

a)  oscillatory 
(recoverable) 

b)  flat  (nonrecover- 
able) 

Rating: 

Resistant 

Susceptible 

Comments: 

Resistant  because  tar¬ 
get  gives  cue;  without 
target  is  Susceptible 

Better  flying  quali¬ 
ties  prior  to  depar¬ 
ture 
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The  maneuver-limiting  factor  for  this  configuration  Is  the  nose 

slice  departure.  The  positive  Cna  and  ui^  throughout  the  achievable 

“dyn  t 

AOA  range  did  allow  it  to  be  controlled  to  considerably  higher  ADA  than 
the  other  configurations.  In  this  respect  it  was  less  susceptible  to 
departure.  However,  at  the  higher  AOA  directional  instability  even¬ 
tually  required  full  opposing  rudder  deflection,  which  was  immediately 
followed  by  departure. 


Overall  Assessment 


The  augmented  flight  control  system  designed  to  enhance  character¬ 
istics  in  normal  flight  improved  pre-departure  control  precision, 
reduced  departure  warning,  and  increased  the  tendency  to  flat  spin.  The 
roll  rate  CAS  and  the  SRI  together  improved  lateral-directional  control, 
which  allowed  higher  AOA  to  be  reached  before  sideslip  became  signifi¬ 
cant.  The  major  portion  of  this  benefit  derived  from  the  SRI.  For 
Configuration  C2  the  roll  rate  damper  suppressed  the  «£g&  contribution 
to  dynamic  stability  and  thus  increased  the  nose-slice  departure  and 
flat  spin  tendencies.  Thus,  this  augmentation  function  actually 
degraded  high  AOA  characteristics. 

C.  SUMMARY 

The  intended  high  AOA  departure  characteristics  were  obtained  in  the 
four  bare-alrframe  and  two  augmented  configurations.  The  results  demon¬ 
strate  the  strong  Influence  changes  in  the  six  key  aerodynamic  coeffi¬ 
cients  (£a>  if  g .  <Xp,  Nq,  Ng,  and  Mg)  can  have  on  departure  warning, 
severity,  and  susceptibility.  The  simulation  plus  literal  approximate 
factors  of  Section  111  show: 

•  The  most  violent  departure  characteristics  and 
flat  spin  tendem.  Lea  to  be  associated  with  large 
negative  Nq,  Ng,  and  Mg. 

•  The  least  violent  departure  characteristics  and 
greatest  spin  resistance  to  be  associated  with 
small  to  moderate  negative  and  Ng  and  positive 
Mg. 
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•  The  tendency  tc  Inconsistent  departure/spin  char¬ 
acteristics  (possibly  several  modes)  to  he  asso¬ 
ciated  with  positive  Cng^yn  brought  about  by 
large  ..negative  Co  Q  and  continued  large  adverse 
y.w  <»J.  LCDP ) 4  ls 

However,  the  observed  high  AOA  stall,  departure,  and  spin  dynamic  char¬ 
acteristics  also  depended  significantly  on  pilot  control  technique. 
This  is  consistent  with  actual  flight  experience  in  which  it  haB  been 
observed,  for  example,  that  spin  characteristics  and  recovery  techniques 
can  be  quite  different,  depending  upon  the  maneuver  from  which  the  spin 
occurred. 

The  final  high  AOA  departure/spin  rating  scale  of  Figure  73  was 
accepted  and  supported  by  the  pilots.  Results  of  this  first  application 
showed  consistent  reflection  of  both  the  airframe  dynamic  characteris¬ 
tics  and  any  sensitivity  to  difference  in  piloting  technique.  However, 
the  ratings  must  be  accompanied  by  detailed  qualitative  description  of 
the  characteristics  observed,  i.e.,  what  is  being  rated. 

In  this  vein,  it  might  be  noted  that  there  are  some  differences 
between  ratings  for  the  initial  familiarisation  and  tracking  runs 
reported  in  Subsection  B  and  for  the  final  tracking  runs  reported  In 
Subsection  A.  These  are  due  to  several  factors.  First,  Pilot  RC  was 
working  with  the  original  open-ended  scale  (Figure  72)  during  the  ini¬ 
tial  runs  and  the  closed-ended  scale  (Figure  73)  during  the  reruns; 
Pilot  JF  used  the  close-ended  scale  for  both  sets  of  runs.  Second, 
during  the  first  sequence  of  runs  Pilot  RC  was  unhurried  in  investi¬ 
gating  the  characteristics  of  each  configuration.  However,  for  the 
repeat  runs  he  was  under  some  pressure  to  complete  the  series  of  maneu¬ 
vers  and  configurations  within  time  constraints.  Thus,  the  rate  at 
which  he  approached  stall/departure  was  more  rapid  and  this  influenced 
his  assessment  of  warning  and  motion  severity,  as  demonstrated  in 
Table  16.  The  greatest  differences  are  for  the  augmented  flight  control 
configurations,  where  RC  saw  considerably  less  warning  and  more  severe 
departure  during  the  reruns.  Results  for  JF,  shown  for  comparison, 
reflect  quite  consistent  ratings  between  the  sets,  since  he  always  used 
a  rapid  pullup  technique. 
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TABLE  16 


COMPARISON  OF  PILOT  RATINGS, 
INITIAL  VS.  RERUNS 


CONFIGU¬ 

RATION 

RC 

JF 

WARNING 

MOTION  SEVERITY 

WARNING 

MOTION  SEVERITY 

INITIAL 

RERUN 

INITIAL 

RERUN 

INITIAL 

RERUN 

INITIAL 

RERUN 

A1 

2.8 

2.9 

2.2 

3.0 

4.8 

5.0 

4.6 

4.3 

B 

4.2 

3.9 

2.7 

3.1 

4.8 

4.6 

4.5 

4.8 

C1 

4.3 

3.6 

1.6 

2.0 

4.0 

4.6 

3.7 

3.9 

D 

4.0 

3.8 

3.5 

3.0 

3.5 

3.7 

1.5 

2.9 

A2 

2.8 

4.3 

1.8 

3.7 

4.8 

5.0 

4.6 

4.5 

C2 

3.2 

4.2 

1.5 

2.7 

4.8 

4.1 

3.8 

3.9 
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SECTIOI  VII 


DEPARTURE  SUSG8PTIBLITT  (XITERIOH 

Pitch  control  techniques  and  aggressiveness  had  obvious  Influences 

on  the  departure  warning  and  susceptibility  assessments  shown  In  the 

previous  section.  When  asked  about  specific  cues  and  personal  criteria 

which  would  prompt  the  pilot  to  acknowledge  departure  and  therefore 

to  Initiate  recovery,  vague  answers  were  obtained,  e.g>,  "persistent 

or  continued  uncommended  roll  or  yaw"  or  "suddon,  rapid  uncommanded 

motion  —  somewhere  around  20-30  deg/sec."  However,  detailed  scrutiny 

of  the  time  traces  failed  to  reveal  correlation  between  release  of  back 

•  ♦ 

stick  and  various  motion  quantities  (p,  r,  $,  ♦  »  6 »  etc.).  These  quan¬ 
tities  were  always  changing  rapidly  and  timing  became  a  big  factor. 
Deviations  of  fractions  of  a  second  in  initiating  recovery  (or  reading 
the  tracao)  produced  a  large  variation  in  the  motion  amplitudes  (e.g., 
see  Figures  57  and  58).  The  greatest  correlation  was  found  in  AOA 
Itself  —  possibly  because  Its  rate  of  change  was  the  slowest  of  all  and 
therefore  could  be  read  and  correlated  with  the  greatest  accuracy  with 
what  the  pilots  aaw. 

Since  the  piloting  task  required  closed-loop  attitude  control,,  and 
there  was  considerable  evidence  of  LCDP  or  b)£  being  a  key  parameter,, 
attention  was  turned  to  the  possibility  of  a  closed-loop  criterion.  A 
simple  analysis  of  theoretical  roll  loop  stability  was  found  to  agree 
quite  well  with  frequencies,  damping  ratios,  and  divergences  obtained 
from  the  time  traces.  Furthermore,  the  results  corroborated  the  differ¬ 
ences  In  pilot  techniques  and  observed  warning  reported  in  the  previous 
section. 

This  section  summarises  the  relationship  obtained  between  the 
closed-loop  anelyals  and  the  piloted  simulation  time  traces  and  presents 
a  closed-loop  criterion  which  is  consistent  with  and  explains  the  seem¬ 
ingly  diverse  departure  susceptibility  assessments  given  by  the  two  Air 
Force  test  pilots. 


155 


A.  CLOSED-LOOP  ANALYSIS 


The  lateral  stick  time  trace  of  Figure  57  (Configuration  A^)  is 
typical  of  most  runs  for  Pilot  RC.  The  phasing  between  6  excursions  and 
stick  deflections  Indicates  little,  if  any,  lead  being  generated  during 
the  PIO  at  a  -  16  deg  and  above.  Therefore,  the  simple  block  diagram  of 
Figure  78  was  chosen  to  model  the  unaugmented  airframe,  control  system, 
and  pilot.  A  Bode-root  locus  closed-loop  survey  plot  for  Configuration 
Aj  at  a  ■  16  deg  and  6-0  deg  is  presented  in  Figure  79.  Closed-loop 
gains  and  frequencies  obtained  from  strip  chart  recordings  of  simulation 
runs  for  Pilot  RC  are  shown  on  the  Bode  amplitude  plot  by  the  eeveral 
horlsontal  lines  (K^)  and  tic  marks  (-4-).  [For  example,  at  16  deg  AOA 
the  traces  of  Figure  57  show  a  pilot  gain  *  -  4.5  (13  dB)  and  fre¬ 
quency  w  *  1.8  rad/sec.]  The  frequency  of  the  undamped  oscillation 
obtained  from  the  strip  chart  is  in  remarkable  agreement  with  the  fre¬ 
quency  for  180  deg  phase  of  the  analytical  predictions,  validating  the 
assumptions  that  this  pilot  is  adopting  neither  leed  nor  lag  in  attempt¬ 
ing  to  control  roll  attitude.  Rather,  he  is  merely  reacting  with  a 
normal  time  delay  with  stick  proportional  and  opposite  to  bank  angle 
excursions. 


Figure  78.  Assumed  Pilot  Lateral  Loop  Closure;  $  +  Sgtk 
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The  measured  gains  range  from  about  10  to  20  dB.  This  ^riatlon. 
Identified  on  the  root  locus  plot  of  Figure  79,  indicates  but  small 
Influence  on  dutch  roll  divergence.  The  13  dB  gain  point  indicates  a 
mild  divergence  which  is  also  in  agreement  with  the  traces  of  Figure  57 
at  16  deg  AOA. 

As  AOA  is  increased  in  Figure  57,  the  lateral  stick  amplitude 
remains  fairly  constant  but  the  4>  excursions  increase.  Thus,  the  loop 
gain  (K^)  decreases  slightly,  the  oscillation  period  increases  to  about 
6  seconds  (w  *  l  rad/sec),  and  a  yaw  rate  divergence  develops.  Again,  a 
closed-loop  survey  plot  for  a  ■  23  deg  (Figure  80)  shows  excellent 
agreement,  l.e.,  uij  *  1  rad/sec  and  a  firat'-order  spiral  divergence 
of  T;  *  2  sec  which  results  from  the  spiral  mode  being  driven  toward  the 
RHP  xero  l/T^. 

Similar  trends  were  obvious  from  Pilot  RC  run  traces  for  moat  con¬ 
figurations.  He  appeared  to  be  following  the  neutral  stability  boundary 
of  the  oscillatory  mode  and  initiating  recovery  when  the  first-order 
divergence  mode  exceeded  soma  as  yet  unidentified  divergence  rate. 

Therefore,  using  the  same  pilot  model,  the  maximum  boundary  for 

stable  modes  was  calculated  and  plotted  as  shown  in  Figure  81.  This 

shown  two  boundaries.  The  upper  is  the  gain  not  to  be  exceeded  at  8  «  0 
and  the  lower  is  the  reduced  limit  gain  for  8  "  5.5  deg.  To  the  right 
(higher  AOA)  of  either  boundary  a  lateral  or  directional  instability 
exists  and  if  AOA  Is  increased  K^,  must  be  decreased  to  maintain  sta¬ 

bility. 

The  roll  loop  gains  employed  by  both  pilots  were  determined  from  the 
time  traces  and  located  on  similar  stability  boundary  plots  drawn  for 
each  vehicle  configuration.  Figure  82  is  the  plot  for  Configuration  k-L 
again,  but  now  showing  gain  bands  for.  each  pilot.  The  band  for  RC  shows 
that  he  tended  to  keep  his  gain  within  10  dB  above  the  maximum  stable 
boundary  and  to  decrease  his  gain  as  AOA  was  increased.  On  the  other 
hand,  JF  started  with  a  lower  gain  and  tended  to  keep  it  constant  during 
the  pull-up. 
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Figure  81.  Maximum  Pilot  Gain  for  Stable  $  ®atk  Closure 

Configuration  Ai 


Figure  82.  Roll  Loop  Closure  Gains  Employed  with 
Configuration 


Figure  S3  Is  Che  gain  plot  foe  Configuration  A2  •  Again  JF  used  a 
lower,  constant  gain  while  RC  tended  to  follow  the  stability  boundary 
until  the  sudden  drop.  Figure  84  is  the  plot  for  Configuration  C^, 
showing  results  similar  to  those  for  Configuration  A^.  This  configura¬ 
tion  is  somewhat  unusual  in  that  a  second  stable  region  exists  for  0^0 
at  very  high  AOA.  Figure  85  for  Configuration  C2  is  similar  to  Configu¬ 
ration  A2  except  there  is  no  dropoff  in  stability  boundary  at  high 
AOA.  The  results  for  Configuration  D  (Figure  86)  are  consistent:  Pilot 
RC  followed  the  stability  boundary  and  hence  should  detect  departure 
onset;  Pilot  JF  used  a  lower,  constant  gain  which  suddenly  crosses  the 
boundaries  as  AOA  is  increased.  He  therefore  obtained  little  warning. 
The  plot  for  Configuration  B  is  presented  in  Figure  87.  This  vehicle  is 
so  well  damped  in  roll  that  it  was  difficult  to  obtain  the  necessary 
data  from  time  traces  except  for  one  run  by  RC. 

At  the  bottom  of  Figures  82-87  negative  values  of  the  open-loop  roll 
numerator  root  1/T^  are  plotted  versus  a  (at  8  ■  0).  Figure  85  does 
not  show  this  extra  plot  because  u>^  is  positive  to  a  >  35  deg  for  this 
configuration.  Also  Identified  on  each  plot  is  the  average  AOA  at  which 
recovery  was  initiated  by  Pilota  RC  (  <>  symbol)  and  JF  (O).  It  will  be 
noted  that: 

•  In  most  instances  a  decrease  in  the  K<j,  boundary 
is  related  to  1/T^  becoming  negative  (RHP  zero). 

•  In  all  but  one  case  recovery  is  initiated  by  both 
pilots  after  exceeding  the  K,  boundary  and  exper¬ 
iencing  negative  1  /T,j, .  (the9  only  exception  is 
Configuration  C,  for  which  the  numerator  roots 
always  lie  in  the  LHP  and  there  is  no  dropoff  in 

y- 

•  Pilot  RC  always  employs  higher  roll  loop  gain, 
always  follows  the  Ka  boundary,  and  generally 
initiates  recovery  at9  lower  AOA.  However,  for 
Configurations  and  C2  he  did  explore  higher 
AOAs  before  recovery. 

•  Pilot  JF  tended  to  adopt  a  lower  roll  loop  gain 
and  kept  it  constant  during  his  more  rapid 
(aggressive)  pull-up.  He  generally  saw  little  or 
no  warning  and  penetrated  further  into  instabili¬ 
ty  regions  before  initiating  recovery. 
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Figure  83.  Roll  Loop  Closure  Gains  Employed  with 
Configuration  A2 
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Figure  84,  Roll  Loop  Closure  Cains  Employed  with 
Configuration  Ci 
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Figure  86.  Roll  Loop  Closure  Cains  Employed  with 
Configuration  D 
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Figure  87.  Roll  Loop  Closure  Gain  Employed 
with  Configuration  B 


The  foregoing  observations  appear  to  form  the  basis  for  each  pilot's 
definition  of  departure  susceptibility  (cf.  Section  VI): 

RC:  Departure  resistant  (R)  if  warning  is  clear  and 
consistent 

Departure  susceptible  (S)  if  warning  is  unclear, 
masked,  or  inconsistent 

JF:  Departure  resistant  (Fi)  ii:  always  recovered 

Departure  susceptible  (S)  if  sometimes  recovered 

Departure  extremely  susceptible  (ES)  if  never 
recovered 

Thus  each  pilot  is  viewing  different  aspects  of  the  departure.  Fortu¬ 
nately,  this  difference  lerves  to  increase  the  information  potential  of 
the  simulation  program.  This  is  probably  typical  of  flight  test,  and 
indicative  of  why  such  stall/spin  programs  should  be  flown  by  more  than 
one  pilot. 
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B.  POSSIBLE  DCPAETOU  SUSCKPTIIILITT  CUTRKIA 


A  comparison  between  departure/spin  susceptibility  predicted  by  the 
Weissman  criterion  end  the  assessments  provided  by  the  two  pilots  Is 
shown  Is  Table  17.  Assessments  substantially  in  agreement  with  predic¬ 
tion  are  shown  in  boxes.  Obviously  the  aggressive  pilot  observed  the 
worst  possible  characteristics  of  each  unaugmented  configuration  as 
predicted  by  the  criterion.  The  less  aggressive  pilot  experienced 
something  quite  different.  As  noted  previously,  the  criterion  does  not 
lend  Itself  to  prediction  of  the  frequency-dependent  augmented  airframe 
charcterletlcs,  but  these  configurations  were  expected  to  be  less  sus¬ 
ceptible  to  departure.  This  influence  was  observed  only  by  the  less 
aggressive  pilot. 


TABLE  17 

COMPARISON  OF  PREDICTED  AMD  ACTUAL  DEPARTURE 
SUSCEPTIBILITY  ASSESSMENTS 


DEPARTURE  SUSCEPTIBILITY 

CONFIGU¬ 

RATION 

PREDICTED 

RC 

JF 

A1 

ES 

R 

ES 

B 

ES 

S 

ES 

Cl 

S 

0 

S-ES 

D 

R 

0 

0 

a2 

NA 

R 

ES 

C2 

NA 

R 

S 

R  -  Resistant  ES  *  Extremely  susceptible 

8  ■  Susceptible  NA  ■  Not  applicable 
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One  Important  difference  between  prediction  and  simulation  was  the 
nature  of  departure.  Figure  63  predicta  predominantly  rolling  departure 
with  no  indication  of  yaw  departure.  Our  Configurations  A,  B,  and  D 
exhibited  initial  yaw  excursions  sometimes  followed  by  roll.  Configura¬ 
tion  C  exhibited  two  or  even  three  different  departure  modes  which  were 
dependent  upon  control  application  at  onset  of  departure.  However,  a 
rolling  type  motion  did  predominate.  Thus,  the  WelBsman  criterion  was 
not  as  useful  in  predicting  the  nature  of  departure  observed  by  the 
pilot. 

The  closed-loop  analysis  to  Identify  causal  factors  behind  the 
widely  differing  pilot  ratings  for  departure/spin  susceptibility  pro¬ 
duced  a  strong  relationship  between  pilot  ratings  for  both  pilots  and 
penetration  of  the  roll  numerator  into  the  RHP.  Figure  88  shows  the 
value  of  the  real  part  of  the  numerator  root  at  the  instant  the  pilot 
decided  he  had  departed,  or  was  about  to  depart,  and  initiated  recovery. 
These  values  are  plotted  against  the  AOA  at  which  recovery  was  started. 
(Note  this  is  not  the  usual  root  locus  ju>  axis.)  The  points  represent 
all  six  vehicle  configurations  as  evaluated  by  both  pilots. 

The  division  between  departure-resistant  (R)  and  susceptible  (S) 
ratings  is  seen  to  lie  at  roughly  -0.5  rad/sec.  This  corresponds  to  a 
time  to  double  amplitude  of  approximately  1.4  sec.  Zeros  which  lie  to 
the  left  of  this  line  apparently  limit  the  first-order  divergence  to  a 
rate  slow  enough  for  pilots  to  respond  and  recover.  Zeros  to  the  right 
of  the  line  apparently  allow  divergence  rates  so  fast  that  the  pilots 
cannot  prevent  departure.  Again  this  is  consistent  with  the  previously 
noted  definition  of  departure  whereby  both  pilots  indicated  a  threshold 
on  rate  of  motion;  however,  the  pilots  were  vague  as  to  the  value  (e.g., 
"maybe  20  or  30  deg/sec").  One  data  point  in  Figure  88  violates  the 
boundary.  This  is  the  augmented  Configuration  C2>  In  this  case  the  SRI 
eliminated  adverse  yaw  and  thus  made  the  vehicle  more  departure  resis¬ 
tant  (as  viewed  by  RC).  However,  it  could  be  departed;  and  then  the 
augmentation  produced  pro-spin  control.  Pilot  JF  rated  spin  suscepti¬ 
bility  (and  recovery)  and  apparently  rated  this  configuration  accord¬ 
ingly. 
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Figure  88.  Departure  Susceptibility  Rating  Versus 
Lateral  Closed-Loop  Divergence  Potential 


The  implication  Is  that  If  the  combined  aerodynamics  and  flight 

control  system  design  Is  such  that  1/T^  never  exceeds  -0.5  throughout 

the  achievable  AOA  range,  the  airplane  will  be  departure-resistant.  It 

should  be  noted  that  this  criterion  places  no  restriction  on  open-loop 

stability.  For  example,  Cn.  can  be  nogatlve  and,  in  fact,  is  nega- 

pdyn 

tlve  for  Configuration  D  (see  Figure  52)  which  is  rated  departure- 
resistant  (see  Table  13)  by  both  pilots. 

Since  u) ^  is  the  dimensional  form  of!  LCDP  and  since,  in  general, 
|1/T^|2  A  | ,  then  one  can  relate  the  above  1/T^  boundary  to  an 
equivalent  LCDP.  For  the  flight  conditlona,  Inertias,  etc.,  employed 
in  this  simulation,  1/T<^  of  -0.5  corresponds  to  LCDP  of  -0.001.  This 
coincides  with  (talesman's  boundary  between  Regions  A  end  B  for  positive 
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Cna  ,  see  Figure  89;  however,  it  is  a  little  more  conservative  at 
pdyn 

negative  Cn^  .  Thus,  the  results  of  our  simulation  are  compatible 
with  and  support  Weissman"s  empirically  derived  LCDP  boundary.  The  key 
difference  in  the  criterion  is  that  1/T^  is  not  restricted  to  airframe 
lateral-directional  static  coefficients  but  can  be  applied  to  the  com¬ 
pletely  coupled  6  DOF  airframe  with  a  full  complement  of  augmentation, 
stlck-to-rudder  crossfeed,  etc.,  throughout  the  aircraft  development 
cycle. 


In  addition,  results  of  our  simulation  tend  to  indicate  that  another 
boundary  might  be  appropriate  in  the  upper  left  quadrant  of  Figure  89  to 
create  two  additional  regions,  E  and  F,  as  suggested  in  Reference  40. 
Region  E  would  be  classified  as  mild  directional  divergence  and  moderate 
spin  tendency.  Region  F  might  be  classified  bb  severe  directional 
divergence  and  strong  spin  tendency.  Note  from  Figure  63  that  our  air¬ 
craft  Configuration  D,  which  had  mild  departure  characteristics  but 
little  or  no  spin  tendency,  would  extend  into  Region  E. 


LCDP 


Figure  89.  Possible  Modifications  to  Weissman  Criterion 
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Additional  support  for  the  limit  on  1/T^  can  be  drawn  from  the 
current  MIL-F-8785C  spiral  divergence  limit*  The  Level  3,*  Category  A 
and  C  Flight  Phases  requirement  is  that  time  to  double  amplitude  be  not 
lees  than  4  sec*  This  limit  came,  in  part,  from  the  Reference  41  inves- 
tigation  in  which  aircraft  configurations  having  various  spiral  diver¬ 
gence  races  were  rated  by  twelve  pilots  in  a  landing  approach  task* 
Results  of  that  investigation,  shown  in  Figure  90,  indicate  that  diver¬ 
gences  having  a  time  to  double  amplitude  of  less  than  2.4  sec  were 
considered  Intolerable.  The  time  to  double  amplitude  (T2)  can  be 
converted  to  an  equivalent  first-order  time  constant  (T)  by  the  rela¬ 
tionship 


T2  *  0.74T 

A  plot  of  pilot  rating  versus  spiral-mode  inverse  time  constant  (1/T)  is 
presented  in  Figure  91.  The  current  -3785C  Level  3  spiral  limit  is 
shown  along  with  the  proposed  1/T^  limit. 

The  two  limits  are  in  good  agreement  when  it  is  remembered  that 
1/Ta  represents  a  pole  (modal  response)  while  1/T^j  represents  a  zero. 
Under  closed-loop  roll  control  the  migration  of  an  airframe  pole  to¬ 
ward  l/T^  would  be  dependent  upon  pilot  gain;  infinite  loop  gain  would 
be  required  to  obtain  a  closed-loop  pole  exactly  at  the  zero  location. 
A  more  realistic  situation  is  shown  in  Figure  92,  which  is  a  Bode-siggle 
plot  for  Configuration  Aj  at  «  ■  21  deg,  8  •  0  deg.  The  dashed  line 
terminating  at  l/T^  ■  -0.6  is  the  locus  of  the  closed-loop  first-order 
divergence  root  for  increasing  gain,  .  Two  typical  gain  lines  for 
Pilot  RC  extracted  from  simulation  traces  at  a  -  21  ±  1  deg  are  shown. 
These  closures  result  in  closed-loop  roots  at  -0.2  and  -0.33  rad/sec. 
Referring  back  to  Figure  91,  it  may  be  observed  that  these  straddle  the 


A 

Level  3:  “Flying  qualities  such  that  the  airplane  can  be  con¬ 
trolled  safely,  but  pilot  workload  is  excessive  or  mission  effectiveness 
is  Inadequate,  or  both.  Category  A  Flight  Phaaes  can  be  terminated 
safely,  and  Category  B  and  C  Flight  Phases  can  be  completed •“ 
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Figure  90.  Variation  of  Pilot  RatingB  with  Time  to  Double  Amplitude 
of  Aperiodic  Motion;  Landing  Approach  Condition  (from  Ref.  40) 
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Figure  91.  Variation  of  Pilot  RatingB  with  Equivalent 
Aperiodic  Inver ne  Time  Constants 
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1/Tg  boundary  and  are  consistent  with  the  Reference  40  division  between 
tolerable  and  intolerable  handlings  Again,  our  Pilot  RC  considered 
Configuration  to  be  departure-resistant  (R). 

The  above  discussion  tends  to  show  that  the  1/T^  limit  selected 
(-0.5)  may  not  be  conservative,  since  a  very  tight  roll  loop  closure 
could  result  In  closed-loop  roots  in  the  intolerable  range  of  Figure 
91.  Further  investigation  will  be  required  to  determine  if  the  limit 
should  be  reduced. 

C.  SUMMARY 


Analysis  and  piloted  simulation  have  Bhown  that: 


The  pilot's  perception  of  departure  suscepti¬ 
bility  was  found  to  be  correlated  with  movement 
of  one  root  of  the  roll  numerator  for  lateral 
stick  control  into  the  RHP  of  the  root  locus, 
i.e.,  a  nonmlnlmum-phaae  zero.  If  pilots  could 
fly  the  aircraft  to  such  AOAs,  they  rated  high 
AOA  configurations  which  produced  a  zero,  1/Tx  , 
more  negative  than  -0.5  departure-susceptible. 
If  this  boundary  was  not  or  could  not  be  ex¬ 
ceeded,  the  aircraft  was  considered  departure- 
resistant.  This  rating  is  not  a  function  of  the 
sign  or  magnitude  of  the  dynamic  stability  param- 

Cter  Cn0dyn* 

A  value  of  1/Tx.  ■  -0.5  corresponds  for  the  air¬ 
frame  tested  to1  an  effective  LCDP  of  -0.001  and 
thus  is  consistent  with  and  supports  the  empiri¬ 
cally  derived  LCDP  departure  boundary  developed 
by  Weissman.  However,  results  of  the  simulation 
were  not  in  agreement  with  the  types  of  departure 
predicted  by  Weissman  in  that  negative  Cngd  re" 
glons  produced  yaw  or  nose-slice  type  departures. 


•  A  value  of  1  /T a .  -  -0.5  also  is  compatible  with 
the  current  MIIHF-8785C  Level  3  limit  for  spiral 
divergence  iu  Category  A  and  C  Flight  Phases. 

It  is  suggested  that  as  the  aircraft  deBign/development  cycle  pro¬ 
gresses  past  the  static  wind  tunnel  phasr  the  Weissman  criterion  for 
departure/spin  susceptibility  be  replaced  by  a  negative  l/T^  limit. 
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SECTION  VIII 


CONCLUSIONS  AND  RKCQMfENDATIONS 


The  goals  of  this  research  program  were  to: 

•  Identify  key  design  parameters  that  limit  high- 
angle-of-attack  maneuverability  for  contemporary 
high-performance  attack  and  fighter-type  air¬ 
craft. 

•  Postulate  fundamental  aerodynamic  and  control 
system  design  methodologies  that  will  alleviate 
the  limiting  conditions. 

•  Formulate  handling  qualities  requirements  for 
high-angle-of-attack  maneuvering  flight  to  be 
incorporated  in  MIL-F-8785C,  "Flying  Qualities  of 
Piloted  Airplanes." 

It  should  be  borne  In  mind  that  the  results  are  based  on  variations  on  a 
single  nonlinear  aerodynamic  model  representing  a  region  dominated  by 
phenomena  that  ore  highly  configuration-dependent.  Therefore,  the  gen¬ 
eralizations  drawn  are  In  need  of  further  substantiation. 

Km  KEY  MAHEUVER-LIMTTIW;  FACTORS 

The  most  crucial  factors  pertain  to  limiting  the  safe  flight  enve¬ 
lope,  i.e.,  departure  from  controlled  flight.  By  varying  static  aerody¬ 
namic  coupling  and  cross-coupling  coefficients,  the  nature  of  high  ADA 
characteristics  and  susceptibility  to  departure  were  changed  to  reflect 
nose  slice,  rolling,  and  pitch-up  departures. 

Nose  slice  (yaw)  departures  predominated  for  configurations  having 
Cngdyn  large,  negative.  Analysis  and  prior  simulation  have  shown  that 
nose  slice  can  be  aggravated  once  sideslip  startB  by  the  static  aerody¬ 
namic  coefficients  Cj  and  Cn  .  These  aerodynamic  coupling  terms  can 

cause  a  cross-coupled  RHP  zero  to  appear  in  the  pitch  numerator,  N?  . 

stab 

Any  simultaneous  attempt  to  control  pitch  attitude  then  produces  or  fur¬ 
ther  aggravates  directional  divergence.  The  large  negative  Cna  con- 

Bdyn 

figurations  were  prone  to  flat  spins. 
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Rolling  departures  predominated  when  a  large  C component  made 
^n3dyn  positive.  These  configurations  also  were  prone  to  wing  rock  and 
lateral  PIO.  They  exhibited  at  least  two  different  PSG  and  spin  modes, 
and  were  sensitive  to  maneuver  and  pilot  activity  at  departure  onset. 


Pitch  due  to  sideslip,  Cmg,  of  either  positive  or  negative  sign 
can  adversely  influence  departure  and  handling  characteristics.  Posi¬ 
tive  Cmp  can  produce  pitch-up  which  tends  to  mask  a  "g-break"  stall 
warning  cue  and  oppose  pilot  attempts  to  prevent  stall/departure  via 
lorward  stick.  Negative  Cm^  tends  "o  augment  "g-break"  stall  warning 
but  reinforces  pilot-commanded  recovery  pitching  moment  to  the  point 
that  pitch-roll-yaw  coupling  may  produce  violent  PSG.  The  presence  of 
appreciable  Cra^  can  be  detected  from  flight  traces  of  motion  about  a 
nominal  zero  Bideslip  via  a  longitudinal  oscillation  at  twice  the  dutch 
roll  frequency. 


Unstable  open-loop  characteristic  modes  (poles)  alone  were  not  found 
to  limit  the  safe  flight  envelope.  If  the  pertinent  control  numerator 
zeros  were  favorably  located  (LHP),  the  pilot  could  prevent  divergence. 
If  the  zeros  were  unfavorably  located  (RHP)  divergence  could  not  be  pre¬ 
vented;  in  this  case  the  presence  of  open-loop  instability  increased 
departure  severity.  Thus,  pole-zero  locations  had  influence  on  depar¬ 
ture  susceptibility  and  severity.  In  turn,  pole-zero  locations  were 
shifted  by  static  aerodynamic  cross-coupling  derivatives  (oCq,  N^,  Mg) 
when  the  magnitudes  of  these  coefficients  approached  those  of  the  common 
static  stability  derivatives  <4;,  Ng,  Mq): 

•  X  'ai  N^,  and  Mg,  in  combination,  influence  the 
high-frequency  modes  of  both  the  lateral  and 
longitudinal  denominators. 
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B.  ADDITIONAL  FINDINGS  AND  OONOUJSIONS 

One  major  contribution  has  been  identification  and  validation  that 
pilot  perception  of  lateral-directional  departure  susceptibility  is 
related  to  one  zero  of  the  numerator  becoming  negative.  Root 

magnitudes  more  negative  than  -0.5  rad/sec  were  consistently  rated  as 
departure-susceptible,  while  those  less  negative  (or  positive)  are  rated 
as  departure-resistant.  This  criterion  reflects  a  closed-loop  diver¬ 
gence  rate  limit  related  to  the  pilot's  threshold  for  uncommanded  motion 
or  ability  to  cope.  As  such  it  is  a  pilot-centered  criterion  which 
should  be  applicable  for  any  flight  situation,  although  it  has  been 
Identified  in  a  low-Mach-number,  fixed-base  simulation.  It  is  consis¬ 
tent  with  the  empirically  established  airframe-alone  departure/spin 
criterion  boundaries  of  Waisaman  and  extends  applicability  of  that 
criterion  to  highly  augmented  airframe  cases.  It  is  also  consistent 
with  previous  in-flight  simulation  of  maximum  controllable  aperiodic 
divergence  rates.  Finally,  it  serves  as  both  a  design  guide  and  a 
flying  quality  specification  item. 

The  "unstable"  (RHP)  zero  generally  results  from  negative  (adverse) 
Cng  and  Cn^g^.  It  therefore  can  be  alleviated  by  aircraft  configura¬ 
tion  and  flight  control  system  modifications  which  principally  alter 
those  coefficients.  As  noted  above,  the  zero  also  can  be  highly  influ¬ 
enced  by  the  static  aerodynamic  cross -coupling  coefficients  and  Cng. 

A  second  major  contribution  has  been  the  development  and  partial 
validation  of  a  flying  quality  rating  scale  applicable  to  the  stall/ 
departure/recovery  flight  regime.  It  is  believed  that  the  scale  has 
potential  for  defining  the  aircraft  permissible  flight  envelope  in  much 
the  same  manner  as  the  Cooper-Harper  scale  is  now  used  in  defining 
flying  quality  bounds. 

The  simulation  demonstrated  the  Importance  of  loss-of-control 
warning  in  pilot  assessment  of  high  AOA  flying  qualities  and  departure 
susceptibility.  A  gradual  and  consistent  degradation  in  stability  and 
controllability  was  found  to  provide  the  best  warning;  however,  the 
realization  of  such  warning  can  be  highly  dependent  upon  pitch  control 
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power.  Oscillations  such  as  wing  rock  or  lateral  PIO  cannot  be  counted 
on  for  departure  warning  in  air  combat  if  the  aircraft  can  be  pulled 
through  the  '•warning"  region  before  the  oscillation  becomes  apparent  to 
the  pilot. 

Low-frequency  dutch  roll  or  lateral  phugoid  oscillation  was  found  to 
mask  departure  onset.  Roll  about  the  velocity  vector  at  high  AOA  pro¬ 
duces  significant  "nose"  yaw.  The  pilot  has  considerable  difficulty 
differentiating  between  nose  slice  onset  and  low-frequency  roll  oscilla¬ 
tion  about  the  velocity  vector.  High  roll  rate  damping  (natural  or 
augmented)  was  also  found  to  reduce  high  AOA  departure  warning  and 
Increase  flat  spin  tendency. 

A  flight  control  or  augmentation  system  mechanization  which  improves 
and  extends  the  aircraft  controllability  range  by  removing  natural 
warning  can  actually  be  a  detriment  if  the  aircraft  can  still  reach  some 
departure  AOA.  To  be  of  real  benefit  such  a  flight  control  system  must 
also  prevent  departure. 

A  stick-to-rudder  interconnect  (SRI)  can  be  beneficial  in  reducing 
adverse  "aileron"  yaw,  thereby  providing  more  favorable  LCDP  or  1/T^ 
values.  However,  if  static  aerodynamic  cross-coupling  is  strong,  the 
analysis  of  Configuration  A2  has  shown  that  even  small  sideslip  can 
shift  the  numerator  roots  in  a  manner  to  negate  the  SRI  contribution. 
These  results  also  imply  that  a  finite-authority  [3-limiter  may  be  of 
little  benefit  as  a  departure  preventer  if  one  or  more  RHP  zeros  occur 
in  control  loops  the  pilot  is  closing.  That  is,  pilot  control  can  drive 
the  vehicle  into  a  divergence  which  may  exceed  the  (S-liraiter's  capa¬ 
bility. 

Reaults  of  this  study  demonstrate  that  linear  or  quasi-llnear 
frozen-point  analysis  is  applicable  to  high  AOA  situations,  but  consid¬ 
erable  caution  must  be  exercised  in  the  interpretation  of  results  where 
aerodynamics  are  strong  functions  of  angles  of  attack  and  sideslip.  Any 
analysis  pertaining  to  such  regions  must  be  based  on  a  6  DOF  model  with 
aerodynamic  and  kinematic  cross-coupling  terms,  must  include  non-zero- 
sldealip  trim  points,  must  view  transfer-function  numerator  as  well  as 
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characteristic-equation  parameters,  and  must  be  supported  by  nonlinear 
simulations. 

In  particular,  the  results  demonstrated  the  power  of  vector  polygon 
and  closed-loop  root  extraction  analytic  techniques  for  identifying 
cause-effect  relationships  in  highly  coupled  dynamic  systems.  The  tech¬ 
niques  facilitated  development  of  transfer-function  literal  approximate 
factors  for  the  7-4  in  high  AOA,  non-zero-sideslip  flight  conditions. 

C.  PROPOSED  MIL-F-8785C  MODIFICATIONS 

Based  upon  results  of  this  program  it  is  recommended  that  modifica¬ 
tions  to  MIL-F-8785C  requirements  be  considered  in  three  areas. 

I.  Departure  from  Controlled  Flight 

The  current  paragraph  (3. 4. 2. 2.1}  simply  requires  that  "all  Classes 
of  aircraft  be  extremely  resistant  to  departure  from  controlled  flight, 
post-stall  gyrations  and  spina.... The  aircraft  shall  exhibit  no  uncom¬ 
manded  motion  which  cannot  be  arrested  promptly  by  simple  application  of 
pilot  control." 

The  "uncommanded  motio  i"  requirement  should  be  strengthened  to 
require  that,  for  roll  control  input  within  the  service  flight  envelope 
(SFE)  the  following  stability  axis  parameter  values  shall  be  no  more 
negative  than: 


Unaugmented  airframe >  LCDP 


>  -0.01 


Augmented  airframe:  1/T^  >  -0.5 

An  alternate  requirement  might  be  that  the  aircraft  should  exhibit 
no  aperiodic  uncommanded  motion  which  exceeds  20  deg/sec  and  cannot  be 
arrested  promptly  by  simple  application  of  pilot  control.  [The  value 
selected  here  is  baaed  on  a  rough  average  of  the  simulation  pilots' 
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commentary  as  to  their  definitions  of  departure,  the  spiral  divergence 
mode  limit  of  Reference  40,  and  the  1/T^  limit  above.] 

Since  the  SFE  specifically  excludes  Btall  and  departure,  the  above 
criteria  would  automatically  cover  high  AOA,  pre-stall  flight.  Other 
potential  locations  for  such  wording  are  paragraphs  3. 4. 2. 1.1,  Stall 
Approach,  and  3. 4. 2. 1.2,  Stall  Characteristics. 

2.  Stall  Definition 

The  current  paragraph  (3.4.2. 1)  allows  definition  of  stall  to  be 
based  on  abrupt  uncontrollable  pitching,  rolling,  or  yawing;  or 

intolerable  buffet.  The  results  of  our  piloted  simulation  indicate  that 
any  abrupt  aperiodic  rolling  or  yawing  motion  which  octurB  without  being 
preceded  by  noticeable  "g-break"  is  considered  to  be  a  departure,  not 
a  stall,  and  results  in  severely  downgraded  flying  qualities.  Thus, 
abrupt  roll  or  yaw  motion  should  be  deleted  as  an  allowable  definition 
of  stall  and  should  not  occur  prior  to  stall.  If  such  characteristics 
cannot  be  achieved  with  the  airframe  alone,  then  the  flight  control 
system  should  prevent  reaching  the  AGA  at  which  the  abrupt  rolling  or 
yawing  motion  is  obtained. 

3.  Departure  Warning 

The  current  paragraphs  (3. 4. 1.1,  3. 4. 2. 1.1.1,  and  3. 4. 2. 1.1. 2) 

covering  warning  pertain  mainly  to  stall.  But,  since  abrupt  uncommanded 
roll  or  yaw  is  also  considered  to  define  stall,  the  Paras.  3. 4. 2. 1.1.1 
and  3. 4. 2. 1.1. 2  requirements,  in  effect,  define  departure  warning. 
These  requirements  are  based  upon  fixed  margins  of  VHtajj  or  ^L8tall  ant* 
do  not  take  into  account  the  severity  of  any  post-stall/departure  mo¬ 
tion,  uncontrolled  altitude  loss,  or  mission  phase. 

Results  of  this  simulation  indicated  that  stall/departure  warning 
margin  should  be  related  to  the  severity  of  any  uncommanded  motion  or 
PSG.  Time  and  altitude  loss  prior  to  recovery  were  heavily  dependent 
upon  the  initial  departure  severity;  pilot  commentary  indicated  the 


178 


overall  departure  ratings  obtained  were  heavily  influenced  by  alti¬ 
tude  loss  and  mission  phase.  Quite  obviously,  consequences  of  stall/ 
departure  in  a  ground  attack  phase  can  be  more  catastrophic  than  in 
high-altitude  air  combat.  Unfortunately,  insuf ficieR';  data  ware  ob¬ 
tained  to  establish  a  specific  departure  warning  cilt  -s/ton. 

4.  Compliance  Demonstration 

The  current  specification  (Para.  4.1)  allows  compliance  with  all 
requirements  of  Section  3  to  be  demonstrated  through  analysis.  Although 
compliance  with  many  of  the  requirements  will  be  demonstrated  by  simula¬ 
tion,  test,  or  both,  this  analysis  and  simulation  program  has  demon¬ 
strated  that  linearised,  frosen-point  analytic  results  may  be'  midleading 
or  extemely  difficult  to  interpret  for  flight  regions  where  the  airframe 
may  exhibit  large  static  aerodynamic  cross-coupling  moments.  Therefore, 
it  is  recommended  that  the  requirement  be  changed  so  that  compliance 
with  paragraphs  pertaining  to  atall/departure  must  be  demonstrated  via 
simulation  and/or  flight. 
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